arxiv: v1 [math-ph] 14 Feb 2014

Transcription

1 AN ELECTROSTATICS PROBLEM ON THE SPHERE ARISING FROM A NEARBY POINT CHARGE JOHANN S. BRAUCHART, PETER D. DRAGNEV, EDWARD B. SAFF arxiv:4.3367v [math-ph] 4 Feb 4 Abstract. For a positively charged insulated d-dimensional sphere we investigate how the distribution of this charge is affected by proximity to a nearby positive or negative point charge when the system is governed by a Riesz s-potential /r s, s >, where r denotes Euclidean distance between point charges. Of particular interest are those distances from the point charge to the sphere for which the equilibrium charge distribution is no longer supported on the whole of the sphere i.e. spherical caps of negative charge appear). Arising from this problem attributed to A. A. Gonchar are sequences of polynomials of a complex variable that have some fascinating properties regarding their zeros.. Introduction For the insulated unit sphere S d in R d+ of total charge + on which particles interact according to the Riesz-s potential /r s, s >, where r is the Euclidean distance between two particles, the equilibrium distribution of charge is uniform; that is, given by normalized surface area measure σ d on S d. However, in the presence of an external field due to a nearby point charge the equilibrium distribution changes. A positive external field repels charge away from the portion of the sphere near the source and may even clear a spherical cap of charge, whereas a negative external field attracts charge nearer to the source, thus thinning out a region on S d opposite to the direction of the source. In the Coulomb case s = and d = this is a well-studied problem in electrostatics cf., e.g., [4]). Here we deviate from this classical setting and show that new and sometimes surprising phenomena can be observed. The outline of the paper is as follows. In Section we introduce and discuss a problem Gonchar s problem) concerning the critical distance from a unit point charge to S d such that the support of the s-equilibrium measure on the sphere for Riesz s-potential /r s, s > ) is no longer all of the sphere when the point charge is at any closer distance. We shall make this more precise below. Date: August 8, 8. Mathematics Subject Classification. Primary 3C, 3B5; Secondary 8A, 3C5, 3B. Key words and phrases. Electrostatics problem; Golden ratio; Gonchar problem; Gonchar polynomial; Plastic Number; Riesz potential, Signed Equilibrium; Sphere. The research of this author was supported, in part, by an APART-Fellowship of the Austrian Academy of Sciences, a Grants-in-Aid program of ORESP at IPFW, and, by the Australian Research Council. The research of this author was supported, in part, by a Grants-in-Aid program of ORESP at IPFW and by a grant from the Simons Foundation no. 87. The research of this author was supported, in part, by the U. S. National Science Foundation under grant DMS-966 as well as by an Australian Research Council Discovery grant.

2 J. S. BRAUCHART, P. D. DRAGNEV, E. B. SAFF Our starting point is the solution to the signed equilibrium problem for a positive charge outside the sphere, which is intimately connected with the solution of this external field energy problem. Using similar methods, we are able to extend the results in [5] to external fields due to a positive/negative point charge inside/outside the sphere. We shall analyze Gonchar s problem for different magnitudes and for both positive and negative charges which, in turn, will provide a more comprehensive picture than presented in [5] and [7], and will reveal some interesting new phenomena. For example, in the case when d s is an even positive integer, a logarithmic term appears in the formula for the critical distance; and for a negative external field due to a source inside the sphere, a crucial issue is whether the Riesz kernel is strictly subharmonic d < s < d) or strictly superharmonic < s < d ). For the logarithmic potential, we are able to provide a complete answer for Gonchar s problem. In the classical harmonic) case s = d and more generally, when d s is an odd positive integer) the critical distance appears as a zero of a certain family of polynomials indexed by the dimension d. The classical case and the associated polynomials are discussed in some detail in [7]. By allowing a negative as well as positive charge q and allowing this charge to be inside or outside the sphere, we derive a family of polynomials for each of the four combinations: i) R < and q >, ii) R > and q >, iii) R < and q <, and iv) R > and q <. In Section 3, we investigate these four families of polynomials arising from Gonchar s problem. Figure 5 illustrates the qualitative patterns of their zeros in the classical case s = d and illustrates how these families complement each other. The last two displays of this figure are for s = d 3. Figure 6 contrasts the cases s = d 3 and s = d 5 for weak q = /), canonical q = ) and strong q = ) external fields. In [6] we discussed the Riesz external field problem due to a negative point charge above the South Pole of a positively charged unit sphere total charge ) and we derived the extremal and the signed equilibria on spherical caps for d < s < d. In Section 4, we provide the details of this derivation and consider also the limiting cases when s d. Section 5 contains the remaining proofs. In the Appendix we study the s-potential of uniform normalized) surface area measure on the d-sphere in more detail.. Potential Theoretic Setting for Gonchar s Problem Let A be a compact subset of S d. Consider the class MA) of unit positive Borel measures supported on A. The Riesz-s potential and Riesz-s energy of a measure µ MA) modeling a positive) charge distribution of total charge on A are defined as U s µ x) := x y s d µy), I s [µ] := x y s d µx) d µy). The Riesz-s energy and the s-capacity of A are given by { } W s A) := inf I s [µ] : µ MA), cap s A) = /W s A). It is well-known from potential-theory cf. Landkof [7]) that if A has positive s-capacity, then there always exists a unique measure µ A,s MA), which is called the s-equilibrium

3 AN ELECTROSTATICS PROBLEM ON THE SPHERE 3 measure on A, such that W s A) = I s [µ A,s ]. For example, σ d is the s-equilibrium measure on S d for each < s < d. A standard argument utilizing the uniqueness of the s-equilibrium measure µ A,s shows that it is the limit distribution in the weak-star sense) of a sequence of minimal s-energy N-point configurations on A minimizing the discrete s-energy E s x,..., x N ) := N N j= k= j k x j x k s, x,..., x N A, over all N-point systems on A. For the discrete s-energy problem we refer to [3].) Weighted Energy and External Fields. We are concerned with the Riesz external field generated by a positive or negative point charge of amount q located at a =, R) on the polar axis with R < or R >. Such a field is given by ) Qx) = Q R,q,s x) := q/ x a s, x R d+. The Riesz-s external field on a compact subset A S d with positive s-capacity is Q restricted to A. For simplicity, we use the same symbol.) The weighted s-energy I Q [µ] associated with such a continuous external field and its extremal value V Q are defined as I Q [µ] := I s [µ] + Qx) d µx), V Q A) := inf {I Q [µ] : µ MA)}. A measure µ Q MA) such that I Q [µ] = V Q A) is called an s-extremal or positive equilibrium) measure on A associated with Q. This measure is unique and it satisfies the Gauss variational inequalities cf. []) ) 3) where U µ Q s x) + Qx) F Q A) everywhere on A, U µ Q s x) + Qx) F Q A) everywhere on support suppµ Q ) of µ Q, F Q A) := V Q A) Qx) d µ Q x). In fact, once suppµ Q ) is known, the equilibrium measure µ Q can be recovered by solving the integral equation U µ s x) + Qx) = everywhere on A for positive measures µ supported on A. In the absence of an external field Q ) and when cap s A) > the measure µ Q coincides with µ A,s. Riesz external fields due to a positive charge on the sphere S d were instrumental in the derivation of separation results for minimum Riesz-s energy points on S d for s d, d) see []). In [5] we studied Riesz external fields due to a positive charge above S d which led to a discussion of a fascinating sequence of polynomials arising from answering Gonchar s problem for the harmonic case in [7]. The least separation of minimal energy configurations on S d subjected to an external field is investigated in [6]. In [5] we also developed a technique for finding the extremal measure associated with more general axis-supported Note that a continuous negative field can be made into a positive one by adding a fixed constant.

