Convergent and divergent series, solutions of the Prolate - PowerPoint PPT Presentation

Convergent and divergent series, solutions of the Prolate Spheroidal differential equation Franc ¸oise Richard-Jung (joint work with F. Fauvet, J.P . Ramis, J. Thomann) LJK, BP 53, 38041 Grenoble Cedex, France e-mail address: Francoise.Jung@imag.fr Convergent and divergent series, solutions of the Prolate Spheroidal differential equation – p.1/26

Context: the D ESIR package (1) A linear homogeneous ODE L ( y ) = a d y ( d ) + a d − 1 y ( d − 1) + · · · + a 0 y = 0 in the neighborhood of the singularities: a d ( t ) = 0 Formal algorithms Series expansions y ( x ) = exp ( Q (1 /x )) x λ ˆ ϕ ( x ) , ˆ ϕ ( x ) ∈ C [[ x ]][ln x ] , t = x n , n ∈ N ∗ Convergent and divergent series, solutions of the Prolate Spheroidal differential equation – p.2/26

Context: the D ESIR package (2) Series expansions (with known Gevrey properties) Numerical algorithms: Stokes Summation of divergent series matrices Graphical visualization Convergent and divergent series, solutions of the Prolate Spheroidal differential equation – p.3/26

Outline • The prolate spheroidal wave functions (definitions and basic properties) • Experimental results on Stokes phenomenon and monodromy • Proof of the conjecture Convergent and divergent series, solutions of the Prolate Spheroidal differential equation – p.4/26

The prolate spheroidal wave functions The prolate spheroidal wave functions, (PSWFs) { ϕ n,σ,τ } , constitute an orthonormal basis of the space of σ -bandlimited functions on the real line, i. e. functions whose Fourier transforms have support on the interval [ − σ, σ ] . The PSWFs are maximally concentrated on the interval [ − τ, τ ] (as precised below) and depend on parameters σ and τ . Convergent and divergent series, solutions of the Prolate Spheroidal differential equation – p.5/26

Properties They can be characterized by one of the following properties: • as the maximum energy concentration of a σ -bandlimited function on the interval [ − τ, τ ] ; that is, ϕ 0 ,σ,τ is the function of total energy 1 ( = || ϕ 0 ,σ,τ || 2 ) such that � τ | f ( t ) | 2 − τ is maximized, ϕ 1 ,σ,τ is the function with the maximum energy concentration among those functions orthogonal to ϕ 0 ,σ,τ , etc...; Convergent and divergent series, solutions of the Prolate Spheroidal differential equation – p.6/26

Properties • as the eigenfunctions of an integral operator with kernel arising from the sinc function S ( t ) = sin( πt ) /πt : � τ σ � σ � ϕ n,σ,τ S π ( t − x ) dx = λ n,σ,τ ϕ n,σ,τ ( t ); π − τ Convergent and divergent series, solutions of the Prolate Spheroidal differential equation – p.7/26

Properties • as the eigenfunctions of an integral operator with kernel arising from the sinc function S ( t ) = sin( πt ) /πt : � τ σ � σ � ϕ n,σ,τ S π ( t − x ) dx = λ n,σ,τ ϕ n,σ,τ ( t ); π − τ • as the eigenfunctions of a differential operator arising from a Helmholtz equation on a prolate spheroid: ( τ 2 − t 2 ) ϕ ′′ n,σ,τ − σ 2 t 2 ϕ n,σ,τ = µ n,σ,τ ϕ n,σ,τ . n,σ,τ − 2 tϕ ′ Convergent and divergent series, solutions of the Prolate Spheroidal differential equation – p.7/26

Expansion of the solutions in the complex plane (1) L σ,τ = ( τ 2 − t 2 ) d 2 dt 2 − 2 t d dt − σ 2 t 2 Solutions of the equation L σ,τ ( ϕ ) − µϕ = 0 . Two finite singularities: τ, − τ . Both singularities are regular one. In the neighborhood of τ , a basis of solutions constituted of: • a regular function f • a solution of the form f log( t − τ ) + g , g is also regular. Convergent and divergent series, solutions of the Prolate Spheroidal differential equation – p.8/26

Expansion of the solutions in the complex plane (2) The point at infinity is also a singularity. This is an irregular one. We obtain the following basis of solutions, with x = 1 t : Iσ ( µ − σ 2 τ 2 ) x 1 − 1 � � y 1 ( x ) = e − Iσ x x + O( x 2 ) σ 2 2 Iσ ( µ − σ 2 τ 2 ) x 1 + 1 � � Iσ x x + O( x 2 ) y 2 ( x ) = e σ 2 2 There is no ramification. The series are a priori divergent, but 1 -summable in each direction but ± Iσ R + . Convergent and divergent series, solutions of the Prolate Spheroidal differential equation – p.9/26

Our goal • to have a more precise idea on the way these series diverge, depending on the parameter µ ; Convergent and divergent series, solutions of the Prolate Spheroidal differential equation – p.10/26

Our goal • to have a more precise idea on the way these series diverge, depending on the parameter µ ; • to compute the Stokes multipliers of the equation in the previous basis... with the hope that they are null for the known values of the eigenvalues. Convergent and divergent series, solutions of the Prolate Spheroidal differential equation – p.10/26

