arxiv: v1 [math.co] 8 Sep 2008

Transcription

1 APPELL POLYNOMIALS AND THEIR ZERO ATTRACTORS ROBERT P. BOYER AND WILLIAM M. Y. GOH arxiv: v [math.co] 8 Sep 2008 Abstract. A polynomial family {p n(x)} is Appell if it is given by ext = P g(t) n=0 pn(x)tn or, equivalently, p n(x) = p n (x). If g(t) is an entire function, g(0) 0, with at least one zero, the asymptotics of linearly scaled polynomials {p n(nx)} are described by means of finitely zeros of g, including those of minimal modulus. As a consequence, we determine the limiting behavior of their zeros as well as their density. The techniques and results extend our earlier work on Euler polynomials. Contents. Introduction 2. The Generalized Szegö Approximations 3 3. Asymptotics Outside the Disk D(0; /r 0 ) 4 4. Asymptotics on the Basic Regions R l 5 5. Dominant Zeros and Szegö Curves 0 6. Asymptotics Inside the Disk D(/a ; δ a ) When a is a Non-Dominant Zero 3 7. Zero Attractor and the Density of the Zeros 4 Appendix A. Density of Zeros 8 A.. Introduction 8 A.2. Contributions over Arcs 9 A.3. Backlund s Method 20 A.4. Main Density Theorem 22 A.5. Special Case 23 References 23. Introduction Let g(t) be an entire function such that g(0) 0. Definition. The Appell polynomials {p n (x)} associated with generating function g(t) are given by () e xt g(t) = n=0 p n (x)t n. Some important examples are: the Taylor polynomials of e x, with g(t) = t; the Euler polynomials, with g(t) = (e t + )/2; and Bernoulli polynomials, with g(t) = (e t )/t; and their higher order analogues. The asymptotics and limiting behavior of the zeros of these families have been investigated by many people; for example, [2], [7], and so on. Date: September 8, Mathematics Subject Classification. Primary 05C38, 5A5; Secondary 05A5, 5A8. Key words and phrases. Appell polynomials, zeros of polynomials, asymptotics.

2 In this paper, we obtained the asymptotics and the limiting behavior of the zeros for all Appell families provided the generating function g(t) satisfies one further condition: that g must have at least one zero. We use the ideas in our earlier paper [2]; furthermore, we simultaneously simplify and generalize some of the techniques there. We found that the asymptotics in the general case are build from the basic example g(t) = t which coincides with the classical work of Szegö on the Taylor polynomials of the exponential function. In our paper [2], we found that the asymptotics for the Euler and Bernoulli polynomials are controlled by certain roots of g(t), the ones of minimal modulus. In the general situation, as expected, the minimal modulus roots of g(t) are needed to describe the asymptotics but there may be finitely many other roots needed to determine the asymptotics. These additional roots are determined through a geometric condition described in the terms of rotated and scaled versions of the Szegö curve: xe x =, x, x C (see Figure 3). We frequently use the following notations. Let Z(g) denote the set of all zeros of g and let r 0 < r < r 2 <... denote the distinct moduli of these zeros in increasing order. Recall that if K and K 2 are two non-empty compact subsets of C, then their Hausdorff distance is the larger of sup{d(x, K ) : x K 2 } and sup{d(x, K 2 ) : x K }. Definition 2. For a family {q n (x)} of polynomials whose degrees are increasing to infinity, their zero attractor is the limit of their set of zeros Z(q n ) in the Hausdorff metric on the space of all non-empty compact subsets of the complex plane C. In the appendix, we discuss how the zero attractor is found in terms of the limsup and liminf of the zero sets Figure. Zeros for degree 000 polynomial, with generating function g(t) = J 0 (t) There is related work on the asymptotics and zeros of the Taylor polynomials for linear combinations of exponentials c j e λ jx where the parameters λ j satisfy a geometric constraint []. The techniques of proof are very different from our approach. 2

3 2. The Generalized Szegö Approximations It is convenient to collect together several results from [2] and some extensions of them concerning the asymptotics of S n (x) = n k=0 xk /k!. The domains of where their asymptotics hold are critical in understanding the behavior for the Appell polynomials. Proposition. (Left-Half Plane) Let /3 < α < /2 and j. On any compact subset K of {w : Rw < }, we have () S n (nw) e nw = (we w ) n ( + O(n 3α ) ), ( w) (2) Dw j (w n S n (nw)) = Dw j (w n e nz (j )! e n ( ) + O(n 3α ( w) j ) ), where the big O constant holds uniformly for x K. The proof of part (a) is in [2]. Part (b) follows from an application of the saddle point method. The following Proposition is also from [2]: Proposition 2. (Outside Disk) Let S be a subset contained in w > with distance δ > 0 from the unit circle, and let α be chosen so /3 < α < /2. Then where the big O term holds uniformly for w S. S n (nw) e nw = (we w ) n (w ) ( + O(n 3α ) ), Proposition 3. (Evaluations of Integrals) If ɛ < w and j, then we have ( ) e () 2πi t =ɛ t t w dt = w n S n (wxn). ( ) e (2) dt = 2πi t =ɛ t (t w) j (j )! Dj w (w n S n (wxn)), where D w is the differentiation operator d dw. Proof. (a) By expanding /(t z) into an infinite geometric series and performing a term-by-term integration, we obtain ( ) e ( ) e dt = 2πi t =ε t t z z2πi t =ε t t dt z = ( ) e ( ) t m dt. z2πi t =ε t z m 0 By the Cauchy integral theorem the terms correspond to m n vanish. Hence ( ) ( e ) dt = 2πi t t z z z m e xtn t n+m dt 2πi t =ε = z n m 0 n m 0 = z z n+ = z n t =ε (xn) n m z m (n m )! n m 0 n j 0 3 (xnz) j j! (xnz) n m (n m )! = ( )z n S n (zxn).