4 4 J. S. BRAUCHART, P. D. DRAGNEV, E. B. SAFF fields. For external fields in the most general setting we refer the reader to the work of Zoriĭ [3, 3, 33] also cf. []). Signed Equilibrium. The Gauss variational inequalities ) and 3) for A = S d imply that the weighted equilibrium potential is constant everywhere on the support of the measure µ Q on S d. In general, one can not expect that the support of µ Q is all of the sphere. A sufficiently strong external field large q > or small R > ) would thin out the charge distribution around the North Pole and even clear a spherical cap of charge. In this insulated sphere setting there is no negative charge that would be attracted to the North Pole.) By enforcing constant weighted potential everywhere on the sphere, in general, one has a signed measure as a solution. This corresponds to a grounded sphere.) In the classical Coulomb case d =, s = ) a standard electrostatic problem is to find the charge density signed measure) on a charged, insulated, conducting sphere in the presence of a point charge q off the sphere see [4, Ch. ]). This motivates the following definition cf. []). Definition. Given a compact subset A R p p 3) and an external field Q on A, we call a signed measure η Q = η A,Q,s supported on A and of total charge η Q A) = a signed s-equilibrium on A associated with Q if its weighted Riesz s-potential is constant on A; i.e., 4) U η Q s x) + Qx) = G A,Q,s everywhere on A. Physicists usually prefer neutral charge η Q A) =. However, for the applications here it is more convenient to have the normalization η Q A) =. It can be shown that if a signed equilibrium η Q on A exists, then it is unique see []). We remark that the determination of signed equilibria is a substantially easier problem than that of finding non-negative extremal measures. However, the solution to the former problem is useful in solving the latter problem. In [5] it is shown for q > and R > that the signed s-equilibrium η Q on S d associated with the Riesz external field ) is absolutely continuous with respect to the normalized surface area measure on S d. Using Imaginary inversion cf. Landkof [7]), the proof can be extended to hold for the class of fields considered here. We remark that in the Coulomb case s = and d = ) this result is well-known from elementary physics cf. [4, p. 6]). Theorem. Let < s < d. The signed s-equilibrium η Q on S d associated with the external field Q of ) is absolutely continuous with respect to the normalized surface area measure on S d ; that is, d η Q x) = η R,q,s x) d σ dx), and its density is given by 5) η R,q,sx) = + qu σ d s a) W s S d ) q R d s W s S d ) x a, x d s Sd. The charge q can be positive or negative and the distance of the charge to the sphere center satisfies R < or R >.) Moreover, the weighted s-potential of η Q on S d equals 6) G S d,q,s = W s S d ) + qu σ d s a). The case d =, s =, where the source is a point on the unit circle, was investigated in [5].

5 AN ELECTROSTATICS PROBLEM ON THE SPHERE 5 In the following we shall use the Pochhammer symbol a) :=, a) n := aa + ) a + n ), n, which can be expressed in terms of the Gamma function Γ by means of a) n = Γn+a)/ Γa) whenever n + a is not an integer, and the Gauss hypergeometric function and its regularized form with series expansions 7) F a, b c ; z ) := n= a) n b) n c) n z n ) n!, F a, b c ; z := We shall also use the incomplete Beta function and the Beta function, 8) Bx; α, β) := x v α v) β d v, and the regularized incomplete Beta function n= 9) Ix; a, b) := Bx; a, b) / Ba, b). a) n b) n z n, z <. Γn + c) n! Bα, β) := B; α, β), The density η R,q,s of 5) is given in terms of the s-energy of Sd 3, ) W s S d Γd) Γd s)/) ) = x y s d σ d x) d σ d y) = s Γd/) Γd s/), and the Riesz-s potential of the uniform normalized surface area measure σ d evaluated at the location of the source of the external field cf. [5, Theorem ]), ) U σ d s a) = R + ) s s/, d/ F ; 4R / R + ) ), R = a. d By abuse of notation we shall also write U σ d s R). From formula 5) we observe that the minimum value of the density η R,q,s is attained at the North Pole p if q >, and at the South Pole p if q <, η R,q,sp) = + qu σ d s a) W s S d ) q R + )d s W s S d ) R, d η R,q,s p) = + qu σ d s a) W s S d ) q R d s W s S d ) R + ). d Proposition 3. Let < s < d. For the external field Q of ) with q and < R < or R >, the signed s-equilibrium is a positive measure on all of S d if and only if a) Positive external field q > ): ) W s S d ) q R + )d s U σ R d d s R). 3 W s S d ) can be obtained using the Funk-Hecke formula []. Also, cf. Landkof [7].

6 6 J. S. BRAUCHART, P. D. DRAGNEV, E. B. SAFF b) Negative external field q < ): 3) In such a case, µ Q = η Q. W s S d ) q R d s R + ) U σ d d s R). Proof. The arguments given in [5] for q > and R > apply. We provide the proof for 3). If suppµ Q ) = S d, then µ Q is a signed equilibrium on S d the Gauss variational inequalities ) and 3) hold everywhere on S d ). By uniqueness of η Q, η Q = µ Q ; hence, it is non-negative and 3) holds. If 3) holds, then η Q is a non-negative measure on S d whose weighted s-potential is constant everywhere on S d ; that is, η Q satisfies the Gauss variational inequalities with F Q S d ) = G S,Q,s. By uniqueness of µ d Q, µ Q = η Q and suppµ Q ) = S d. Note that R = satisfies ) and 3) with strict inequality for any choice of q. Given a charge q, any R, ), ) for which equality holds in ) or in 3) is called critical distance. At a critical distance R, the density η R,q,s of 5) assumes the value at one point on S d and η R,q,s is strictly positive away from this unique minimum). Interestingly, in the case of negative external fields due to a source inside the sphere, there can be more than one critical distance as discussed below. The critical distances) anchor the subintervals of radii R in, ), ) for which η R,q,s > everywhere on Sd. Remark Positive external fields). For every fixed positive charge q, there is a unique critical distance R q such that for R R q > R R q < ) the signed s-equilibrium is a positive measure on the whole sphere S d. This follows from the fact that the right-most part of ) is a strictly decreasing increasing) function of R if R > < R < ). The technical details for this and the next remark will be postponed until Section 5. Remark Negative external fields). The subtleties of the right-hand side of 3), R d s fr) := R + ) U σ d d s R), gives rise to a multitude of different, even surprising, cases cf. Theorems 6, 7, 8, and 9). The function f is continuous on [, ), negative on, ), and bounded from below. Therefore, 3) is trivially satisfied for all charges q with W s S d )/q < min R fr); otherwise, at least one critical distance exists. For an exterior field source that is, on, )) the function f has the same qualitative behavior for all < s < d in the sense that it is strictly monotonically increasing with lower bound f + ) = W s S d ) and a horizontal asymptote at level. Consequently, there is a unique critical distance if q <, which is the least distance R such that η Q on S d, and none if q. For an interior field source that is, on, )) the qualitative behavior of f changes with the potential-theoretic regime superharmonic, harmonic, subharmonic s). Figure illustrates the typical form of f. In particular, one can have more than one critical value of R for which equality is assumed in 3) as demonstrated for the case d = 4 and s = when 3) reduces to 4) W s S d )/q R) 3 + R) 4 + R /5, R <.

7 AN ELECTROSTATICS PROBLEM ON THE SPHERE R R R Figure. The typical behavior of the function fr) right-hand side of 3)) in the strictly superharmonic, harmonic, and strictly subharmonic case shown here for d = 3 and s = d )/, s = d, and s = d /. The right-hand side above is convex in, ) with a unique minimum at R , so that equality will hold in above relation at two radii R q, and R q, in, ) near this minimum and the above relation will hold on [, R q, ] and [R q,, ). We remark that the complementary problem of fixing the distance R and asking for the critical charge q R is trivial. A simple manipulation yields the unique q R for which equality holds in ) or 3) such that for q > q R > positive external field) or for q < min{q R, } negative external field) the signed s-equilibrium measure is a positive measure on all of S d. Gonchar s Problem. A. A. Gonchar asked the following question cf. [8]): A positive unit point charge approaching the insulated unit sphere carrying the total charge + will eventually cause a spherical cap free of charge to appear. 4 What is the smallest distance from the point charge to the sphere where still all of the sphere is positively charged? For the classical harmonic Newtonian potential s = d ) we answer this question in [5] and discuss it in detail in [7]. In this particular case the s-energy W s S d ) of the d-sphere equals see )), and the mean-value property for harmonic functions implies that the s-potential of the normalized surface area measure σ d at a simplifies to U σ d s R) = R d R ), cf. ) below. So, when requiring that R > and q >, Proposition 3 yields that suppµ Q ) = S d if and only if the following rational relation is satisfied: 5) /q R + ) R ) d R d. The critical distance R q from the center of S d ) is assumed when equality holds in 5). Curiously, for d = classical Coulomb case) the answer to Gonchar s problem is the Golden ratio φ, that is R the distance from the unit sphere) equals + 5)/ 5 ; and for d = 4, the answer is the Plastic constant P defined in Eq. 4) below). For general dimension d, the critical distance R q = R q d) is for positive exterior external fields a solution of the following algebraic equation [ ] 6) G d, q; R) := R ) d /q R R d + R ) d =, 4 On a grounded sphere a negatively charged spherical cap will appear. 5 An elementary physics argument would also show that R φ for q = ) implies suppµ Q ) = S.