Our goal • to have a more precise idea on the way these series diverge, depending on the parameter µ ; • to compute the Stokes multipliers of the equation in the previous basis... with the hope that they are null for the known values of the eigenvalues. From now on: → ( t 2 − 1) y ′′ + 2 ty ′ + ( t 2 − µ ) y . σ = τ = 1 and D µ : y �− Convergent and divergent series, solutions of the Prolate Spheroidal differential equation – p.10/26

Neighborhood of 1 f µ ( t ) , f µ ( t )ln( t − 1) + ϕ µ ( t ) , where f µ ( t ) = 1+( − 1 2+ µ 2) ( t − 1)+( 1 16 µ 2 − 1 4 µ − 1 16) ( t − 1) 2 +O(( t − 1) 3 ) , and ϕ µ ( t ) = (1 2 − µ ) ( t − 1)+(1 2 µ + 1 16 − 3 16 µ 2 ) ( t − 1) 2 +O(( t − 1) 3 ) . Convergent and divergent series, solutions of the Prolate Spheroidal differential equation – p.11/26

Neighborhood of − 1 g µ ( t ) , g µ ( t )ln( t + 1) + ψ µ ( t ) , where g µ ( t ) = 1+(1 2 − µ 2) ( t +1)+( 1 16 µ 2 − 1 4 µ − 1 16) ( t +1) 2 +O(( t +1) 3 ) , and ψ µ ( t ) = ( − 1 2+ µ ) ( t +1)+(1 2 µ + 1 16 − 3 16 µ 2 ) ( t +1) 2 +O(( t +1) 3 ) All the series f µ , ϕ µ , g µ , ψ µ are convergent, with a radius of convergence at least 2 . Convergent and divergent series, solutions of the Prolate Spheroidal differential equation – p.12/26

Neighborhood of ∞ 1 − 1 2 I ( − 1 + µ ) x + ( − 1 8 µ 2 + 1 2 µ + 1 � − I x x 8) x 2 y 1 ( x ) = e ˆ � +O( x 3 ) 1 + 1 2 I ( − 1 + µ ) x + ( − 1 8 µ 2 + 1 2 µ + 1 � I 8) x 2 y 2 ( x ) = e ˆ x x � +O( x 3 ) Convergent and divergent series, solutions of the Prolate Spheroidal differential equation – p.13/26

Stokes phenomenon (1) − I I Let ˆ y 1 and ˆ f 1 = e f 2 = e x x y 2 . x ˆ x ˆ The series ˆ f 1 is 1 -summable in all directions but − I R + . The series ˆ f 2 is 1 -summable in all directions but I R + . θ ∈ ] − π 2 , π � − ζ t B ˆ f − i ( x ) = e f i ( ζ ) dζ, 2 [ d θ θ ∈ ] π 2 , 3 π � − ζ t B ˆ f + i ( x ) = e f i ( ζ ) dζ, 2 [ d θ 1 = f + f − 1 Convergent and divergent series, solutions of the Prolate Spheroidal differential equation – p.14/26

Stokes phenomenon (2) π 0 2 π − π θ θ Domain of y + Domain of y − i i Convergent and divergent series, solutions of the Prolate Spheroidal differential equation – p.15/26

Stokes phenomenon (3) i and y + The functions y − i are defined together on the sector ] − π, π [ . The Stokes matrix associated to the Stokes ray π 2 is: π 2 ) = ( y + y + 2 . ( y − y − 2 ) S 1 1 � � 1 α π 1 = y + 2 = As y − 1 , S . 0 1 For similar reasons: � � 1 0 − π 2 = S . α 1 Convergent and divergent series, solutions of the Prolate Spheroidal differential equation – p.16/26

Stokes matrices, in practice It is possible to compute: 1. a linear ODE satisfied by ˆ f 1 (formal) 2. a linear ODE satisfied by B ˆ f 1 (formal) 3. a basis of formal solutions of the previous ODE near each singularity ω , on the example ω = − 2 I (formal) 4. the expression of B ˆ f 1 in this basis (numerical) We deduce a numerical value of α . Convergent and divergent series, solutions of the Prolate Spheroidal differential equation – p.17/26

Stokes matrices: for which equations ? • of single rank, k � = 1 ; • the Borel transform of the divergent series can have any polar, ramified or logarithmic singularities; • the Borel transform of the divergent series don’t have irregular singularities; • the Borel transform of the divergent series don’t have many singularities aligned on a half line issued from the origin. Convergent and divergent series, solutions of the Prolate Spheroidal differential equation – p.18/26

A first Stokes matrix For the first eigenvalue µ 0 = 0 . 319 : > StokesMatrices(subs(mu=319/1000,resu), [8 , 8]); � � [[ − π 1 0 2 , ] , 0 . 1095087321 10 − 23 − 0 . 4342357172 10 − 6 I 1 − 0 . 1095087321 10 − 23 − 0 . 4342357172 10 − 6 I � � [ π 1 2 , ]] 0 1 Convergent and divergent series, solutions of the Prolate Spheroidal differential equation – p.19/26

A curve α ( µ ) Numerically : α ∈ I R . We draw the curve ℑ ( α ( µ )) , µ ∈ ]0 , 7] . 2 1 0 1 2 3 4 5 6 7 –1 –2 µ 0 = 0 . 319 , µ 1 = 2 . 593084 , µ 2 = 6 . 533471 . Convergent and divergent series, solutions of the Prolate Spheroidal differential equation – p.20/26