4 Part (b) follows from differentiating (a) j times with respect to z. 3. Asymptotics Outside the Disk D(0; /r 0 ) Theorem. Let K be any compact subset in the annulus A(/r 0 ; ). We have holds uniformly for x K. (xe) n / = ( + O(/n)). g(/x) Proof. We shall find an asymptotic approximation to in the region A(/r 0 ; ) = Use the generating relation equation () to get p n (x) = 2πi t =ɛ e xt dt, g(t)tn+ {x : x > r0 }. where ɛ < r 0. Since both sides of the above equation are entire functions of x, by analytic continuation this representation for p n (x) is valid for all x C. Hence we can replace x by nx to get (2) = 2πi t =ɛ ( e xt t ) n dt tg(t). The above expression is valid for 0 < ɛ < r 0 and is the starting point of the analysis in the sequel. Let K be an arbitrary compact subset {x : x > r 0 } and let x K. We can certainly choose ɛ small enough so that for all x K, ɛx <. By a change of variables, we get = xn 2πi t =ɛ x ( e t t ) n dt tg(t/x). Observe that the zeros of g(t/x) have the form ax where a Z(g). Moreover, they must lie outside the closed unit disk since x > /r 0, so we can deform the integration path from the circle with radius ɛ x to the unit circumference. Thus = xn 2πi = xn 2πi t = t = x ) ( e t t ) n dt tg(t/x) n(t ln t) dt e tg(t/x). It can be easily seen that t = is the saddle point of the integral and the classical saddle point method is applicable here [4]. Hence = ( (ex)n + O( n ) g( ), where the implied O constant holds uniformly for x K. The last equation can be written as (ex) n / = g( x ) ( + O( n ) ), x > /r 0. We have the: Corollary. (a) On the complement of the disk D(0; /r 0 ), lim n n ln (ex) n / = 0 where the limit holds uniformly on compact subsets. (b) The zero attractor must be contained in the closed disk D(0; /r 0 ). 4

5 Note that part (b) follows easily from (a) since g(x) never vanishes on the disk D(0; /r 0 ). 4. Asymptotics on the Basic Regions R l Let r 0, r,... denote the distinct moduli of the zeros a of the generating function g. Fix an integer l. We fix ρ > 0 so it is not equal to any zero modulus {r 0, r,... }. For each zero a Z(g) with a = r l, we consider the circle x = / a and the disk D(/a, δ a ). Now the tangent line T a to the circle x = / a at the point /a determines the half-plane H a, which contains 0; that is, R(ax) <. We choose ɛ l > 0 to be less than the distance from the portion of the tangent line T a that lies outside the disk D(/a; δ a ) to the circle x = / a for any a = r l+ ; that is, ɛ l < /rl 2 + δ2 a /r l. Finally, we make the requirement the disks D(/a; δ a ) be mutually disjoint for all a Z(g) with a < ρ. Definition 3. With these conventions, the region R l is described in terms of the half-planes H a and disks as (3) R l = { Ha \ D( a ; δ a) : a = r l } \ D(0; r l+ + ɛ l+ ) We note that the regions R l are not disjoint; in fact, by construction, its inner boundary which consists of the portion of the circle x = r l+ + ɛ that lie outside the disks D(/a; δ a ), a = r l+, actually lies inside the region R l+. It is convenient to introduce a region that contains all of the R l s. Definition 4. Let R ρ be the domain given as (4) R ρ = [ ] {H a : a Z(g), a = r 0 } \ {D(/a; δa ) : a Z(g), a < ρ} D(0; /ρ). Note the order of dependence: first we can given the cut-off modulus ρ > 0 for the moduli of the zeros; next, δ a > 0 for each a Z(g) is given and is a function of ρ [see later section], then finally, ɛ l is determined relative to each zero moduli r l which is a function of δ a. For any a Z(g) with r 0 a < ρ, let s a (t) be the singular part of at its pole a. Next we set g (t) to be tg(t) (5) g (t) = tg(t) {s a (t) : a Z(g), r 0 a < ρ} we see that g (t) is analytic for t < ρ. We develop the asymptotics for {} on the regions R l where r 0 r l < ρ. Now we saw already that we can write as = 2πi t =ɛ ( e xt where s(t) = {s a (t) : a Z(g), r 0 a ρ}. t ) n g (t) dt + 2πi t =ɛ ( ) e s(t) dt, Lemma. With g (t) given above in equation 5, we have ( ) e g (t) dt = xn e n g (/x) ( + O(/n)) 2πi t t =ɛ uniformly on compact subsets of the annulus A(/ρ; ). 5 t

6 Proof. Let x K A(/ρ, ). By a change of variables, we write ( ) e ( ) g (t) dt = xn e t n g (t/x) dt. 2πi t 2πi t t =ɛ t =ɛ x By construction, g(t/x) is analytic on a disk of radius greater than. So the contour in the last integral can be deformed to the unit circle x = without changing its value. Finally, by an application of the saddle point method we find that x n ( ) e t n g (t/x) dt = xn e n g (/x) ( + O( 2πi t n )). t = Figure 2. Generic Plot of Polynomial Zeros and Zero Attractor When g Has Two Roots; Tangent Lines and Circles Displayed To state the next two lemmas, we need to introduce special polynomials I n (z) in z and J(a; z) in z. The polynomial I n (z) comes from expanding the derivative of Dz m (z n e nz ). Consider m ( ) m Dz m (z n e nz ) = (D p p zz n )(Dz m p e nz ) = p=0 m p=0 ( m p m = z n e nz n m ) ( n)( n ) ( n p + )z n p (n m p e nz ) p=0 p=0 ( m p ) ( n)( n ) ( n p + )(nz) p m ( )( ) m n + p = z n e nz n m ( ) p p! (nz) p p p = z n e nz n m I m (nz), 6

7 where I m (z) is given in Definition 5. m ( )( ) m n + p (6) I m (z) = ( ) p p! z p. p p p=0 For a Z(g), we define J(a; z) which are also polynomials in z. We write out the singular part s a (t) of the function tg(t) at its nonzero pole a by (7) s a (t) := β a m= where β a is the order of a as a zero of g(t) so b a,βa 0. b a,m (t a) m, Definition 6. For a Z(g), let J(a; z) be the polynomial in z given as (8) J(a; z) = β a m= b a,m (m )! zm I m (az). Lemma 2. Let a Z(g) and let x K, a compact subset of the half-plane H a, R(ax) <. Then ( ) e s a (t) dt = a n e nax J(a; nx) + en x n s a (/x) ( + O(n 3α ) ) 2πi t t =ɛ where s a (t) is the singular part of /(tg(t)) at the zero a of g(t). Proof. We first write out the integral in terms of the singular part s a (t) ( ) e β a b a,m ( s a (t) dt = 2πi t (m )! Dm a a n S n (nax) ) t =ɛ m= where the coefficients b a,m are given in equation (7). We now study the asymptotics of the typical term Da m (a n S n (nax)). We may use the generalized half-plane Szegö asymptotics with 3 < α < 2 because of the restriction that a Z(g) with a r l to obtain D m a (a n S n (nax)) = x n+m Dax m ((ax) n S n (nax)) { = x n+m Dz m (z n e nz ) e n z=ax Combining these estimates we obtain (9) D m a (m )! (a n S n (nax)) = x n+m {(ax) n e nax n m I m (nax) (m )! e n ( + O(n 3α ( ax) m ) ) } = a n e nax (nx) m I m (nax) Hence after summation we obtain ( ) e s a (t) dt 2πi t = (0) t =ɛ β a m= ( + O(n 3α ( ax) m ) ) }. (m )! e n x n+m ( + O(n 3α ( ax) m ) ). b a,m (m )! Dm a (a n S n (nax)) = a n e nax J(a; nx) + en x n s a ( x ) ( + O(n 3α ) ). 7

8 Corollary 2. For a Z(g), a r l, we have ( ) e s a (t) dt = a n e nax J(a; nx) + en x n s a (/x) ( + O(n 3α ) ) 2πi t t =ɛ uniformly on the compact subsets of R l, where s a (t) is the singular part of /(tg(t)) at the zero a of g(t). Lemma 3. Let a Z(g) and let x K, where K compact subset of the disk-complement A(/ a ; ). Then ( ) e s a (t) dt = en x n s a (/x) ( + O(n 3α ) ). 2πi t =ɛ t Proof. We will use the disk-complement generalized Szegö asymptotics. A( + c, ), for any c > 0, we have S n (nz) = zn en(ζ ln ζ) dζ 2πi ζ z ζ = By Dividing z n and taking derivatives up to order m, we get Dz m (z n (m )! en(ζ ln ζ) S n (nz)) = 2πi (ζ z) m dζ In the above, replace z by ax to obtain = (m )! ζ = e n ( z) m ( + O(n 3α ) ). Da m (a n S n (nax)) = x n+m Dax m ((ax) n S n (nax)) (m )! e n x n+m ( () = + O(n 3α ( ax) m ) ). By summation, we obtain the asymptotics for the original integral: ( ) e β a b a,m s a (t) dt = 2πi t (m )! Dm a (a n S n (nx)) (2) t =ɛ m= = en x n s a ( x ) ( + O(n 3α ) ). Corollary 3. For a Z(g) with r l+ a < ρ, we have ( ) e s a (t) dt = en x n s a (/x) ( + O(n 3α ) ), 2πi t t =ɛ uniformly on the compact subsets of R l. For z in the annulus The remaining case for the above integration involving s a (t) on the disk D(/a; δ) will be handled in a later section. Proposition 4. For x R l, we have (xe) n / = x g (/x) { } J(a; nx) φ(ax) n : a Z(g), a r l + { ( ) } x s a : a Z(g), a < ρ ( + O(n 3α )) x 8

9 uniformly on the compact subsets of R l, where φ(x) = xe x and /3 < α < /2. Proof. Putting the last two corollaries into Equation (2) and simplifying, we have = xn e n g ( ( x ) + O( ) n ) { a n e nax J(a; nx) : a Z(g), a < ρ } ( + O(n 3α ) ). Proposition 5. For x R l, we have (xe) n / = g(/x) {J(a; nx) φ(ax) n : a Z(g), a r l} + O(n 3α ) uniformly on the compact subsets of R l, where φ(x) = xe x and /3 < α < /2. Proof. By the definition of g (t) (see equation (5)), we see that ( ) (3) x g = x g(/x) { x s a( } ) : a Z(g), a < ρ. x Insert this into the above. Since the s a (/x) term cancels, we have uniformly for x R l : (4) (ex) n / = g(/x) { (axe ax ) n J(a; nx) : a Z(g), a < ρ } + O ( n 3α). Lemma 4. If a Z(g) with a < ρ and x R l, then J(a; nx) = b ( ) a,β a ax βa (β a )! (nx)βa ( + o()). ax Proof. Recall that J(a; nx) = β a m= It is easy to see that that is, as n Hence m b a,m I m (nax), I m (nax) = (m )! ( n + p p p=0 p=0 )(nax) p = (ax) p p! ( ) p ( m p ( + o()), )( n + p m ( ) ( ) m ax m I m (nax) ( ) p (ax) p =. p ax (5) J(a; nx) = b a,β a (β a )! (nx)βa ( ax ax ) βa ( + o()). p ) p! (nax) p. Since the coefficient b a,βa in the definition of the singular part s a (t) is nonzero, we find for fixed x that the precise order of J(a; nx) as a polynomial in n is n βa. We note the following 9

10 Corollary 4. a r l. lim n (xe) n / = g(/x), x R l provided φ(ax) > for all a Z(g) with Corollary 5. lim n n ln (xe) n / = ln φ(a 0x), for all x R l satisfying φ(a 0 x) > for a 0 Z(g) and φ 0 (ax) φ(ax) for all a a 0 such that a r l. By construction, R l R ρ for all l chosen so r l < ρ. Consequently, we have the asymptotics: Theorem 2. On R ρ, we have the following uniform asymptotics where /2 < α < /3. (ex) n / = g(/x) ( + O(/n)) { (axe ax ) n J(a; nx) : a Z(g), a < ρ } + O ( n 3α), It remains to develop the asymptotics in the disks D(/a; δ a ) and well as determining domination among a Z(g) of φ(ax). 5. Dominant Zeros and Szegö Curves Let φ(x) = xe x which is an entire function that is conformal on the unit disk. The standard Szegö curve S is the portion of the level curve φ(x) = that lies inside the closed unit disk or equivalently, inside the closed left-hand plane R(x). S is a closed simple closed convex curve; in fact, it has the form t = ± e 2(x s) s 2 where x = s + it and s [ W (e ), ] and W is the principal branch of the Lambert W -function. Definition 7. Let a be a nonzero complex number. We call any curve of the form a S a Szegö curve x Figure 3. Szegö Curve: ze z = and z 0

11 Lemma 5. Let a, b be non-zero distinct complex numbers. Then the intersection a S b S has at most 2 points. Proof. The intersection of the two curves a S b S must satisfy φ(ax) = φ(bx). This modulus condition determines a line; so, the intersection of the two curves lie on the as and a line. Since as is convex, the intersection contains at most two points. Write x = s + it and b a = α + iβ. Then φ(ax) = φ(bx) reduces to the line: ae ax = be bx, e (b a)x = b/a, e R[(b a)x] = b/a, R[(b a)x] = ln b/a, αs βt = ln b/a. Definition 8. Call a Z(g) a dominant zero if either a = r 0 or /a does not lie in the interior of any of Szegö curves b S where b is a minimal modulus zero, b = r 0. If /a lies on one of these Szegö curves, call a an improper dominant zero; otherwise, a is proper. Note: if a is an improper dominant zero, then there exists a minimal modulus zero b so the intersection a S b S consists of a single point, namely /a. In [], the behavior of the zeros of the Taylor polynomials of m j= c je λjz is determined if the numbers λ j satisfy a condition in terms of a convex m-gon. For the Appell polynomials, the geometric condition is more subtle. Let W denote the principal value of the Lambert W -function. Since the radius of the largest circle centered 0 that lies in the interior of the standard Szegö curve S is W (e ), we have the following Lemma 6. If a is a zero of g such that a > r 0 /W (e ), then a must be a non-dominant zero. Hence, there are at most finitely many dominant zeros. Definition 9. If a and b are two dominant zeros such that b S as consists of two points, then the line φ(ax) = φ(bx) determines two half-planes. Let E a +,b denote the closed half-plane that contains /a. Let a, a 2,..., a n be the dominant zeros of g. Recall that if f(x) is any analytic function on a domain D then its modulus f(x) is a subharmonic function on D. We need two basic properties of subharmonic functions: they satisfy the maximum modulus principle; and the maximum of two subharmonic functions is still subharmonic. Definition 0. Set Φ(x) = max( φ(a x),, φ(a n x) ) so Φ(x) is subharmonic on the punctured complex plane. We work with φ(ax) rather than φ(ax) so we may apply the Maximum Principle since φ(ax) vanishes at x = 0. We observe that the level curve Φ(x) = divides the complex plane into finitely many connected components just as the original curve φ(x) = divides the complex plane into three connected components. We make the: Definition. Consider the connected components of the complement of the level curve Φ(x) =. Let D 0 be the closure of the connected component that contains 0. Note that if x D 0 and is x 0, then Φ(x) with strict inequality when x lies in the interior of D 0. We assume that ρ > 0 is chosen so large that the closed disk D(0; /ρ) is a subset of the interior of D 0. Furthermore, the singularity in log φ(ax) log φ(bx) always cancels so this difference is always a harmonic function on C. Let a be a non-dominant zero of g such that a < ρ. Then we know that a S is a subset of the interior of D 0 and that the disk D(0; /ρ) lies inside of a S. The following Proposition follows easily from the definition of D 0 and dominant zeros:

12 Proposition 6. Uniformly on the compact subsets of R ρ \ D 0, we have lim n (ex) n / = g(/x) To understand the asymptotics inside D 0 requires subharmonic function theory. By construction, this means that φ(a x) < for x D 0, while Φ(x) = for x D 0. Before we can apply the Maximum Modulus Principle for subharmonic functions, we need to deal with the common singularity at 0. However, this is easily dealt with by multiplying both Φ(x) and φ (x) by x which shows that we can remove this singularity relative to the inequality. In particular, Φ(x) is strictly larger than φ (x) for all x D 0, with x 0. So there exists a positive constant, say α, so that Φ(x) > φ(a x) + α, We state the above discussion formally: x D 0 \ D(0; /ρ). Proposition 7. Let a be a non-dominant zero of g. On the domain D 0 \ D(0; /ρ), we have the order estimates φ(a x) n = O((Φ(x) α) n ) = o(φ(x) n ). Definition 2. Given a proper dominant zero a of g, let D a = {x D 0 : φ(ax) φ(bx), for all dominant zeros b a, φ(ax) }. Note that the definition of D a is independent of the choice of ρ sufficiently large. Further, D a has the alternate description in terms of the half-planes E a +,b: a Interior(S) {E a +,b : b is a proper dominant zero of g, b a}. It is easy to state formally the basic properties of D a : Lemma 7. Let a be a proper dominant zero of g. Then D a is a non-empty compact connected subset of D 0. Proof. Since a S is a convex curve and the intersection of half-planes is connected, the set D a must be a connected convex set. We now restate Proposition 7 relative to a domain D a : Proposition 8. Let a be a non-dominant zero of g with a < ρ. Then a must lie in a domain D a for some proper dominant zero a of g. For x D a \ D(0; /ρ), the following holds uniformly φ(a x) n = O( φ(ax) α n ) = o( φ(ax) n ). Finally, we now have a refinement of Theorem 2 as Theorem 3. Let ρ be chosen greater than / a where a is any proper dominant zero of g. Then on R ρ, we have the following uniform asymptotics (ex) n / = g(/x) ( + O(/n)) { φ(ax) n J(a; nx) : a Z(g) and dominant } + O ( n 3α) + o(φ(x)), where /3 < α < /2 and Φ(x) = max{ φ(ax) : a Z(g) and dominant }. 2

13 6. Asymptotics Inside the Disk D(/a ; δ a ) When a is a Non-Dominant Zero We first state an easy consequence of a previous Proposition. Proposition 9. Let a Z(g). Then on the disk D(/a ; δ a ), the normalized polynomials have the asymptotics (xe) n / = ( ) ( ( )) x g + O + { } J(a; nx) x n φ(ax) n : a Z(g), a a, a a + { ( ) } x s a : a Z(g), a a, a < ρ + σ a (x), x where σ a (x) = β a m= b a,m (m )! Dm a ( (a ) n S n (na x) ). Proposition 0. Let a be a non-dominant zero of g with a < ρ. Then there exists a choice of δ a > 0 such that σ a (x) = O (e 6nδ ρ) a where ρ is the cut-off bound for the zeros of g. Proof. To estimate σ a, we make use of the elementary estimate: If f(z) is analytic function of z, then for any ɛ > 0, we have D z j f(z) (j )! ɛ j max f(ζ). ζ z =ɛ By the definition of σ a (x), we find β a σ a b a,m (m )! Dm a ( a n S n (na x)) m= β a b a,m max ζ n S n (nζx) ζ a =δ a m= δ m a K δa max ( ζ n S n ( ζx n)) ζ a =δ a where K δa > 0 is a constant that depends on the zero a and the radius δ a. To go further we observe for x D( a, δ a ) and ζ a = δ a : ζx ( a + δa ) x a x + x δa + a δ a + x δ a = + δ a ( a + x ). Since a < ρ by assumption, ζx + 2ρδ a. But ζ a δ a and x a δ a, so we get ζx ( a ) ( ) ( δ a a δ a δ a a + a ) 2δ a ρ. Collecting these two inequalities, we get 2δ a ρ ζx + 2δ a ρ. 3

14 Now use that S n (nt) e nt : max eζx n S n ( ζx n) e n 2δ a ρ n e n(+2δ a ρ) ζ a =δ a = 2δ a ρ n e 2nδ a ρ For 0 x /2, /( x) e 2x ; if we choose δ a such that 2δ a ρ /2, then we have 2δ a ρ n e 4δ a ρ. With this choice of δ, we obtain the desired bound max ( eζx n S n ( ζx n)) e 4δ a ρ e 2nδ a ρ = e 6nδ a ρ ζ a =δ a 7. Zero Attractor and the Density of the Zeros In our paper [2], we determined the limit points of the zeros of the Euler polynomials by means of the asymptotics and the zero density. Here, we separate out first the question of find the support of the zero density measure, which is, of course, the zero attractor. Then we determine the zero density by applying our general result in the appendix. Proposition. Let f n (x) = /(xe) n. Then the following limits hold uniformly on compact subsets of the indicated domains: () On the domain A(/r 0 ; ), lim n n ln[f n(x)] = 0. (2) On the domain R ρ, lim n n ln[f n(x)] = 0. (3) On the domain D a A(/ρ; ) where a is any dominant zero of g, lim ln φ(ax). n n ln[f n(x)] = Proof. We use the asymptotic expansions for developed in the previous sections. For x > /r 0, we noted already that the indicated limit must be 0. We observe that if a is a nondominant zero of g with a < ρ, then for δ a > 0 sufficiently small, the disk D(/a ; δ a ) will lie in the domain D a for some dominant zero a, then on D(/a ; δ a ), lim n n ln[f n(x)] = ln φ(ax). To describe the zero attractor requires a closer examination of the boundary of each domain D a where a is a proper dominant zero. The boundary D a where a is a proper dominant zero of g has several natural families: D a D 0 which is an outer boundary and a polygonal curve consisting of the line segments D a D b where b is another dominant zero of g. Note that D a D b is a subset of D 0. It will be useful to subdivide D a D 0 into two connected components denoted by D a ± that come from deleting {/a} from [ D a D 0 ]. Lemma 8. The zero attractor of the Appell polynomials {} must lie inside the compact set { Da : a is a proper dominant zero of g}. Proof. First, we let x let in the infinite exterior of D 0. Recall that lim n pn (nx)/(xe) n = /g(/x) uniformly on compact subsets. If x nk is a zero of p nk (n k x) and x nk x, then appealing to this limit we find that the limit must be 0 while the right-hand side is /g(/x ) 0. Secondly, suppose x lies in the interior of D 0 but not on any boundary set D a, where a is a dominant zero. By construction, x will lie in the interior of one of the domains D b, where b is a dominant zero. Then lim n /(xe) n /n = φ(bx) uniformly on compacta in the interior of D a. By the same reasoning as before, x cannot be a limit of zeros. 4

15 The following Theorem is an immediate consequence of the above lemma together with the result of Sokal in section A.4 of the Appendix. Theorem 4. The zero attractor of the Appell polynomials {} is given by { Da : a is a proper dominant zero of g}. Proof. Let a be any proper ominant zero of g and let x D a ±. Let ɛ > 0 be given. Then we find that { lim ln n (xe) n / = 0, x D(x ; ɛ) \ D 0, ln φ(ax), x D(x ; ɛ) Int(D 0 ) holds uniformly on compact subsets. Next suppose that x is nonzero and lies on one of the line segments of the form D a D b where b is another proper dominant zero. Again, we find that { lim ln n (xe) n / ln φ(ax), x D(x = ; ɛ) Int(D a ), ln φ(bx), x D(x ; ɛ) Int(D b ) which also holds uniformly on compact subsets. By Sokal s result [6] which is described in the appendix, we conclude that x is in lim sup Z(p n ) since there can be no harmonic function v(x) on the disk D(x ; ɛ) that satisfies the inequalities lim inf n ln (xe) n / v(x) lim sup n ln (xe) n /. On the disks D(/a ; δ a ) where a is a non-dominant zero which must lie inside D 0, the contribution of the dominant zero dominants. This reasoning handles all but finitely many points: /a where a is a dominant zero of g. However, since the zero attractor must be a compact set and points in D(/a; ɛ) [ D a D 0 ] lie in the zero attractor, we conclude that /a also lie in the attractor. Theorem 5. Let a be a proper dominant zero of g. (a) The zero density measure on any proper subcurve of D a D 0 is the pull-back of the normalized Lebesgue measure on the unit circle under the conformal map φ(ax). (b) Let b be a proper dominant zero of g so b a. Then the zero density measure on any proper line segment of D a D b is Lebesgue measure. Proof. For both parts, we can use the asymptotics given in Theorem 3. For part (a), let f n (x) = g(/x)/(xe) n. Let a be a dominant zero of g, and let C be a proper subcurve of D 0 D a ±. Then there exists a neighborhood U of C such that U R ρ [(C \ D 0 ) D a ] so that the asymptotics in Theorem 3 can be written as (xe) n / = g(/x) ( + O(/n)) J(a; nx) φ(ax) n + O(n 3α ) + o(φ n,a(x)), where Φ,a (x) = max{, φ(ax) }. Hence, by dividing by g(/x)), we find that f n (x) has the form: where f n (x) = + a n (x)φ(ax) n + e n (x), { e n (x) = a n (x) = g(/x)j(a; nx), o(), x U (C \ D 0 ), o(φ(ax) n ), x U (D a R ρ ). Since φ(ax) is conformal in the disk D(0; / a ), we may apply Theorem 6 from the Appendix Section A.4 on the density of zeros. 5

16 Let a and b be two distinct proper dominant zeros of g such that D a D b is nonempty. On D a D b R ρ, the asymptotics in Theorem 3 can be written as (xe) n / = g(/x) ( + O(/n)) ( J(a; nx) φ(ax) n + J(b; nx) φ(bx) n + ) {J(a ; nx) φ(a x) n : a proper dominant zero, a a, b } + O(n 3α ) + o(φ(x) n ) = g(/x) ( + O(/n)) ( ) J(a; nx) φ(ax) n + J(b; nx) φ(bx) n + O(n 3α ) + o(φ n a,b (x)), where Ψ a,b (x) = max{/ φ(ax), / φ(bx) }. Let L be a proper line segment of the intersection D a D b. Let U be a neighborhood of L so both φ(ax) < and φ(bx) < for x U. On the intersection U R ρ, we work with a different normalization than before: φ(ax) n T n (x) = (xe) n J(a; nx) p n(nx). Note that in this normalization the term that contains φ(ax) n becomes the constant for T n (x). Of course, this new normalization has exactly the same zeros as in U so the zero density is unchanged. Then we find that where and T n (x) = + a n (x)ψ(x) n + e n (x), ψ(x) = φ(ax) φ(bx) = a b e(b a)x, a n (x) = J(b; nx) J(a; nx), φ(ax) n ( e n (x) = O(n 3α ) + o(φ n a,b (x))). J(a; nx) On U, we have that φ(ax) n Φ n a,b (x) = max{, ψ(x) n }; while on D a U, ψ(x) < and on D b U, ψ(x) >. This allows us to write e n (x) as o(ψ(x) n ), x D a U, e n (x) =. o(), x D b U. By construction, φ(ax)/φ(bx) = a b e(b a)x is a conformal map on U R ρ that maps L onto an arc of the unit circle. By Corollary 7 in the Appendix section A.4, the result follows. We close with several examples that illustrate the main constructions in the paper. Example. Let g(t) be an entire function whose minimal modulus zero a = such that all its other zeros a satisfy / a < W (e ) Then the zero attractor for the associated Appell polynomials coincide with the classical Szegö curve in Figure 3. Example 2. The higher order Euler polynomials E n (m) (x), where m Z +, have generating function g(t) = (e t +) m /2 m ; while the higher order Bernoulli polynomials B n (m) (t) have generating function g(t) = (e t ) m /t m. Then their zero attractors are independent of m and coincide with a scaled version of the zero attractor for the Taylor polynomials for cos(x), see Figure 4. Example 3. The zero attractor for the Appell polynomials associated with generating function g(t) = J 0 (t), where J 0 (t) is the zero-th order Bessel function, is a scaled version as the zero attractor for the Taylor polynomials for cosh(x), see Figure, since the minimal modulus zeros of J 0 (t), a = , are the only dominant zeros since all the zeros of J 0 (t) lie on the real axis. 6

17 Figure 4. Zero Attractor for Taylor polynomials of cos(x) Example 4. Let g(t) = (t ) (t 2 + 2). See Figure 5 for its zero attractor and zeros for degree 400. Figure 5. (a) Zeros for degree 400 polynomial with generating function g(t) = (t ) (t 2 + 2); (b) Zero Attractor with polynomial zeros Example 5. Consider the Appell polynomials with generating function g(t) = (t a)(t b)(t c) with a =.2e i3π/6, b =.3e i7π/6, and c =.5. In this case, all three roots of g(t) are dominant. See Figures 6 and 7. These last two examples both illustrate the following general fact. We assume that the generating function g(t) has exactly three proper dominant zeros a, b, and c. Then the three lines determined by φ(ax) = φ(bx), φ(ax) = φ(cx), and φ(bx) = φ(cx) have a common intersection point, a 7

18 Figure 6. (a) Zero Attractor only, for generating function g(t) = (t a)(t b)(t c), a =.2e i3π/6, b =.3e i7π/6, c =.5; (b) Boundary of the Domain D 0. Figure 7. Zeros for degree 400 polynomial together with the Zero Attractor, for generating function g(t) = (t a)(t b)(t c), a =.2e i3π/6, b =.3e i7π/6, c =.5 so-called triple point. This follows by interpreting the lines as the boundary between the change of asymptotics of the Appell polynomial family; that is, the boundaries of the domains D a, D b, and D c. Appendix A. Density of Zeros A.. Introduction. We generalize the density result for the zeros of the Euler polynomials in [2] to highlight how the asymptotic structure of the polynomial family may determine the density of its zeros. 8

19 Let ψ(x) be an analytic function on a domain D C that is conformal on D. We write ζ = ψ(x). We sometimes write x(ζ) for x = ψ (ζ). We assume that there exists ɛ 0 > 0 and 0 α < β 2π so that the annular sector (6) S = {ρe iθ : ρ [ ɛ 0, + ɛ 0 ], θ [α, β]} lies in the image ψ(d). Next we define two subsectors of S as S + = {ρe iθ : ρ [ ɛ 0, ), θ [α, β]} S = {ρe iθ : ρ (, + ɛ 0 ], θ [α, β]}. Let C be the unimodular curve ψ ({e iθ : θ [α, β]}), so φ(x) = for x C. By construction, C is smoothly parametrized as x(e iθ ) for θ [α, β]. Of course, we have ψ (S) = ψ (S ) C ψ (S + ) as a disjoint union. Let {T n (x)} be a sequence of analytic functions on ψ (S) where we assume that the analytic functions satisfy the basic asymptotic relation: (7) T n (x) = + a n (x)ψ(x) cn + e n (x), where {c n } is a increasing unbounded sequence of positive numbers, δ > 0 is a constant so that a n (x) δ, and a n (x) = exp[o(c n )], uniformly on ψ (S). The term e n (x) satisfy the following estimates uniformly: { o(ψ(x) c n ), x S e n (x) = +, o(), x S. In the sequel, we may assume either form for e n (x) if x lies on the common boundary C of the two regions S ± that is, ψ(x) =. Let Z n be the set of all zeros of T n that lie in ψ (S), which we assume is finite for all n. For [γ, γ 2 ] (α, β), let (8) N n (γ, γ 2 ) = #{x Z n : arg x [γ, γ 2 ]}. Choose ɛ > 0 so 3ɛ < ɛ 0. By the Argument Principle, we find that N n (γ, γ 2 ) = d dζ T n(x(ζ)) dζ 2πi T n (x(ζ)) where Γ is the boundary of the sector {ρe iθ : ρ [ ɛ, + ɛ], θ [γ, γ 2 ]}. The closed contour Γ naturally has four parts of the form Γ ±ɛ and Γ γj, j =, 2 where A.2. Contributions over Arcs. Lemma 9. lim n Γ ɛ Γ ±ɛ = {( ± ɛ)e iθ : θ [γ, γ 2 ]}, d dζ T n(x(ζ)) dζ = 0. T n (x(ζ)) Γ γ = {ρe iγ : ρ [ ɛ, + ɛ]}. Proof. Let ζ Γ ɛ. Then x(ζ) ψ (S ) and e n (x(ζ)) 0 uniformly on Γ ɛ. Furthermore, we find a n (x(ζ))ψ n (x(ζ)) exp[o(c n )] ζ cn = e cn ln(/( ɛ))+o(cn)) = o(). From the basic asymptotic expression T n (x(ζ)) = + a n (x(ζ))ζ cn + e n (x(ζ)), we find that ln[t n (x(ζ))] converges uniformly to 0 on Γ ɛ. Hence their derivatives must also converge uniformly to 0 and so the desired integrals converge to 0. 9 Γ

20 ( Lemma 0. lim I n c n Γ +ɛ d dζ T ) n(x(ζ)) dζ = γ 2 γ. T n (x(ζ)) Proof. For ζ Γ +ɛ, ζ = + ɛ and x(ζ) ψ (S + ). By the basic expansion T n (x(ζ)) = + a n (x(ζ))ζ cn + e n (x(ζ)), we have that T n (x(ζ)) a n (x(ζ))ζ = + + e e n(x(ζ)) cn a n (x(ζ))ζ = + ζ cn + n(x(ζ)) ζ cn. cn a n (x(ζ)) We recall our assumptions that a n (x(ζ)) δ and e n (x(ζ))/ζ cn = o(). Hence, we find that + e n (x(ζ)) a n (x(ζ))ζ cn = ζ cn + en(x(ζ)) ζ cn a n (x(ζ)) ( + ɛ) cn + o() δ = o() T n(x(ζ)) uniformly on Γ +ɛ. In particular, a n(x(ζ))ζ converges to uniformly on Γ cn +ɛ so ln[ converges uniformly to 0 there as do their derivatives. In other words, we know that d dζ T n(x(ζ)) T n (x(ζ)) d dζ a n(x(ζ)) c n a n (x(ζ)) ζ 0 Note that a n (x) = exp[o(c n )] implies that a n(x)/a n (x) = o(c n ). Hence, we find We now conclude that The lemma now follows easily. d dζ T n(x(ζ)) T n (x(ζ)) c n = = d dζ a n(x(ζ)) + c n a n (x(ζ)) ζ + o() d dx a n(x) dx a n (x) dζ + c n ζ + o() = c n ζ + o(c n). d dζ T n(x(ζ)) T n (x(ζ)) ζ uniformly on Γ +ɛ. T n(x(ζ)) a n(x(ζ))ζ ] cn A.3. Backlund s Method. Our estimates for the integrals over Γ γ, for γ [γ, γ 2 ], are inspired by the 98 method of R. Backlund s proof of the Riemann-von Mangoldt asymptotic formula for the number of zeros of the Riemann zeta function. We follow the exposition of Chandrasekharan [3] (pages 35-38). Recall that Γ γ is parametrized as ρe iγ for ρ [ ɛ, + ɛ]. Let l be the number of zeros of R[T n (x(ζ))] for ζ Γ γ exclusive of endpoints. Then the contour integral along Γ γ can be written as a sum of integrals over line segments C ab paramerized as ρe iγ, ρ [a, b] where a and b are two consecutive zeros of R[T n (x(ζ))]. In particular, R[T n (x(ρe iθ ))] has constant sign for ρ [a, b]. Then ( d dζ I T ) ( ) n(x(ζ)) dξ dζ = I, T n (x(ζ)) C n ξ C ab where C n is the image of the line segment ρe iγ, for ρ [a, b], under the map ζ T n (x(ζ)). By assumption, C n can only intersect the imaginary axis ir only at its two endpoints; in particular, C n must lie either in the left and right half-plane. By Cauchy s theorem, we can deform C n into a 20

21 semicircle K ab that lies in the same half-plane and has the same endpoints on ir so that the value of the contour integral is unchanged. This allows us to make the estimate I dξ C n ξ = dξ I K ab ξ π. We summarize this discussion as: Lemma. I ρ [ ɛ, + ɛ]. Γ γ d dζ T n(x(ζ)) T n (x(ζ)) dζ (l + )π where l is the number of zeros of R[T n(x(ρe iθ ))] for We now use Jensen s formula to make a useful estimate for l. For an analytic function h(z) on some domain E, define h(z) on E c as h(z) = h(z) where E c = {z : z E} which will be analytic on E c. For ξ D( ɛ; 2ɛ), which is symmetric about the real axis, let [ ˆT n (ξ) = 2 T n (x(ξe iγ )) + T ] n ( x(ξe iγ )). Then ˆT n (ξ) is analytic on D( ɛ; 2ɛ) and ˆT n (ξ) = R[T n (x(ξe iγ ))], ξ [ ɛ, + ɛ] R. For convenience, we recall Jensen s inequality. Let h(z) be an analytic function on the closed disk D(a; R), and let 0 < r < R. Suppose h(a) 0. Let m be the number of zeros of h(z) in the closed disk D(a; r) counted according to their multiplicity. Then ( R r ) m max{ h(z) : z a = R}. h(a) Now each zero of R[T n (x(ζ))] for ζ Γ γ corresponds to a zero of ˆT n (ξ) for ξ [ ɛ, + ɛ]. Let ˆl be the number of zeros of ˆT n (ξ) for ξ D( ɛ; 2ɛ). Then we have at once the inequality l ˆl. We will apply Jensen s inequality to the disk D( ɛ; 3ɛ) and r = 2ɛ to obtain ( 3 )ˆl max{ ˆT n (ξ) : ξ ( ɛ) = 3ɛ} 2 ˆT. n ( ɛ) Lemma 2. (a) ˆT n ( ɛ) = + o(). (b) max{ ˆT n (ξ) : ξ ( ɛ) = 3ɛ} = O (e o(cn) ( + 2ɛ) cn ). Proof. (a) Since ( ɛ)e iγ S, we have the estimate T n (x(( ɛ)e iγ )) + e o(cn) ( ɛ) cn + o() = + o(). In particular, ˆT n ( ɛ) = + o(). (b) To estimate the maximum of ˆT n (ξ) for ξ ( ɛ) = 3ɛ, we observe that ˆT [ n (ξ) 2 T n (x(ξe iγ )) + T ] n ( x(ξe iγ )). Let ζ lie in the closed disk D( ɛ; 3ɛ) so ζ + 2ɛ. By the basic asymptotic relation, we find that whether ζ S + or S : T n (x(ζ)) + a n (x(ζ))ζ cn + e n (x(ζ)) + e o(cn) ( + 2ɛ) cn + o(( + 2ɛ) cn ) = O 2 (e o(cn) ( + 2ɛ) cn ),

22 where the big-oh constant holds uniformly on D( ɛ; 3ɛ). In particular, this estimate holds for ξ that lie on the circle ξ ( ɛ) = 3ɛ. A similar estimate holds for T n ( x(ξe iγ )). We sum up this discussion as max{ ˆT ) n (ξ) : ξ ( ɛ) = 3ɛ} = O (e o(cn) ( + 2ɛ) cn. The last lemma together with Jensen s inequality allows us to make an estimate for ˆl: ( 3 )ˆl max{ ˆT n (ξ) : ξ ( ɛ) = 3ɛ} ) = O (e 2 ˆT o(cn) ( + 2ɛ) cn. n ( ɛ) Recalling that l ˆl, we have the bound Hence, we have the estimate for the integral d I dζ T n(x(ζ)) dζ T n (x(ζ)) c n Γ γ l ln[o ( e o(cn) ( + 2ɛ) cn) ]. ln(3/2) ( ) π ln( + 2ɛ) + O. ln(3/2) c n Since these bounds hold for all ɛ > 0 sufficiently small, we have shown the Lemma 3. For all ɛ > 0 sufficiently small, d I dζ T n(x(ζ)) dζ T n (x(ζ)) c n where Γ γ is the line segment ρe iγ, ρ [ ɛ, + ɛ]. γ 2 γ 2π Γ γ Finally, this last inequality allows us to make the estimates ( ) ln( + 2ɛ) + O 2 ln(3/2) c n ( ) π ln( + 2ɛ) + O ln(3/2) c n N n(γ, γ 2 ) γ 2 γ + 2πc n 2π ( ) ln( + 2ɛ) + O. 2 ln(3/2) c n A.4. Main Density Theorem. By combining the above lemmas and noting that these bounds hold for all ɛ > 0 sufficiently small, we obtain our main density result. Theorem 6. Let α < γ < γ 2 < β, and let N n (γ, γ 2 ) denote the number of zeros of T n (x) whose arguments lie in [γ, γ 2 ], given in equation (8). Then N n (γ, γ 2 ) lim = γ 2 γ ; n c n 2π that is, the image of the zero density under ψ is Lebesgue measure on an arc of the unit circle. We need to recall the notions of lim sup and lim inf of a sequence {X n } of compact sets in the complex plane. Now x lim sup X n if for every neighborhood U of x, there exists a sequence x nk X nk U that converges to x while x lim inf X n if for every neighborhood U of x, there exists an index n and a sequence x n X n U, for n n that converges to x. It is known that if the lim inf X n and lim sup X n agree and are uniformly bounded, then the sequence {X n } converges in the Hausdorff metric. When the density result holds, then the lim inf Z(T n ) must agree with lim sup Z(T n ). Hence, we have the following: Corollary 6. As compact subsets of ψ (S), Z(T n ) converges to the unimodular curve C in the Hausdorff metric. 22

23 Although we can determine the zero attractor and the zero density completely in the above framework. it is conceptually useful to have the result of Sokal that gives a description of the support of the zero density measure. [Sokal] [6]: Let D be a domain in C, and let z 0 D. Let {g n } be analytic functions on D, and let {a n } be positive real constants such that { g n an } are uniformly bounded on the compact subsets of D. Suppose that there does not exist a neighborhood V of z 0 and a function v on V that is either harmonic or else identically such that lim inf a n ln g n (z) v(z) lim sup a n ln g n (z) n n for all z V. Then z 0 lim inf Z(g n ). Remark: We can state the asymptotic form for T n (x) in a more symmetric form as: N T n (x) = ψ 0 (x) + a n,k (x)ψ k (x) cn + e n (x) k= where N is fixed and the error term has the form e n (x) = o(max{ψ k (x) cn, 0 k n}) This version explains the asymmetry in our first result where we have ψ 0 (x) = and the zeros accumulate along the curve ψ 0 (x) = ψ(x). A.5. Special Case. Theorem 6 shows that the images of zeros under the conformal map ψ are uniformly distributed along the corresponding circular arc. It can be applied to many cases that arise in a broad spectrum. A special case is worthy of attention; namely, the analytic arc C is a straight line segment and ψ(x) (see Lemma 5) has the form e ax+b, where a and b are constants. Corollary 7. If the analytic arc C is a straight line segment and ψ(x) is of the form e ax+b, where a and b are constants, then the zero density along the line segment C is a multiple of Lebesgue measure. Proof. Let c and c 2 be the endpoints of the line segment so it is parametrized as x(t) = c ( t)+c 2 t, 0 t. Then ψ(x(t)) becomes ψ(x(t)) = e a[c ( t)+c 2 t]+b Since ψ(x(t)) =, ψ(x(t)) can be written as ψ(x(t)) = e 2πiθ(t), where θ(t) is a linear function of t. By Theorem 6 the density of images of zeros under ψ(x) along the corresponding circular arc is Lebesgue measure. Hence, its pull-back under ψ is also Lebesgue measure since θ(t) is linear. References [] Pavel Bleher and Robert Mallison,Jr., Zeros of Sections of Exponential Sums, International Journal of Mathematics 2006 [2] Robert P. Boyer and William M.Y. Goh, On the Zero Attractor of the Euler Polynomials, Advances in Applied Math [3] K. Chandrasekharan, Arithmetic Functions, Springer-Verlag, 970. [4] E. T. Copson, Asymptotic Expansions, Cambridge University Press 965. [5] Peter Duren, Univalent Functions, Springer-Verlag, 983. [6] Alan Sokal, Chromatic roots are dense in the whole complex plane, Combin. Probab. Computing. 3 (2004), [7] G. Szegö, Über eine Eigenschaft der Exponentialreihe,Sitzungsber. Berl.Math.Ges.,23 (924), R. Boyer: Department of Mathematics, Drexel University, Philadelphia, PA address: rboyer@drexel.edu W. M. Y. Goh: Department of Statistics and Finance, University of Science and Technology of China (USTC), Hefei ,China address: wmygoh@hotmail.com 23