8 8 J. S. BRAUCHART, P. D. DRAGNEV, E. B. SAFF which follows from 5) and gives rise to the family of polynomials studied in [7]. In fact, it is shown in [7] that Rq is the uniqe real) zero in, + ) of the Gonchar polynomial G d, q; z). Asymptotical analysis see [7, Appendix A]) shows that Rq = + [log3q)] /d + O/d ) as d. The answer to Gonchar s problem for the external field of ) for general parameters < s < d, R >, and q > relies on solving a, in general, highly non-algebraic equation for the critical distance Rq, namely the characteristic equation cf. relation )) Ws Sd ) R + )d s = Usσd R). d q R ) 7) Figure displays the graphical solution to Gonchar s problem for the dimensions d =, 4. Taking into account that R >, we can rewrite 7) as Figure. Surfaces equality in )) representing the answer to Gonchar question for d =, 4 for the selected ranges for R, q, and s. The dots indicate the solution to Eq. 6) Newtonian case, q = ). 8) G d, s, q; R) =, where we define the Gonchar function of the first kind Ws Sd ) d d s 9) G d, s, q; R) := R ) R + ) Rd + Rd R )d Usσd R). q Answering Gonchar s problem for < s < d, q >, and R > amounts to finding the unique cf. remark after Proposition 3) real)6 zero in, ) of the function G d, s, q; R). 6Depending on the formula used for Usσd R) e.g., )), G d, s, q; R) can be analytically continued to the complex R-plane.

9 AN ELECTROSTATICS PROBLEM ON THE SPHERE 9 Moreover, the density η R,q,s of Theorem evaluated at the North Pole can be expressed as η R,q,sp) = q W s S d ) G d, s, q; R), R >, q >, < s < d. R d d R ) The quadratic transformation formula for Gauss hypergeometric functions [, Eq ] applied to the formula of U σ d s R) of ) yields ) U σ d s R) = /R) s d s)/, s/ F ; /R) ), R >. d + )/ The hypergeometric function simplifies to if s = d harmonic case) and reduces to a polynomial of degree m in /R for d s = m. Thus R q will be an algebraic number if W s S d )/q is algebraic, which is interesting from a number-theoretic point of view. The following result generalizes [7, Theorem 5]. Theorem 4. For the external field Q of ) with < s < d, q >, and R > the signed s-equilibrium is a positive measure on all of S d if and only if R R q, where R q is the unique real) zero in, ) of the Gonchar function G d, s, q; z). If s = d m for m a non-negative integer, then G d, d m, q; z) is a polynomial. In particular, the solution to Gonchar s problem is given by ρd, s) = R q. We remark that integer values of s give rise to special forms of the Gonchar function G. Indeed, for s = d, d 3,..., the s-potential U σ d s and therefore the Gonchar function G reduces to a polynomial. If d s)/ is a positive integer and d is an even dimension, then successive application of the contiguous function relations in [8, 5.5ii)] to ) lead to a linear combination of cf. [8, Eq. 5.4., 5.4.6]) ) /, F 3/ ; z = z log + z z = ) /, z atanh z, F / ; z = z ) with unique coefficients that are rational functions of z, where here z = /R. For the convenience of the reader, we record here that for d = 4 and s =, the s-potential in ) reduces to which can be verified directly by using, for example, MATHEMATICA) U σ 4 R) = 3 R R ) R ) log 8 R 3 R 3 R + ) which yields the Gonchar function of the first kind, ) 3/4 G 4,, q; R) = q R )4 R + ) R R )4 R R + ) R )4 R ) log R ) R + ),

10 J. S. BRAUCHART, P. D. DRAGNEV, E. B. SAFF and for d = 6 and s = 4 or s = one has U σ 6 4 R) = 5 R 4 3 R R + ) R ) R ) log 3 R 4 8 R 5 R + ), U σ 6 R) = 5 R + ) R R + ) 5 R ) 4 R ) log 8 R 4 5 R 5 R + ). This representations hold, in fact, for all R, ), ).) Further analysis shows that for even d and s = d, d 4,..., the Gonchar function of the first kind reduces to G d, d m, q; R) = W ss d ) q R ) d R d + R P R) + QR) log R ) R + ) for some polynomials P and Q. A curious fact is that for odd dimension d 3 and d s)/ a positive integer that is, s = d m), the s-potential U σ d d m R) is a linear combination of a complete elliptic integral of the first kind, π/ d θ Km) := = π ) m sin θ) F /, / ; m, and a complete elliptic integral of the second kind, π/ Em) := m sin θ) d θ = π ) F /, / ; m, with m = /R and with coefficients that are rational functions of /R. 7 This follows by applying contiguous function relations for hypergeometric functions to ) cf. [8, 5.5ii)]). For example, for d = 3 and s = d =, one has U σ d R) = 4 3π + R R ER ) 4 3π R R KR ), R <. If d s)/ is not an integer, then the linear transformation [8, Eq.s 5.8.4] applied to ) followed by the linear transformation [8, last in Eq. 5.8.] gives the following formula valid for all R, ), ), ) U σ d s R) = W s S d )R d R + + ) d s d/, d + s/ F ; d s)/ Γd + )/) Γs d)/) R + R d s R d π Γs/) ) d/, s/ R ) F ; + d s)/ R + ). ) d ) R ) R + ) 7 A similar relation holds if < R <. Then m = R.

11 AN ELECTROSTATICS PROBLEM ON THE SPHERE Both hypergeometric functions reduce to a polynomial if d is an even positive integer. In this case the Gonchar function reduces to G d, s, q; R) = W ss d ) q for some polynomials P and Q. R ) d R d + W s S d ) R + ) d s P R) + Cd, s) R ) d s QR) Gonchar s Problem for Interior Sources. A positive unit) point charge is placed inside the insulated unit sphere with total charge +. What is the smallest distance from the point charge to the sphere so that the support of the extremal measure associated with the external field due to this interior source is just the entire sphere? The trivial solution is to put the field source at the center of the sphere. Then the signed s-equilibrium η Q on S d and the s-extremal measure µ Q on S d associated with the external field Qx) = q/ x s coincide with the s-equilibrium measure σ d on S d. We are interested in non-trivial solutions. First, we answer this question for the classical Newtonian case s = d ). The maximum principle for harmonic functions implies that the s-potential of the s-equilibrium measure σ d is constant on S d and this extends to the whole unit ball Faraday cage effect); that is, U σ d d a) = W d S d ) = for all a R d+ with a. Assuming < R < and q >, by Proposition 3, suppµ Q ) = S d if and only if ) /q + R) R) d. It follows that for a positive external field q > ) induced by an interior source < R < ) the critical distance R q = R q d) to the center of S d ) is a solution of the following algebraic equation 3) G d, q; R) := [ + /q)] R) d R =. For d = classical Coulomb case) the answer to Gonchar s problem is 9 + 8/q R q =, + /q) which reduces to 7 )/4 for q =. This number seems to have no special meaning. 8 The answer to Gonchar s problem for the external field of ) for general parameters < s < d, < R <, and q > relies on solving the characteristic equation cf. )) 4) W s S d ) q = + R)d s U σ R) d d s R). Taking into account that R <, we can rewrite 4) as 5) G d, s, q; R) =, 8 Trivia: The digit sequence of R = 3 + 7)/4 is sequence A88485 of Sloane s OEIS []. One feature is the periodic continued fraction expansion [,, 3].

12 J. S. BRAUCHART, P. D. DRAGNEV, E. B. SAFF where we define the Gonchar function of the second kind, 6) G d, s, q; R) := W ss d ) R) d + R) d s + R) d U σ d s R). q Answering Gonchar s problem for < s < d, q >, and < R < amounts to finding the unique cf. remark after Proposition 3) real) zero in, ) of the function G d, s, q; R). Moreover, the density η R,q,s of Theorem evaluated at the North Pole can be expressed as η R,q,sp) = q G d, s, q; R), < R <, q >, < s < d. W s S d ) R) d The quadratic transformation formula for Gauss hypergeometric functions [, Eq ] applied to the formula of U σ d s R) of ) yields 7) U σ d d s)/, s/ s R) = F ; R ), d + )/ R <. The hypergeometric function simplifies to if s = d harmonic case) and reduces to a polynomial of degree m in R for d s = m. Thus R q will be an algebraic number if W s S d )/q is algebraic. Theorem 5. For the external field Q of ) with < s < d, q >, and < R < the signed s-equilibrium is a positive measure on all of S d if and only if R R q, where R q is the unique real) zero in, ) of the Gonchar function G d, s, q; z). If s = d m for m a non-negative integer, then G d, d m, q; z) is a polynomial. In particular, the solution to Gonchar s problem is given by ρd, s) = R q and the trivial solution ). We remark that in the case of s d m, m a non-negative integer, one can use alternative representations of U σ d s R) similar to those derived after Theorem 4. Connecting Interior and Exterior External Fields. The Riesz-s external fields of the form ) induced by an interior < R <, q > ) and an exterior R >, q > ) point source giving rise to signed s-equilibria η Q and η Q on S d with the same weighted s-potential, U η Q s x) + Q x) = U η Q s x) + Qx) everywhere on S d, are connected by the following necessary and sufficient condition cf. 6)) 8) q U σ d s R ) = q U σ d s R). One way to realize this condition is known as the principle of inversion for S d. For a fixed s, d), the principle states that to a charge q > at distance R > from the center of the sphere there corresponds a charge q = qr s at distance R = /R so that the s-equilibria η Q and η Q coincide and thus have the same weighted s-potential on S d. Indeed, the hypergeometric function in ) is invariant under inversion R /R. The adjustment of the charge follows from 8). An inspection of 5) shows that η Q = η Q.) Now, if we require 8) but do not assume that η Q = η Q, then new phenomena emerge. The potential-theoretic regime superharmonic, harmonic, subharmonic) determines for

13 AN ELECTROSTATICS PROBLEM ON THE SPHERE R R R Figure 3. The typical behavior of the s-potential U σ d s R) for a formula see )) in the strictly superharmonic, harmonic, and strictly subharmonic case shown here for d = 3 and s = d )/, s = d, and s = d /. what ratios q /q, radii R and R exist so that 8) can be satisfied; cf. Figure 3. Basic calculus 9 shows that in the strictly superharmonic case < s < d ) the continuous s-potential U σ d s R) is strictly monotonically decreasing on, ). This implies that for any positive charges q and q with < q /q <, each R, ) uniquely determines an R > such that 8) holds. In the harmonic case s = d ), 8) reduces to q = q R d, since U σ d d R ) = on [, ] cf. 7)), each R, ) is mapped to R = q /q) / d) provided q < q. In the strictly subharmonic case d < s < d) the s-potential U σ d s R) is strictly monotonically increasing on, ) and strictly monotonically decreasing on, ). This implies that for any positive charges q and q with < q /q < W s S d ), each R, ) uniquely determines an R > such that 8) holds. As an example we consider the case d = and s = 3/ and assume that q = q. Then each R, ) determines a unique R > satisfying U σ 3/ R ) = U σ 3/ R), or equivalently, R + R + R + =. R R Curiously, a positive charge at the center of the sphere that is, R = in above relation) would require that R is a solution of the equation R + R )/R = ; that is, one can show that R is the smaller of the two positive zeros of the minimal polynomial x 4 4x Remark. The equations 6) and 3) characterizing the corresponding critical distance R q are related by the principle of inversion; that is, application of the transformations R /R and q qr d to either of the two reproduces the other one. Negative External Fields. A Gonchar type question can also be asked for Riesz external fields Q of the form ) induced by a negative point source. The density η R,q,s of the signed equilibrium on S d associated with Q, given in 5), assumes its minimum value at the South Pole provided R, ), ). Thus, Gonchar s problem concerns the distance from the North Pole such that η R,q,s is zero at the South Pole. Answering this question for negative Riesz external fields is more subtle than for positive external fields. As discussed in the 9 Using in 7) if R < ) and ) if R > ).

14 4 J. S. BRAUCHART, P. D. DRAGNEV, E. B. SAFF second remark after Proposition 3, a solution of the characteristic equation cf. 3)) 9) W s S d ) q R d s = fr) := R + ) U σ d d s R) exists for q satisfying W s S d )/q min R fr) cf. Figure ). For exterior negative point sources R >, < s < d), we have fr) > f) = W s S d ) on, ). Thus, a critical distance R q and therefore an answer to Gonchar s problem) exists if and only if q <. If R q exists, then it is unique.) Defining the Gonchar function of the third kind, Ws S d ) 3) G 3 d, s, q; R) := q ) R + ) d R ) d s R d + R d R + ) d U σ d s R), we can rewrite 9) as G 3 d, s, q; R) =. Given q <, the critical distance R q is the unique real) zero of G 3 d, s, q; R) in, ). Moreover, the density η R,q,s of Theorem evaluated at the South Pole can be expressed as η R,q,s p) = q W s S d ) G 3 d, s, q; R), R >, q <, < s < d. R d d R + ) The following result generalizes [7, Theorem 5] to exterior negative Riesz external fields. Theorem 6. For the external field Q of ) with < s < d, q <, and R > the signed s-equilibrium is a positive measure on all of S d if and only if one of the following conditions holds i) q [, ) and R >. Gonchar s problem has no solution no critical distance). ii) q < and R R q, where R q is the unique real) zero in, ) of the Gonchar function G 3 d, s, q; z). The solution to Gonchar s problem is R q. If s = d m for m a non-negative integer, then G 3 d, d m, q; z) is a polynomial. For an external field with negative point source inside the sphere < R < ), one needs to differentiate between the i) strictly superharmonic < s < d ), ii) harmonic s = d ), or iii) strictly subharmonic d < s < d) case; cf. Figure. We define the Gonchar function of the fourth kind, 3) G 4 d, s, q; R) := W ss d ) q and have η R,q,s p) = + R) d R) d s + + R) d U σ d s R) q G 4 d, s, q; R), < R <, q <, < s < d. W s S d ) R + ) d In case i), the function fr) in 9) has a single minimum in, ), say, at R with fr ) < W s S d ). Set q := W s S d )/fr ). Then 9) has no solution in, ) if q q, ), two solutions in, ) if q, q ) as demonstrated for the special case d = 4 and s = in 4)) which degenerate to one as q q, and one solution in, ) if q.

15 AN ELECTROSTATICS PROBLEM ON THE SPHERE 5 Theorem 7. For the external field Q of ) with < s < d, q <, and < R < the signed s-equilibrium is a positive measure on all of S d if and only if one of the following conditions holds i) q q, ) and R, ). Gonchar s problem has no solution as there exists no critical distance. ii) q = q and R, ). The solution to Gonchar s problem is R as η R,q,s p) = at R = R but η R,q,s p) > for R, R ) R, ). iii) q, q ) and R, R q, ) R q,, ), where R q, and R q, are the two only real) zeros in, ) of the Gonchar function G 4 d, s, q; z). Gonchar s problem has two solutions R q, and R q,. iv) q and R, R q ], where R q is the unique real) zero in, ) of the Gonchar function G 4 d, s, q; z). The solution to Gonchar s problem is R q. v) q < and R = trivial solution). Gonchar s problem has no solution. If s = d m for m a non-negative integer, then G 4 d, d m, q; z) is a polynomial. In case ii) the function fr) is strictly monotonically decreasing and convex on, ) with = W d S d ) fr). Hence, 9) has no solution in, ) if q [, ) and one solution in, ) if q <. Theorem 8. For the external field Q of ) with s = d, q <, and < R < the signed s-equilibrium is a positive measure on all of S d if and only if one of the following conditions holds i) q [, ) and R, ). Gonchar s problem has no solution no critical distance). ii) q < and R, R q ], where R q is the unique real) zero in, ) of the Gonchar function G 4 d, s, q; z). The solution to Gonchar s problem is R q. iii) q < and R = trivial solution). Gonchar s problem has no solution. If s = d m for m a non-negative integer, then G 4 d, d m, q; z) is a polynomial. In case iii) the function fr) is strictly monotonically decreasing on, ) with W s S d ) fr) like in case ii) but neither convex nor concave on all of, ), since f + ) > and f R) as R. Equation 9) has a solution in, ) if and only if q <. Theorem 9. For the external field Q of ) with d < s < d, q <, and < R < the signed s-equilibrium is a positive measure on all of S d if and only if one of the following conditions holds i) q [, ) and R, ). Gonchar s problem has no solution no critical distance). ii) q < and R, R q ], where R q is the unique real) zero in, ) of the Gonchar function G 4 d, s, q; z). The solution to Gonchar s problem is R q. iii) q < and R = trivial solution). Gonchar s problem has no solution. If s = d m for m a non-negative integer, then G 4 d, d m, q; z) is a polynomial.

16 6 J. S. BRAUCHART, P. D. DRAGNEV, E. B. SAFF We remark that in the harmonic case, the analogues of the algebraic equations 6) and 3) characterizing the critical distances) either in, ) or, ) are given by [ ] 3) G 3 d, q ; R) := R + ) d /q R + R d + R + ) d =, R >, q <, 33) G 4 d, q ; R) := [ + /q )] + R) d + R =, R <, q <. Both equations are related by the principle of inversion: simultaneous application of the transformations R /R and q q R d changes one equation into the other. Gonchar s Problem for the Logarithmic Potential. The logarithmic potential log/r) follows from the Riesz s-potential /r s s ) by means of a limit process: log/r) = d r s r s = lim, r >. d s s s s This connection allows us to completely answer Gonchar s problem for the logarithmic potential and the logarithmic external field 34) Q log,a,q x) := q log x a, x R d+, where a S d. We remark that the logarithmic potential with external field in the plane is treated in [6] and [5, 6] deal with the logarithmic case on S which can be reduced to an external field problem in the plane using stereographic projection as demonstrated in [9, 7]. The general case for higher-dimensional spheres seems to have not been considered yet. Theorem. The signed logarithmic equilibrium η Qlog on S d associated with the external field Q log = Q log,a,q of 34) is absolutely continuous with respect to the normalized) surface area measure on S d ; that is, d η Qlog x) = η log x) d σ dx), and its density is given by 35) η logx) = + q q R d x a d, x Sd. Proof. Differentiating the signed equilibrium relation 4) for A = S d with respect to s, i.e. { d η } R,q,s y) d s x y s d σ d y) + d { } q d s x a s = d U η Q s x) + d Qx) = d G S d,q,s, d s d s d s and letting s go to zero yields [ log η x y R,q,s y) ] d σ s dy) + q log x a 36) = d G [ S d,q,s d η ] R,q,s y) d σ d y). d s d s s Using ) and ), one can easily verify that η log y) = lim s η R,q,s y). As the right-hand side above does not depend on x S d, the measure µ with d µy) = η log y) d σ dy) has constant weighted logarithmic potential on S d ; that is, by the uniqueness of the signed

17 AN ELECTROSTATICS PROBLEM ON THE SPHERE 7 equilibrium cf. [5, Lemma 3]), the measure µ is the signed logarithmic equilibrium on S d associated with Q log = Q log,a,q. From 36) and 6) it follows that the weighted logarithmic potential of η Qlog everywhere on S d the constant [ G S d,q log,log = W log S d ) + qu σ d d η ] log a) R,q,s y) d σ d y), d s s equals where the logarithmic energy of S d is explicitly given by cf., e.g., [3, Eq.s.4),.6)]) 37) W log S d ) := {I log [µ] : µ MS d )} = log + [ψd) ψd/)] / here ψs) := Γ s)/ Γs) denotes the digamma function), and the logarithmic potential of σ d can be represented as 38) U σ d log a) = log R + + d/) k 4R) k d) k k R + ). k On the other hand, direct computation of the weighted logarithmic potential at the North Pole gives R G S d,q log,log = + q) W log S d d ) q x a log d x p d σ dx) + q log R. The density η log assumes its minimum value at the North South) Pole if q > q < ), η logp) = + q q k= d R + )d R, R d η log p) = + q q R + ). d From that one can obtain necessary and sufficient conditions for when suppµ Qlog ) = S d. Theorem. Let a R d+ with R = a and q. Set α := [ + q)/q] /d. Then suppµ Qlog ) = S d and thus µ Q = η Q ) if and only if 39) i) for positive fields q > ) + q q R + )d R d, i.e. R α + α R α α + for R > and q >, for < R < and q >, ii) or for negative fields q < ) + q R R + α for R > and q <, d 4) q R + ), i.e. α d R α for < R < and q <, + α whereas the first inequality is trivially satisfied for all R if q <.

18 8 J. S. BRAUCHART, P. D. DRAGNEV, E. B. SAFF Remark. While the critical distance R q,log is given by equality in above relations for q < or q >, in a weak negative logarithmic field q < ) there exists no critical distance. Observe that for fixed charge q, the sequence of critical distances {R q,log } d goes to + ) monotonically as d for q > q < ). Proof of Theorem. By differentiating the expressions in ) and ) with respect to s and letting s, we arrive at 37) and 38). The equivalence in Theorem and µ Qlog = η Qlog ) follows in a similar way as in the proof of Proposition 3. The relations in 39) and 4) are obtained by simple algebraic manipulations. Beyond Gonchar s Problem. This problem arises in a natural way when studying the external field problem on the d-sphere in the presence of a single positive point source above S d exerting the external field. The answer to Gonchar s problem pinpoints the critical distance R q of a charge q > from the center of S d such that the support of the s-extremal measure µ Q associated with Q is all of S d for R R q but suppµ Q ) is a proper subset of S d for < R < R q. Finding the s-extremal measure µ Q when suppµ Q ) S d turns out to be much more difficult. Given d s < d, a convexity argument cf. [5, Theorem ]) shows that S Q := suppµ Q ) is connected and forms a spherical cap centered at the pole opposite to the charge q; in fact, S Q minimizes the F s -functional 4) F s A) := W s A) + Q d µ A, A S d compact with cap s A) >, where W s A) is the s-energy of A and µ A is the s-extremal measure on A. Remarkably, if the signed s-equilibrium on a compact set A S d associated with Q exists, then F s A) = G A,Q,s cf. 4)). This connection to signed equilibria is exploited in [5] when determining the support S Q and, subsequently, the s-extremal measure µ Q on S d associated with Q as the signed equilibrium η Q on a spherical cap Σ t := { u x, u) : u t, x S d } with critical size t = t. We remark that either the variational inequality ) would be violated on S d \ Σ t in case of a too small Σ t t < t ) or the density of η Q would be negative near the boundary of Σ t in case of a too large Σ t t > t ). Somewhat surprisingly it turns out that for s = d the signed equilibrium η Q on Σ t has a component that is uniformly distributed on the boundary of Σ t which vanishes if t = t. It should be noted that η Q on Σ t can be expressed in terms of the s-balayage measures onto Σ t ) Bal s σ d, Σ t ) and Bal s δ a, Σ t ) by means of η Q = c Bal s σ d, Σ t ) q Bal s δ a, Σ t ), where c is chosen such that η Q has total charge. In fact, the balayage method and the restricted if d s < d ) principle of domination play a crucial role in the derivation of these results cf. [5]). Those techniques break down when < s < d d 3). For this gap new ideas are needed. It is the Riesz analog of the Mhaskar-Saff functional from classical logarithmic potential theory in the plane see [] and [6, Ch. IV, p. 94]). Given a measure ν and a compact set A S d, the s-balayage measure ˆν := Bal s ν, A) preserves the Riesz s-potential of ν onto the set A and diminishes it elsewhere on the sphere S d ).

19 AN ELECTROSTATICS PROBLEM ON THE SPHERE 9 Padovan Sequence and the Plastic Number. The Padovan sequence,,,,, 3, 4, 5, 7, 9,, 6,... sequence A93 in Sloane s OEIS []) is named after architect Richard Padovan cf. [8]). These numbers satisfy the recurrence relation P n = P n +P n 3 with initial values P = P = P =. The ratio P n+ /P n of two consecutive Padovan numbers approximates the Plastic constant P as N : P n+ 4) P := lim = 9 69) / ) /3 = n P n /3 3 /3 Padovan attributed [4] its discovery to Dutch architect Hans van der Laan, who introduced in [3] the number P as the ideal ratio of the geometric scale for spatial objects. The Plastic number is one of only two numbers x for which there exist integers k, m > such that x + = x k and x = x m see []). The other number is the Golden ratio φ. Figure 4 shows one way to visualize the Padovan sequence as cuboid spirals, where the dimensions of each cuboid made up by the previous ones are given by three consecutive numbers in the sequence. Further discussion of this sequence appears in [5]. Figure 4. Cuboid spiral visualizing the Padovan sequence. 3. The Polynomials Arising from Gonchar s Problem The discussion of Gonchar s problem for positive/negative Riesz-s external fields with interior/exterior point source leads to the introduction of four kinds of functions. These functions reduce to polynomials if the Riesz-s parameter is given by s = d m for m a non-negative integers. The last two displays in Figure 5 and Figure 6 illustrate zero patterns of the four families of polynomials when m =, which already indicate intriguing features. For example, one notices that the Gonchar polynomials of the second kind have an isolated group of m zeros indicated by in Figures 5 and 6) inside of the left crescent-shaped region, whereas the other zeros are gathered near the right circle. Numerics indicate that the m zeros coalesce at as d increases and the other zeros go to the right circle. In general, fixing m, the zeros of all four polynomials seem to approach One origin of the name is Dutch: plastische getal.

20 J. S. BRAUCHART, P. D. DRAGNEV, E. B. SAFF the three circles and the two vertical line segments connecting the intersection points of the circles as d increasing cf. Figure 5 versus Figure 6). In the following we investigate the four families of Gonchar polynomials G d, q; z), G d, q; z), G 3 d, q ; z), and G 4 d, q ; z) given in 6), 3), 3), and 33). We shall assume that q+q = and q >. Aside from the solution to Gonchar s problem, these polynomials are interesting in themselves and their distinctive properties merit further studies. We studied the family of polynomials G d, q; z) for q = in [7]. Two of the conjectures raised there have been answered in [6]. Interrelations and Self-Reciprocity. The four kinds of Gonchar polynomials are connected. The Gonchar polynomials of the first and second kind are related through 43) 44) z d G d, q /z) d ; /z) = G d, q; z), z d G d, q /z) d ; /z) = G d, q; z). Similarly, the Gonchar polynomials of the third and fourth kind are related through 45) 46) z d G 3 d, q /z) d ; /z) = G 4 d, q ; z), z d G 4 d, q /z) d ; /z) = G 3 d, q ; z). The Gonchar polynomials of the second and fourth kind have in common that for q, 47) G d, q; z) = + ) + z) d + z = G 4 d, q ; z). q For even dimension d and canonical charges q = and q =, the Gonchar polynomials of first and third kind are connected by means of ] 48) z) d G d, ; /z) = ) [ d + z) d + z z d + + z) d = G 3 d, ; z), whereas for odd dimension d and conical charges one has 49) z) d G d, ; /z) + G 3 d, ; z) = + z) d. A polynomial P with real coefficients is called self-reciprocal if its reciprocal polynomial P z) := z deg P P /z) coincides with P and it is called reciprocal if P z) = ±P z). That means, that the coefficients of z k and of z deg P k in P z) are the same. The polynomial G d, ; z) is self-reciprocal for even d, since [ ] 5) z d G d, ; /z) = z) d z z d + z) d. The Gonchar polynomial of third kind is reciprocal for every dimension d; i.e., 5) z d G 3 d, ; /z) = + z)d z d + z d + z d + z) d = G 3 d, ; z). Consequently, if ζ is a zero of G 3 d, ; z), then so is /ζ for any d. The Gonchar Polynomials of the First Kind. We refer the interested reader to [7].

21 AN ELECTROSTATICS PROBLEM ON THE SPHERE d, q d, q [ ] z ) d /q z z d + z ) d = +) [ + /q)] z) d z = ) [ ] z + ) d /q z + z d + z + ) d = ) [ + /q ) ] + z) d + z = ) d, q d, q d, q d, q, m d, q, m Figure 5. Roots of the given polynomial equations for selected d and q q = q). The bottom row shows the roots for the case s = d 3.

22 J. S. BRAUCHART, P. D. DRAGNEV, E. B. SAFF d 4, q, m d 4, q, m d 4, q, m d 4, q, m d 4, q, m d 4, q, m Figure 6. Zeros of the polynomials analogue to the one given in Figure 5 but for the case s = d m; that is, s = d 3 left column) and s = d 4 right column).

23 AN ELECTROSTATICS PROBLEM ON THE SPHERE 3 The Gonchar Polynomials of the Third Kind. The correspondence 48) implies that for even dimension d, the number ζ is a zero of G 3 d, ; z) if and only if /ζ is a zero of G d, ; z). Conjecture. Let Γ be the set consisting of the boundary of the union of the two unit disks centered at and and the line-segment connecting the intersection points Then, as d, all the zeros of G 3 d, q ; z) tend to Γ, and every point of Γ attracts zeros of these polynomials. The Gonchar Polynomials of the Second and Fourth Kind. Both polynomials can be derived from the trinomial 5) P d q; w) := + ) w d + w q by means of a linear transformation of the argument; that is, G d, q; z) = P d q; z), G 4 d, q ; z) = P d q ; + z). Applying results of the Hungarian mathematician Egerváry on the distribution of zeros of trinomials delightfully summarized in [9] and otherwise difficult to come by in the English speaking literature), we derive several results for the Gonchar polynomials of the second and fourth kind. An interesting observation central to Egerváry s work is that the zeros of a trinomial polynomial can be characterized as the equilibrium points of an external field problem in the complex plane of unit point charges that are located at the vertices of two regular concentric polygons centered at the origin. Proposition. Let d and q. Then P d q; w) has simple zeros. Proof. This is clear for d =. Let d. Then, according to [9, Comment after Theorem ], the polynomial equation Az n+m + Bz m + C = has a root with higher multiplicity if and only if ) n+m Am C n B = m m n n n+m n + m) n+m. Taking n = d, m =, A = + /q, B =, and C =, it follows that the left-hand side is negative and the right-hand side positive. Hence P d q; w) has no zero with higher multiplicity. Corollary 3. All the zeros of the Gonchar polynomials of the second and fourth kind are simple. Proposition 4. Let d and q. Assuming that the force is inversely proportional to the distance, the zeros of P d q; w) are the equilibrium points of the force field of positive) unit point charges at the vertices of two concentric polygons in the complex plane. The vertices are given by [ ζ,k := ) ] /d ) d + e i α+k+)π)/d ), k =,..., d, q

24 4 J. S. BRAUCHART, P. D. DRAGNEV, E. B. SAFF and ζ,k := [ + ) d ] /d + e i π α+k+)π)/d, k =,..., d, where α = if q, ), ) and α = π if q, ). Proof. This follows from [9, Theorem ]. q Corollary 5. Let d. If q >, then the zeros of G d, q; z) are the equilibrium points of the force field of positive unit point charges at the vertices ζ,k k =,..., d ) and ζ,k k =,..., d). If q < q ), then the zeros of G 4 d, q ; z) are the equilibrium points of the force field of positive unit point charges at the vertices ζ,k k =,..., d ) and ζ,k k =,..., d). The paper [9] and the more recent [9] discuss annular sectors as zero inclusion regions. Proposition 6. Let d 3. Suppose s and t are the unique positive roots of x d x = and x d + x =, respectively. Choose radii ρ and ρ with /3 ρ t < < s ρ 3/) /d ). Then each of the annular sectors { } 53) z C : ρ z ρ, πk argz) d 3/)/d ), k =,..., d, d contains exactly one of the d zeros of P d ; w) = w d + w. Proof. The substitution w = C/A) /d ζ for any d-th root of C/A transforms the equation Aw d + Bw + C = into ζ d aζ =, where a = B/C) C/A) /d. Hence P d q; w) = if and only if ) 54) ζ d aζ =, where a = /d ) /d + and w = + ζ. q q Let s and t, with < t < < s, be the unique positive roots of x d a x = and x d + a x =, respectively. Then, by [9, Thm. 3.], every root of 54) lies in the annulus {ζ C : t ζ s}. It can be readily seen that a + a ) /d ) < for d 3 and q >. Let d 3 and q >. Then each of the disjunct sectors { ζ C : πk argζ) d θ d }, k =,..., d, where sin θ = a + a ) /d ), contains exactly one root of 53) by [9, Thm. 5.3]. For q =, the coefficient a reduces to / and both P d ; w) = w d +w and ζ d aζ have the same zeros. The result follows. The bounds for ρ and ρ follow from [9, Lem..6]. Proposition 6 can be obtained for P d q; w) with general q. Figure 7 illustrates the force field setting and the zero inclusion regions for the canonical case q = and d = 6).

25 AN ELECTROSTATICS PROBLEM ON THE SPHERE 5 d 6, q Figure 7. The zeros ) of P 6 ; w) = w 6 + w, lying in the annular sectors of 53) with ρ = t and ρ = s, are equilibrium points of the force field due to the vertices of the concentric pentagon and hexagon cf. Prop. 4). The resultant of the force stream plot) is 5 k= z ζ,k) + 6 k= z ζ,k). Theorem 5 and Theorem 8, respectively, imply the following properties of the Gonchar polynomials of the second and fourth kind. Proposition 7. If q >, then G d, q; z) has a unique real zero in the interval, ). Proposition 8. If q, ), then G 4 d, q ; z) has no zero in, ). If q =, then G 4 d, q ; z) = z. If q <, then G 4 d, q ; z) has a unique real zero in the interval, ). From the fact that P d q; ζ) = implies ζ = + q ζ d, where the right-hand side is changing exponentially fast as d when the zeros avoid an ε-neighborhood of the unit circle, we get that the zeros of P d q; w) approach the unit circle as d. This in turn implies that the zeros of G d, q; z) for q > ) approach the circle Γ := {z C : z = } and the zeros of G 4 d, q ; z) for q < and q ) approach the circle Γ := {z C : z + = } as d.

26 6 J. S. BRAUCHART, P. D. DRAGNEV, E. B. SAFF 4. Negatively Charged External Fields Signed Equilibrium on Spherical Caps We are interested in the external field due to a negative charge below the South Pole, 55) Q b,s x) = q x b s, x Sd, b = Rp, R >, q <, that is sufficiently strong to give rise to an s-extremal measure on S d that is not supported on all of the sphere. In [6] we outline how to derive the signed equilibrium on spherical caps centered at the South Pole and, ultimately, the s-extremal positive) measure on S d associated with the external field 55). Here we present the details. For the statement of the results we need to recall the following instrumental facts: We assume throughout this section that s d. 3 The signed s-equilibrium on a spherical cap Σ t := {x S d : p x } associated with Q b,s can be represented as the difference 56) η t = Φ st) W s S d ) ν t qε t, in terms of the s-balayage measure 4 onto Σ t of the positive unit point charge at b and the uniform measure σ d on S d given by 57) ɛ t = ɛ t,s := Bal s δ b, Σ t ), ν t = ν t,s := Bal s σ d, Σ t ). Furthermore, the function 58) Φ s t) := W ss d ) + q ε t ), d < s < d, ν t where ɛ t = S d d ɛ t and ν t = S d d ν t, plays an important role in the determination of the support of the s-extremal measure on S d. In particular, one has 59) U ηt s x) + Q b,s x) = G Σt,Q b,s,s = Φ s t) = F s Σ t ). Here, F s is the functional of 4) for the field Q b,s. First, we provide an extended version of [6, Theorem 9] that includes asymptotic formulas for density and weighted potential valid near the boundary of the spherical cap. Proposition 9. Let d < s < d. The signed s-equilibrium η t on the spherical cap Σ t S d, < t <, associated with Q b,s in 55) is given by 56). It is absolutely continuous in the sense that for x = u x, u) Σ t, 6) d η t x) = η tu) ω d ω d u ) d/ d u d σd x), 3 When d =, s = d is understood as the logarithmic case s = log. 4 Given a measure ν and a compact set K of the sphere S d ), the balayage measure ˆν := Bal s ν, K) preserves the Riesz s-potential of ν onto the set K and diminishes it elsewhere on the sphere S d ).

27 AN ELECTROSTATICS PROBLEM ON THE SPHERE 7 where with R = b and r = R + Rt + ) 6) η tu) = ) d/ Γd/) t t u W s S d ) Γd s/) u t, d/ {Φ s t) F d s)/ ; t u ) u q R )d s F r d ) s d)/, d/ R + ) ; d s)/ r )} t u. u The density η t is expressed in terms of regularized Gauss hypergeometric functions. As u approaches t from below, we get ) s d)/ η tu) Γd/) t u = W s S d ) Γd s/) Γ d s)/) t 6) { } R )d s Φ s t) q + Rt) t u) + Ot u) ) as u t, r d where 63) Rt) = d 64) 65) d s d s Φ s t) q R )d s r d t d R )d s q d s r d R r. Furthermore, if z = ξ z, ξ) S d, the weighted s-potential is given by U ηt s z) + Q b,s z) = Φ s t), z Σ t, U ηt s z) + Q b,s z) = Φ s t) + q ρ s I R ) ξ t Φ s t) I r + ξ ; d s, s ξ t + ξ ; d s, s ) ), z S d \ Σ t, where ρ = R + Rξ + and Ix; a, b) is the regularized incomplete Beta function. As ξ approaches t from above, we get ) { d s)/ ξ t Γd/) Us ηt z) + Q b,s z) = Φ s t) + + t Γ + d s)/) Γs/) 66) [ ] } R )d s q Φ r d s t) + Oξ t) as ξ t +. The next remarks, leading up to Proposition, emphasize the special role of Φ s t). Remark. The behavior of the density η t near the boundary of Σ t inside Σ t determines if η t is a positive measure; namely, the signed equilibrium η t on Σ t associated with Q b,s is a

28 8 J. S. BRAUCHART, P. D. DRAGNEV, E. B. SAFF positive measure with support Σ t if and only if R ) d s 67) Φ s t) q R + Rt + ). d/ Indeed, relations 6) and 63) show that 67) is necessary and sufficient for η tu) > in a sufficiently small neighborhood t ε, t), ε >, and this inequality extends to all of [, t) as shown after the proof of Proposition 9. Note that the density η tu) has a singularity at u = t if Φ s t) qr ) d s /R + Rt + ) d/ and η tu) approaches as u t when equality holds in 67). In that case, however, d η t d u = Γ + d/) Γd s/) Γ d s)/) t u t ) s d)/ { Φs t) + Ot u) W s S d ) Remark. The weighted s-potential of the signed equilibrium η t on Σ t associated with Q b,s exceeds the value Φ s t) assumed on Σ t strictly outside of Σ t but on S d ) if and only if R ) d s 68) Φ s t) q R + Rt + ). d/ Indeed, the expansion 66) shows that U ηt s z) + Q b,s z) > Φ s t) in a small neighborhood of the boundary of S d \ Σ t if and only if 68) holds. In addition, if Φ s t) satisfies 68), then 69) Us ηt z) + Q b,s z) > Φ s t) everywhere on S d \ Σ t. This inequality is shown after the proof of Proposition 9.) The weighted s-potential of η t tends to Φ s t) when the boundary of S d \ Σ t is approached from the outside ξ t + ). There is a vertical tangent if Φ s t) qr ) d s /R + Rt + ) d/ which turns into a horizontal one at a t for which equality holds in 69). In such a case d { } ) d s)/ { } ξ t Γ + d/) R Φ Us ηt s t) z)+q b,s z) = + Oξ t). d ξ + t Γ + d s)/) Γs/) R + Rt + For η t to coincide with the s-extremal measure on S d associated with Q b,s with support Σ t both 68) and 67) have to hold. The next result is [6, Theorem ]. Proposition. Let d < s < d. For the external field 55) the function Φ s t) given in 58) has precisely one global minimum t c, ]. This minimum is either the unique solution t c, ) of the equation Φ s t) = q R ) d s / R + Rt + ) d/, or t c = when such a solution does not exist. In addition, Φ s t) is greater than the righthand side above if t, t c ) and is less than if t t c, ). Moreover, t c = max{t : η t }. The extremal measure µ Qb,s on S d is given by η tc see 6)), and suppµ Qb,s ) = Σ tc. Figure 8 illustrates the typical behavior of the signed equilibrium using density and weighted potential) on spherical caps that are too large, too small, and have just the right size. The right column shows which conditions are violated when the spherical cap is too small or too large. This figure should be compared with [5, Fig. ]. }.

29 AN ELECTROSTATICS PROBLEM ON THE SPHERE 9.8 t t c ' Η t Η t Qb,s U s t > t c, Us ηt z) + Q b,s z) F s Σ t ) on S d \ Σ t, Us ηt z) + Q b,s z) = F s Σ t ) on Σ t, η t on Σ t. c.8.84 Η t Qb,s U s.86 t = tc, Η t ' Us ηt z) + Q b,s z) F s Σ t ) on S d \ Σ t, Us ηt z) + Q b,s z) = F s Σ t ) on Σ t, η t on Σ t. c.8 Η t '.84 Η t Qb,s U s.86 t < t c, t c Us ηt z) + Q b,s z) F s Σ t ) on S d \ Σ t, Us ηt z) + Q b,s z) = F s Σ t ) on Σ t, η t on Σ t. c Figure 8. The weighted s-potential of η t for t > t c, t = t c, and t < t c versus altitude ξ of z = ξ z, ξ) S d for d =, s =, q = 5, and R = + φ φ the Golden ratio), cf. Propositions 9 and. Insets show the respective density η t. The horizontal line Red in colored version) indicates Φ s t c ) = F s Σ tc ) = G Qb,s,s. Observe the vertical tangent at the graph of η t as t t c in the middle display see first remark after Proposition 9). In the limiting case s = d with s > it can be shown that the s-balayage measures 7) ɛ t := ɛ t,d = Bal d δ b, Σ t ), ν t := ν t,d = Bal d σ, Σ t ) exist and both have a component that is uniformly distributed on the boundary of Σ t. Moreover, unlike the case d < s < d, the density for µ Qb,s, where s = d, does not vanish on the boundary of its support. We introduce the measure β t x) := δ t u) σ d x), x = u x, u) S d.

30 3 J. S. BRAUCHART, P. D. DRAGNEV, E. B. SAFF Using similar methods as in [5], one can show the following. Theorem. Let d 3. The signed s-equilibrium η t on the spherical cap Σ t associated with Q b,d x) = q x b d, q < and b =, R) R > ), is given by η t = [ Φ d t)/w d S d ) ] ν t qɛ t, Φ d t) := W d S d ) + q ɛ t ) / ν t, where ν t and ɛ t are given in 7). More explicitly, for x = u x, u) S d, 7) d η t x) = η tu) d σ d Σt x) + q t d β t x), where the density with respect to σ d restricted to Σ t takes the form η tu) = Φ d t) W d S d ) q R ) W d S d ) R + Ru + ) d/+ and the boundary charge uniformly distributed over the boundary of Σ t is given by q t = t ) [ ] t d/ q R ) Φ d t). R + Rt + ) d/ Furthermore, for any fixed t, ), the following weak convergence holds: 7) ν t,s ν t, ɛ t,s ɛ t, as s d ) +. The function Φ d t) has precisely one global minimum t c, ]. This minimum is either the unique solution t c, ) of the equation Φ d t) = q R ) / R + Rt + ) d/, or t c = when such a solution does not exist. Moreover, t c = max{t : η t }. The extremal measure µ Qb,d on S d with suppµ Qb,d ) = Σ tc is given by [ ] 73) d µ Qb,d x) = d η t x) = Φ d t ) R + ) R + Rt c + ) d/ d σ Σt W d S d ) R + Ru + ) d/+ d x). 74) 75) Furthermore, if z = ξ z, ξ) S d, the weighted d )-potential is given by U η t d z) + Q b,d z) = Φ d t), z Σ t, [ U η t d z) + Q b,d z) = Φ d t) + q ρ d [ Φ d t) R ) R + Rt + ξ t + ξ ξ t + ξ ) d/ ] ) d/ ], z S d \ Σ t, where ρ = R + Rξ +. As ξ approaches t from above, we get { [ ] 76) U η t d z) + Q b,d z) = Φ d t) + ξ t d R ) q Φ + t r d d t) A similar result holds for the logarithmic case s = log on S. + Oξ t) }.

31 AN ELECTROSTATICS PROBLEM ON THE SPHERE 3 5. Proofs 5.. Proofs and Discussions for Section. First, we show the result for the signed s-equilibrium on S d. Proof of Theorem. For R > this result has been proven in [5] cf. [7] for the harmonic case). Let R <. By linearity of the s-potential, we can write [ U η Q s x) = + qu ] σ d s a) U σ q y a s R ) d d σ d d y) W s S d s x) ) W s S d ) R ) s y x s y a d y a. d Using Imaginary inversion i.e., utilizing y a y a = R and d σ d y)/ y a d = d σ d y )/ y a d ), the integral reduces to W s S d )/ x a s. Since U σ d s x) = W s S d ) on S d, one gets that the weighted potential of η Q is constant, U η q Q s x) + x a s = W s S d ) + q U σ d s a) everywehre on S d. This shows 6). In a similar way, d η Q = + qu σ d s a) W s S d ) q W s S d ) y a s R ) d d σ d y) R ) s y a d y a =, d where the integral reduces to U σ d s a) under Imaginary inversion. Thus η Q, indeed, satisfies Definition for A = S d. Furthermore, the signed measure η Q is absolutely continuous with respect to σ d. The result follows. Next, we provide the technical details for the discussion in the two remarks after Proposition 3. First remark after Proposition 3. For R >, we can use term-wise differentiation in the following series expansion of the right-hand side of ) to show monotonicity, [ s/) ] k d/)k 4R) k d) k k! R + ). s+k k= For R, ), we show that the following representation of the right-hand side of ) + R) d s [ R) d F d s/, d/ ; 4R )], d R) obtained by applying the linear transformation [8, last of Eq.s 5.8.] to ), is strictly monotonically increasing on, ) for each s, d). It is easy to see that the ratio + R) d s / R) d has this property for each s, d) and using the differentiation formula [8, Eq. 5.5.], the square-bracketed expression above has a positive derivative on, ) for each s, d). The result follows.

32 3 J. S. BRAUCHART, P. D. DRAGNEV, E. B. SAFF Second remark after Proposition 3. The continuity of f is evident by continuity of the s-potential U σ d s a), which is a radial function depending on R = a, in the potentialtheoretical case < s < d. Note that, since σ d is the s-equilibrium measure on S d, f) = U σ d s ) = W s S d ) <. The negativity of f in, ), ) follows from [ ] U σ d R ) d s)/ s R) = R + ) d F d s/, d/; d ) 4R > R + ) R d s R + ) d, where the right-hand side is derived from ) using the last transformation in [8, Eq. 5.8.]. Clearly, fr) U σ d s R). From ) it follows that U σ d s R) is strictly monotonically decreasing on, ). Since σ d is the s-equilibrium measure on S d, fr) U σ d s ) = W s S d ) = f) on [, ). Since the continuous function f is bounded on the compact interval [, ], f is bounded from below on [, ). Verification of monotonicity of f on, ) by direct calculation of f seems to be futile, but by using ) and the differentiation formula [8, Eq. 5.5.], we get s R + ) s fr) + R + ) s f R) = {R + ) s fr)} ) d s R = d s) R + R + ) + s R ) R + ) 3 F + s/, + d/ 4R ; >. + d R + ) Since it has already been established that fr) <, it follows that f R) > on, ). It is easy to see that fr) as R, thus f has a horizontal asymptote at level. Next, we prove results for Gonchar s problem for negative external fields. The following proof concerns in particular Theorem 7. Proof. We consider interior sources, that is R <. The right-hand side in 3) is fr) := gr) U σ d s R), where gr) := R) d s + R) d. We show that the equation f R) =, or equivalently, g R)/R = {U σ d s R)} /R has only one solution in the interval, ) in the superharmonic regime < s < d. Expanding the last equation using formula 7) for U σ d s R), we get [d s + R)] R)d s + R) d+ R = sd s) d + d s)/, + s/ F ; R ). + d + )/ The function hr) := g R)/R at the left-hand side satisfies hr) + as R and h) =. Since {R + R) d+ hr)} = {[d s + R)] R) d s } <, it follows that h R) <. The right-hand side of the equation assumes the positive value sd s)/d+) at R = and is strictly decreasing if < s < d 3, constant if s = d 3, or strictly increasing if d 3 < s < d. In either case for < s < d there is exactly one solution

33 AN ELECTROSTATICS PROBLEM ON THE SPHERE 33 in the interval, ). That is, fr) has a single minimum in the interval, ) because f ) = s d < and [8, Eq. 5.4.]) ) f sd s) d s)/, + s/ ) = F ; d + + d + )/ sd s) Γ + d + )/) Γd s) = d + Γd s/) Γd + s)/) = s Γd + )/) Γd s) Γd s/) Γd + s)/) = s W ss d ) >, < s < d. The last step follows by applying the duplication formula for the gamma function and ). In particular, the minimum from above is strictly less than f ) = W s S d ). The function gr) is strictly monotonically decreasing on, ) for all < s < d. From g R) = R)d s R + ) d+ { [d + ) s + )R + )] [d sr + )] + s R ) } one infers that gr) changes from convex to concave as R and g R) as R if d < s < d. Since the s-potential U σ d s R) is a strictly increasing and convex function on, ) in the subharmonic range d < s < d, it follows that fr) is strictly decreasing on, ) and neither convex nor concave on all of, ) for d < s < d. Moreover, f ) > and f R) as R for d < s < d. 5.. Proofs and Discussions for Section 4. In the following we make use of a Kelvin transformation spherical inversion) of points and measures that maps S d to S d. Let K R denote the Kelvin transformation stereographic projection) with center b =, R) and radius R ; that is, for any point x R d+ the image x := K R x) lies on a ray stemming from b, and passing through x such that 77) x b x b = R. The image of x = u x, u) S d is again a point x = u ) x, u ) S d, where the formulas 78) + u = R ) R + Ru + u), u = R + ) + u), R + Ru + relating the heights u and u, follow from similar triangle proportions. From this and the formula x y = x y for x, y S d it follows that the Euclidean distance of two points on the sphere transforms like 79) x y = R ) x y x a y a, x, y Sd. Geometric properties include that the Kelvin transformation maps the North Pole p to the South Pole q and vice versa, K R S d ) = S d, and K R sends the spherical cap

34 34 J. S. BRAUCHART, P. D. DRAGNEV, E. B. SAFF Figure 9. Kelvin transform with center at a and radius R. A R := { u x, u) : /R u, x S d } to B R := { u x, u) : u /R, x S d } and vice versa, with the points on the boundary being fixed. We note that the uniform measure σ d on S d transforms like 8) x b d d σ d x ) = x b d d σ d x). Further, given a measure λ with no point mass at b, its Kelvin transformation associated with a fixed s) λ = K b,s λ) is a measure defined by 8) d λ x ) := R ) s/ x b s d λx), where the s-potentials of the two measures are related as follows e.g. [, Eq. 5.)]) d λ 8) Us λ x y ) x b s s d λy) x b ) = x y s = = R ) s/ s x y R ) U λ s/ s x). The Kelvin transformation has the duality property K b,s λ x )) = λx).) For convenience, we recall the specific form of the s-balayage onto Σ t of the uniform measure see [5, Lemma 4]) and its norm see [5, Lemma 3]). Proposition. Let d < s < d. The measure ν t = Bal s σ d, Σ t ) is given by 83) d ν t x) = ν tu) ω d ω d u ) d/ d u d σd x), x Σ t, where the density ν tu) is given by 84) ν tu) := Γd/) Γd s/) t u ) d/ ) s d)/ t u t, d/ F d s)/ ; t u u ).