arxiv: v2 [nlin.si] 4 Oct 2017

Transcription

1 ASYMPTOTIC BEHAVIOUR OF THE FIFTH PAINLEVÉ TRANSCENDENTS IN THE SPACE OF INITIAL VALUES NALINI JOSHI AND MILENA RADNOVIĆ arxiv: v [nlinsi] 4 Oct 017 Abstract We study the asymptotic behaviour of the solutions of the fifth Painlevé equation as the independent variable approaches zero and infinity in the space of initial values We show that the limit set of each solution is compact and connected and, moreover, that any solution with an essential singularity at zero has an infinite number of poles and zeroes, and any solution with an essential singularity at infinity has infinite number of poles and, moreover, takes the value unity infinitely many times 1 Introduction Following the construction of initial value spaces of Painlevé equations by Okamoto (1979), and more recent asymptotic analysis of the solutions of the first, second and fourth Painlevé equations in such spaces, we investigate the solutions of the fifth Painlevé equation: ( d y 1 P V : dx = y + 1 )( ) dy 1 ( dy (y 1) + αy + β ) y 1 dx xdx x y (11) + γy δy(y +1) + x y 1, in an asymptotic limit in its initial value space Complete information about the limit sets of transcendental solutions and their behaviours near the infinity set are found Unlike earlier asymptotic investigations of P V, we do not impose any reality constraints; here y is a function of the complex variable x, and α, β, γ, δ are given complex constants Noting that P V has (fixed) essential singularities only when the independent variable x takes the values 0 and, we investigate the behaviour of the solutions near x = 0 A similar analysis can be carried out near and we also include an outline of the main results for this limit We show that each solution that is singular at x = 0 has infinitely many poles and zeroes in every neighbourhood of this point Similarly, each solution singular at x = has infinitely many poles and, moreover, takes the value 1 infinitely many times in each neighbourhood of infinity The starting point for our analysis is the compactification and regularisation of the initial-value space To make explicit analytic estimates possible, we calculate detailed information about the Painlevé vector field after each resolution (or blow-up) of this space A similar approach was carried out for the first, second, and fourth Painlevé equations respectively by Duistermaat and Joshi (011), Howes and Joshi (014), and Joshi and Radnović (016) However, the construction of the initial-value spaces in each This research was supported by an Australian Laureate Fellowship# FL from the Australian Research Council The research of MR was partially supported by by the Serbian Ministry of Education and Science (Project no 17400: Geometry and Topology of Manifolds and Integrable Dynamical Systems) 1

2 NALINI JOSHI AND MILENA RADNOVIĆ of these earlier works consisted of exactly nine blow-ups, while in the present paper, we will see that eleven blow-ups are needed, followed by two blow-downs The initial-value space is then obtained by removing the set, denoted by I, of points which are not attained by any solution The main results obtained in this paper fall into four parts: Existence of a repeller set: Theorem 48 in Section 4 shows that I is a repeller for the flow The theorem also provides the range of the independent variable for which a solution may remain in the vicinity of I Numbers of poles and zeroes: In Corollary 49, we prove that each solution that is sufficiently close to I has a pole in a neighbourhood of the corresponding value of the independent variable Moreover, Theorem 55 shows that each solution with essential singularity at x = 0 has infinitely many poles and infinitely many zeroes in each neighbourhood of that point The complex limit set: We prove in Theorem 5 that the limit set for each solution is non-empty, compact, connected, and invariant under the flow of the autonomous equation obtained as x 0 Asymptotic behaviour as x : We show in Section 6 that each solution with an essential singularity at x = has infinitely many poles and takes the value unity infinitely many times in each neighbourhood of that point The asymptotic analysis of the fifth Painlevé transcendent has been studied by many authors, including Andreev and Kitaev(1997a, 1997b, 000), Zeng and Zhao (016), Bryuno and Parus and Qin and Shang (006), Lu and McLeod (1999a, 1999b), McCoy and Tang (1986a, 1986b, 1986c), and Jimbo (198) However, the literature on the asymptotic behaviours of the fifth Painlevé transcendent concentrates on behaviours on the real line, often focusing on special behaviours or solutions, while we consider all solution behaviours for x C For other mathematical results related to P V, see (Boelen et al, 013; Shimomura, 011; Kaneko and Ohyama, 007; Sasaki, 007; Clarkson, 005; Lu and Shao, 004; Gordoa et al, 001a; 001b), while for applications in physics see(jimbo et al, 1980; Dyson, 1995; Schief, 1994), and references therein This paper is organised as follows In Section, we construct and describe the space of the initial values for equation (11), with complete details of all the necessary calculations provided in Appendix A In Section 3, we consider the special solutions of P V Section 4 contains the analysis of the behaviours of solutions near the infinity set in the space of initial values Results on the complex limit sets of solutions when the independent variable approaches 0 are provided in Section 5 The behaviours of the fifth Painlevé transcendent in the limit x are outlined in Section 6 A summary of the notation used in this paper is given in Appendix B The space of initial values Since the fifth Painlevé equation is a second-order ordinary differential equation, solutions are locally defined by two initial values Therefore, the space of initial values is two complex-dimensional However, standard existence and uniqueness theorems only cover values of y that are not arbitrarily close to 0, 1 or infinity (where the second derivative given by P V becomes ill-defined) In this section, we explain how to construct a regularized, compactified space of all possible initial values that overcomes these issues

3 ASYMPTOTICS OF P V 3 We start by formulating P V as an equivalent system of equations in Section 1 and describing its autonomous limiting form obtained as x 0 in Section The mathematical construction of the space of initial values is then given in Section 3 Where y is arbitrarily close to the singular values 0, 1,, the solutions have singular power series expansions, which become regularized Taylor expansions in corresponding domains of the initial value space These regular expansions are provided in Section 4 1 A system equivalent to P V Withthechangeoftheindependentvariablet = logx, Equation (11) becomes: ( d y 1 dt = y + 1 )( ) ( dy +(y 1) αy + β ) (1) +e t γy +e tδy(y+1) y 1 dt y y 1 We rewrite Equation (1) in the following way: () dy dt =y(y 1) z (θ 0 +η)y +(θ 0 +η θ 1 e t )y θ 0, dz dt = (y 1)(3y 1)z + ( (θ 0 +η)y θ 0 η +θ 1 e t) z 1 ǫ(θ 0 +η θ ), where θ = α, θ 0 = β, θ 1 = δ (θ 1 0), η = γ θ 1 1, and ǫ = 1 (θ 0 +θ +η) Remark 1 Here, we assumed that δ 0, which is a generic case of the fifth Painlevé equation Whenδ = 0, the P V is equivalentto the third Painlevéequation (Ohyama and Okumura, 006) The system () is Hamiltonian: with Hamiltonian function: dy dt = H z, dz dt = H y, (3) H = y(y 1) z (θ 0 +η)y z +(θ 0 +η θ 1 e t )yz θ 0 z + 1 ǫ(θ 0 +η θ )y The autonomous equation The autonomous equation corresponding to (1) is: ( d y 1 (4) dt = y + 1 )( ) ( dy +(y 1) αy + β ), y 1 dt y which is equivalent to the autonomous version of (): dy (5) dt = (y 1) (yz θ 0 )+ηy(y 1), dz dt = (y 1)z(η +θ 0 +z 3yz)+ηz 1 ǫ(θ 0 +η θ ) System (5) is also Hamiltonian: with Hamiltonian: dy dt = E z, dz dt = E y, (6) E = y(y 1) z (θ 0 +η)y z +(θ 0 +η)yz θ 0 z + 1 ǫ(θ 0 +η θ )y

4 4 NALINI JOSHI AND MILENA RADNOVIĆ Using the first equation of (5) to express z, and using the fact that E is constant along solutions, we get: ( ) dy = (y 1) (4Cy +θ0 dt θ 0(η +θ 0 )y +θ y ), C = const It is worth observing that the constant function y 1 is the only solution of this equation taking the value 1 From (5), the corresponding function z is the solution of That is, we have dz dt = ηz 1 ǫ(θ 0 +η θ ) z = c 1 e ηt + ǫ(θ 0 +η θ ), η where c 1 is a constant The flow (5) has four fixed points: ( 1, ǫ(θ ) 0 +η θ ) η (y,z) = ( Y i, θ ) 0 η +,i {1,,3}, Y i (1 Y i ) where Y 1, Y, Y 3 are the roots of the following cubic polynomial in Y: ( ) θ 6η +8ηθ 0 +θ0 Y 3 + ( ) η 1ηθ 0 θ0 +θ Y +θ 0 (4η θ 0 )Y +θ0 3 Resolution of singularities In this section, we explain how to construct the space of initial values for the system () First, we motivate the reason for this construction before introducing the notion of initial value spaces in Definition 3, which is based on foliation theory We then explain how to construct such a space by carrying out resolutions or blow-ups, which are described in Definition 4 The system () is a system of two first-order ordinary differential equations for (y(t),z(t)) Given initial values (y 0,z 0 ) at t 0, local existence and uniqueness theorems provide a solution that is defined on a local polydisk U V in C C, where t 0 U C and (y 0,z 0 ) V C Our interest lies in global extensions of these local solutions However, the occurence of movable poles in the Painlevé transcendents acts as a barrier to the extension of U V to the whole of C C The first step to overcome this obstruction is to compactify the space C, in order to include the poles We carry this out by embedding C into CP and explicitly represent the system () in the three affine coordinate charts of CP, which are given in Sections A1 A3 The second step in this process results from the occurence of singularities in the Painlevé vector field () These occur in the (y 0,z 0 ) and (y 03,z 03 ) charts The appearence of these singularities is related to irreducibility of the solutions of Painlevé equations as indicated by the following theorem, due to Painlevé Theorem ((Painlevé, 1897)) If the space of initial values for a differential equation is a compact rational variety, then the equation can be reduced either to a linear differential equation of higher order or to an equation for elliptic functions

5 ASYMPTOTICS OF P V 5 It is well known that the solutions of Painlevé equations are irreducible (in the sense of the theorem) Since CP is a compact rational variety, the theorem implies CP cannot be the space of initial values for () By the term singularity we mean points where(dy/dt, dz/dt) becomes either unbounded or undefined because at least one component approaches the undefined limit 0/0 We are led therefore to construct a space in which the points where the singularities are regularised The process of regularisation is called blowing up or resolving a singularity We now define the notion of initial value space Definition 3 [(Gérard, 1976), (Gérard and Sec, 197; Gérard, 1983; Okamoto, 1979)] Let (E,π,B) be a complex analytic fibration, Φ its foliation, and a a holomorphic differential system on E, such that: the leaves of Φ correspond to the solutions of ; the leaves of Φ are transversal to the fibres of E; for each path p in the base B and each point X E, such that π(x) p, the path p can be lifted into the leaf of Φ containing point X Then each fibre of the fibration is called a space of initial values for the system The properties listed in Definition 3 imply that each leaf of the foliation is isomorphic to the base B Since the fifth Painlevé transcendents (in the t variable) can be globally extended as meromorphic functions on C(Joshi and Kruskal, 1994; Hinkkanen and Laine, 001), we are searching for the fibration with the base equal to C In order to construct the fibration, we apply the blow-up procedure, defined below, (Hartshorne, 1977; Griffiths and Harris, 1978; Duistermaat, 010) to the singularities of the system () that occur where at least one component becomes undefined of the form 0/0 Okamoto (1979) showed that such singular points are contained in the closure of infinitely many leaves Moreover, these leaves are holomorphically extended at such a point Definition 4 The blow-up of the plane C at point (0,0) is the closed subset X of C CP 1 defined by the equation u 1 t = u t 1, where (u 1,u ) C and [t 1 : t ] CP 1, see Figure 1 There is a natural morphism ϕ : X C, which is the restriction of the projection from C CP 1 to the first factor ϕ 1 (0,0) is the projective line {(0,0)} CP 1, called the exceptional line Remark 5 Notice that the points of the exceptional line ϕ 1 (0,0) are in bijective correspondence with the lines containing (0,0) On the other hand, ϕ is an isomorphism between X\ϕ 1 (0,0) and C \{(0,0)} More generally, any complex two-dimensionalsurface can be blown up at a point(hartshorne, 1977; Griffiths and Harris, 1978; Duistermaat, 010) In a local chart around that point, the construction will look the same as described for the case of the plane Notice that the blow-up construction separates the lines containing the point (0, 0) in Definition 4, as shown in Figure 1 In this way, the solutions of () containing the same point can be separated Additional blow-ups may be required if the solutions have a commont tangent line or a tangency of higher order at such a point The explicit resolution of the vector field () is carried out in Appendix A As mentioned above, the process requires 11 resolutions of singularities, or, blow-ups

6 6 NALINI JOSHI AND MILENA RADNOVIĆ ϕ X Figure 1 The blow-up of the plane at a point L -1 L -3 L 3 - L 8 - L 6-1 L 5-1 L 0 - L 7-1 L 4 - L 10 (t) -1 L 1-4 L 9 (t) - Figure The fibre D(t) is obtained from CP by 11 blow-ups The resulting surface D(t) is shown in Figure Note that in this figure, L is the proper preimage of the line at the infinity, while L 0, L 1, L, L 3, L 4, L 8, L 9(t) are proper preimages of exceptional lines obtained by blow ups at points a 0, a 1, a, a 3, a 4, a 8, a 9 respectively and L 5, L 6, L 7, L 10 (t) are exceptional lines obtained by blowing up points a 5, a 6, a 7, a 10 respectively The self-intersection number of each line (after all blow-ups are completed) is indicated in the figure Okamoto described so called singular points of the first class that are not contained in the closure of any leaf of the foliation given by the system of differential equations At such points, the corresponding vector field is infinite (For example, in chart (y 0,z 0 ) from Section A such a singular point is given by y 0 = 0 with non-zero z 0 ) In the surface D(t), all remaining singular points are of the first class, and the fibre of the initial

7 ASYMPTOTICS OF P V 7 value space is obtained by removing them: ( 4 ) E(t) = D(t)\ L j L 8 L 9 L j=0 In D(t), each line with self-intersection number 1 can be blown down again Blowing down L and then the projection of L 0, we get the surface F(t), which is shown in Figure 3 The projection of each remaining line from D(t) is denoted by the same index but now with superscript p Notice that the self-intersection numbers of L p 1 and L p are no longer the same as of the corresponding pre-images L 1 and L In this space, we denote by I the set of all singular points of the first class in F(t): 4 I = L p j Lp 8 L p 9 The fibre E(t) of the initial value space can be identified with F(t)\I j=1 L p - L p 3 - L p 8 - L p 6-1 L p 5-1 L p 7-1 L p 4 - L 10 (t) p -1 L p 1 - L 9 (t) p - Figure 3 The fibre F(t) is obtained from CP by 11 blow-ups and two blow-downs If S is a surface obtained from the projective plane by a several successive blow-ups of points, then the group of all automorphisms of the Picard group Pic(S) preserving the canonical divisor K is generated by the reflections X X +(Xω)ω, with ω Pic(S), Kω = 0, ωω = That group is an affine Weyl group and the lines of self-intersection are its simple roots Representing each such line by a node and connecting a pair of nodes by a line only if they intersect in the fibre, we obtain the Dynkin diagram shown in Figure 4, which is of type D (1) 5 For detailed expositions on the topic of surfaces and root systems, see (Demazure, 1976; Harbourne, 1985) and references therein In the limit Ret, the resulting Okamoto space is compactified by the fibre F, corresponding to the autonomous system (5), see Figure 5 Its infinity set is given by I = L a 1 Lp La 3 F is constructed by a sequence of blow-ups, carried out along the lines of the construction of the space of the initial values for the non-autonomous system () For those readers who may wish to carry this out separately, we note that the resolution at point a 4 (found in Section A5) needs to be done only for the non-autonomous system

8 8 NALINI JOSHI AND MILENA RADNOVIĆ L p L p 4 L p 1 L p 8 L p 3 L p 9 Figure 4 The Dynkin diagram of D (1) 5 L a - L a 3 - L a 6-1 L a 5-1 L a 7-1 L a 1-1 Figure 5 The fibre F corresponding to the autonomous system (5) 4 Movable singularities in Okamoto s space Here, we consider neighbourhoods of exceptional lines where the Painlevé vector field () becomes unbounded The construction given in Appendix A shows that these are given by the lines L 5, L 6, L 7 and L 10 Movable pole with residue θ The set L 5 \I is given by y 6 = 0 in the (y 6,z 6 ) chart, see Section A9 Suppose y 6 (τ) = 0, z 6 (τ) = B, for arbitrary complex numbers τ, B Solving the system of differential equations for y 6, z 6 in Section A9 in a neighbourhood of t = τ, we obtain y 6 = θ (t τ)+ θ ( B θ η θ 1 e τ) (t τ) +O((t τ) 3 ) Since the transformation from (y 6,z 6 ) to (y,z) (also given in Section A9) is y = 1 y 6 and z = y 6 (y 6 z 6 +1), we obtain series expansions for (y,z) given by y = θ t τ B θ η θ 1 e τ θ +O(t τ), z = θ (t τ)+ θ (B(1+θ ) η θ θ 1 e τ )(t τ) +O((t τ) 3 ) Clearly, y has a simple pole with residue θ, while z has a simple zero at t = τ

9 ASYMPTOTICS OF P V 9 Movable pole with residue θ At the intersection with L 6 \I, y has a simple pole with the residue θ, while z has a simple zero This case is analogous to the previous one, see Sections A9 and A10 Movable zero with coefficient θ 0 The set L 7 \I is given by z 81 = 0 in the (y 81,z 81 ) chart, see Section A11 Suppose y 81 (τ) = B and z 81 (τ) = 0 Then integration of the vector field gives z 81 (t) = (t τ)+ η θ 0 θ 1 e τ with F = 1 ǫ(θ 0 +η θ ) Since we obtain (t τ) + F ηθ 0 +θ0 θ 1e τ 4B (t τ) 3 +O((t τ) 4 ), 3 y = z 81 (y 81 z 81 +θ 0 ), z = 1 z 81, y = θ 0 (t τ)+ B +θ 0(η θ 0 θ 1 e τ ) (t τ) +O((t τ) 3 ), z = 1 t τ η θ ( 0 θ 1 e τ (η θ0 θ 1 e τ ) + 4 +O((t τ) ) Thus, y has a simple zero and z a simple pole at t = τ F ηθ 0 +θ 0 θ 1 e τ 4B 3 ) (t τ) Movable points where y becomes unity ThesetL 10 \I isgivenbyz 111 = 0inthe(y 111,z 111 ) chart, see Section A14 Suppose y 111 (τ) = B and z 111 (τ) = 0 Then z 111 = (t τ)+ θ 1e τ η 1 (t τ) Since we have ( 1+ 5 η +η B θ 1 e τ +θ 1e τ (1 θ 0 η y = 1+θ 1 e t z 111 +(1+η)θ 1 e t z111 +y 111z111 3, 1 z = z111(θ 1 e t +(1+η)θ 1 e t z 111 +y 111 z111), )+ 1 θ1 eτ ) (t τ) 3 +O((t τ) 4 ) it follows that y(τ) = 1 and z has a double pole at t = τ Their expansions around this point are given by y = 1+θ 1 e τ (t τ)+ θ 1e τ (θ 1 e τ +η +3) (t τ) + θ ( 1e τ 4η 4η B 6 z = (θ 1e τ ) 1 (t τ) + 1+(θ 1e τ ) 1 + t τ ) θ 1 e +(11+5η θ 0)θ τ 1 e τ +θ1 eτ (t τ) 3 +O((t τ) 4 ), ( B θ 1 e + ηθ ) 1e τ +θ τ 0 θ 1 e τ + θ 1 eτ 1+3η +O(1) 4 4

10 10 NALINI JOSHI AND MILENA RADNOVIĆ 3 Special solutions In this section, we consider the pencil of curves corresponding to the Hamiltonian (6) of the autonomous system (5) We show that by a birational equivalence, this pencil can be transformed to a pencil of conics In Section 31, we analyse the solutions of the fifth Painlevé equation corresponding to singular conics from the pencil In Section 3, we give a summary of other special solutions 31 Special solutions and singular conics The pencil of curves arising from the Hamiltonian of the autonomous system (5) is given by the zero set of the one-parameter family of polynomials (parametrised by c): (31) h c (y,z) = y(y 1) z (θ 0 +η)y z +(θ 0 +η)yz θ 0 z + 1 ǫ(θ 0 +η θ )y c For each c, the curve h c (y,z) = 0 is birationally equivalent to h 1 c (y 1,z 1 ) = z 1 (θ 0 + cy 1 +θ y 1 ) 4y 1 = 0, where c = 4c θ 0 (η +θ 0 ) and the birational equivalence is given by y 1 = y, z 1 = y(y 1)z (θ 0 +η)y +θ 0 The family of level curves h1 c = 0 forms a pencil of conics The singular conics in the pencil are: c = θ 0 θ, (z 1 +θ y 1 +θ 0 )(z 1 θ y 1 θ 0 ) = 0, c = θ 0 θ, (z 1 θ y 1 +θ 0 )(z 1 +θ y 1 θ 0 ) = 0, c =, y 1 w 1 = 0, where [y 1 : z 1 : w 1 ] are homogeneous coordinates Moreover, the base points of this pencil are: [0 : θ 0 : 1], [0 : θ 0 : 1], [1 : θ : 0], [1 : θ : 0] Recalling that the Hamiltonian E of the autonomous system, given by Equation (6), has the time derivative E = θ 1 e t( z(θ 0 (η +θ 0 )y +y(y 1)z)+E ), we can also search for conditions under which all successive derivatives of E are zero This happens if and only if z = 0 or θ 0 (η +θ 0 )y +y(y 1)z = 0 These cases are investigated in further detail below Casez = 0 Fromthesecondequationof()wehaveǫ(θ 0 +η θ ) = 0,ieθ 0 +η = ±θ In this case z 1 = ±θ y 1 +θ 0, and so this case corresponds to lines in the pencil of conics containing the base point [0 : θ 0 : 1] The first equation of () is then a Riccati equation: which is equivalent to dy dt = θ y +(θ 0 ±θ θ 1 e t )y θ 0, (3) x dy dx = θ y +(θ 0 ±θ θ 1 x)y θ 0

11 ASYMPTOTICS OF P V 11 For θ θ 1 0, the solutions of this equation can be expressed in terms of the Whittaker functions (Olver et al, 010) Note that θ 1 0 is equivalent to δ 0, when the P V can be renormalised to δ = 1/ Then, the solutions of (3) are given by where y = zφ (x) θ φ(x), φ(x) = C 1M κ,µ (ξ)+c W κ,µ (ξ) ξ κ e ξ/, with ξ = ±x, κ = ( θ θ 0 +1)/, µ = θ +θ 0, and C 1, C being arbitrary constants Solutions from this class intersect only the pole lines L 5 and L 6 Case θ 0 (η +θ 0 )y +y(y 1)z = 0 The first equation of () is then: dy dt = θ 0y (θ 0 +θ 1 e t )y +θ 0, which is again a Riccati equation that can be solved analogously to the previous case of (3) Since in this case z 1 = θ 0 y θ 0, the condition on the constants is θ 0 = ±θ The solutions belong to the lines from the pencil of conics that contain the base point [0 : θ 0 : 1] Solutions from this class intersect only the pole line L 7 3 Other special solutions of P V In this section, we give a brief account on special solutions and Bäcklund transformations of P V For a detailed exposition, we refer readers to (Lukaševič, 1968; Airault, 1979; Okamoto, 1987; Gromak, 1999; Schief, 00; Forrester and Witte, 00; Masuda et al, 00; Clarkson, 005; Olver et al, 010) and the references therein All rational solutions of P V are of the form y = λx+µ+ P(x) Q(x), where λ, µ are constants and P, Q polynomials of degrees n 1 and n respectively Such solutions can be obtained by applying Bäcklund transformations: S 1 : (y(x);α,β,γ,δ) (y( x);α,β, γ,δ), ( ) 1 S : (y(x);α,β,γ,δ) y(x) ; β, α, γ,δ, ( T ε1,ε,ε 3 : y(x);α,β,γ, 1 ) ( Φ ε1 xy ;α,β,γ, 1 ), Φ where ε 1,ε,ε 3 { 1,1} and α = 1 ( )) ; (γ +ε 1 1 ε 3 β ε α 8 β = 1 ( )) ; (γ ε 1 1 ε 3 β ε α 8 ( ) γ = ε 1 ε 3 β ε α ; Φ = xy ε αy +ε 3 β +(ε α ε3 β +ε1 x)y

12 1 NALINI JOSHI AND MILENA RADNOVIĆ to so called seed solutions of P V, which are given by (33) y = κx+µ, for α = 1, β = 1 µ, γ = κ( µ), δ = 1 κ ; κ x+κ, for α = 1, β = κ µ, γ = κµ, δ = µ; κ+x κ x, for α = 1 8, β = 1 8, γ = κµ, δ = µ, where κ, µ are arbitrary constants In addition, P V has the following elementary (non-rational) solutions: (34) y = 1+κ x, for α = µ, β = 1 8, γ = µκ, δ = 0; y = κe µx, for α = β = 0, γ = µ, δ = 1 µ, where κ, µ are arbitrary constants For δ = 0, algebraic solutions are all rational functions in x and can be obtained by consecutive application of the Bäcklund transformations S 1 and S to the first function of (34) 4 The solutions near the infinity set In this section, we study the behaviour of the solutions of the system () near the set I, where the vector field is infinite We prove that I is a repeller for the solutions and that each solution which comes sufficiently close to I at a certain point t will have a pole in a neighbourhood of t Lemma 41 For every ǫ 1 > 0 there exists a neighbourhood U of L p E E +θ 1e t < ǫ 1 in U such that Proof In the corresponding charts near L p (see Section A6), the function is equal to: r = E E +θ 1e t θ 1 e t y 31 z 31 (ηy 31 (θ 0 y 31 1)(y 31 z 31 1)) r 31 = ηy 31 (y 31 z 31 1)+Fy31 θ 0y 31 (y 31 z 31 1) +y31 z 31 y 31z 31 +1, θ 1 e t y 3 z 3 (η (y 3 1)(θ 0 z 3 )) r 3 = (y 3 1)z 3 (η (y 3 1)(θ 0 z 3 ))+F The statement of the lemma follows immediately from these expressions, since L p is given by z 31 = 0 and y 3 = 0 Lemma 4 For each compact subset K of L p 1 \Lp 8 there exists a neighbourhood V of K and a constant C > 0 such that: e te E < C in V for all t

13 ASYMPTOTICS OF P V 13 Proof Near L p 1, in the respective coordinate charts (see Section A5), we have: θ 1(y 1 +1) e te E y 1 1, θ 1 (z +1) z 1 Since the projection of L p 8 to these two charts is the point on Lp 1 given by coordinates y 1 = 1 and z = 0, the statement is proved The modulus d of the function d from the next lemma will serve as a measure for the distance from the infinity set I Throughout the paper, we denote by J n the Jacobian of the coordinate change from (y,z) to (y n,z n ): J n = y n y z n z y n z n z y Lemma 43 There exists a continuous complex valued function d on a neighbourhood of the infinity set I in the Okamoto space, such that: 1, in a neighbourhood of I E \(Lp 3 L p 4 L p 8 L p 9), J 8, in a neighbourhood of L p 3 d = \Lp 1, J 10, in a neighbourhood of (L p 4 L p 8)\L p 1, J 11, in a neighbourhood of L p 9 \Lp 8 Proof From Section A11, the line L p 3 is given by z 8 = 0 in the (y 8,z 8 ) chart Thus as we approach L p 3, ie, as z 8 0, we have EJ 8 1 From Section A13, the line L p 8 is given by z 10 = 0 in the (y 10,z 10 ) chart Thus as we approach L p 8, ie, as z 10 0, we have EJ 10 1 θ 1e t y 10 In the same chart, the line L p 4 is given by y 10 = θ 1 e t From Section A14, the line L p 9 is given by z 11 = 0 in the (y 11,z 11 ) chart We have: J 11 J 10 = 1+ (1+η)θ 1e t y 11 Lemma 44 (Behaviour near L p 3 \ L p 1) If a solution at a complex time t is sufficiently close to L p 3 \Lp 1, then there exists unique τ C such that (y(τ),z(τ)) belongs to the line L 7 In other words, the solution has a pole at t = τ Moreover t τ = O( d(t) y 8 (t) ) for sufficiently small d(t) and bounded y 8 For large R 3 > 0, consider the set {t C y 8 R 3 } Its connected component containing τ is an approximate disk D 3 with centre τ and radius d(τ) R 3, and t y 8 (t) is a complex analytic diffeomorphism from that approximate disk onto {y C y R 3 } Proof Forthe study of the solutions near L p 3\L p 1, we use coordinates (y 8,z 8 ), see Section A11 In this chart, the line L p 3\L p 1 is given by the equation z 8 = 0 and parametrised by y 8 C Moreover, L p 7 is given by y 8 = 0 and parametrised by z 8 C

14 14 NALINI JOSHI AND MILENA RADNOVIĆ Asymptotically, for z 8 0 and bounded y 8, e t, we have: (41a) y 8 1, z 8 (41b) z 8 4y8z 3 8, (41c) J 8 = z 8, (41d) (41e) J 8 = η θ 0 θ 1 e t 4y 8 +O(z 8 ) = η θ 0 θ 1 e t 4y 8 +O(J 8 ), J 8 EJ 8 1 Integrating (41d) from τ to t, we get J 8 (t) = J 8 (τ)e K(t τ) e θ 1(e t e τ) (1+o(1)), K = η θ 0 4y 8 ( τ), where τ is on the integration path Because of (41b), z 8 is approximately equal to a small constant, and from (41a) follows that: y 8 y 8 (τ)+ t τ z 8 Thus, if t runs over an approximate disk D centered at τ with radius z 8 R, then y 8 fills and approximate disk centered at y 8 (τ) with radius R Therefore, if z 8 (τ) τ, the solution has the following properties for t D: z 8 (t) z 8 (τ) 1, and y 8 is a complex analytic diffeomorphism from D onto an approximate disk with centre y 8 (τ) and radius R If R is sufficiently large, we will have 0 y 8 (D), ie the solution of the Painlevé equation will have a pole at a unique point in D Now, it is possible to take τ to be the pole point For t τ τ, we have: d(t) d(τ) 1, ie z 8 (t) d(τ) J 8(t) d(τ) 1, y 8 (t) t τ z 8 t τ d(τ) Let R 3 be a large positive real number Then the equation y 8 (t) = R 3 corresponds to t τ d(τ) R 3, which is still small compared to τ if d(τ) is sufficiently small Denote by D 3 the connected component of the set of all t C such that {t y 8 (t) R 3 } is an approximate disk with centre τ and radius d(τ) R 3 More precisely, y 8 is a complex analytic diffeomorphism from D 3 onto {y C y R 3 }, and d(t) d(τ) 1 for all t D 3 We have E(t)J 8 (t) 1 when z 8 1 Thus E(t)J 8 (t) 1 for the annular disk t D 3 \ D 3, where D 3 is a disk centered at τ with small radius compared to radius of D 3

15 ASYMPTOTICS OF P V 15 Lemma 45 (Behaviour near L p 9 \ Lp 8 ) If a solution at a complex time t is sufficiently close to L p 9 \Lp 8, then there exists unique τ C such that (y(τ),z(τ)) belongs to the line L 10 In other words, the solution has a pole at t = τ Moreover t τ = O( d(t) y 11 (t) ) for sufficiently small d(t) and bounded y 11 For large R 9 > 0, consider the set {t C y 11 R 9 } Its connected component containing τ is an approximate disk D 9 with centre τ and radius d(τ) R 9, and t y 11 (t) is a complex analytic diffeomorphism from that approximate disk onto {y C y R 9 } Proof For the study of the solutions near L p 9 \ Lp 8, we use coordinates (y 11,z 11 ), see Section A14 In this chart, the line L p 9 \ Lp 8 is given by the equation z 11 = 0 and parametrised by y 11 C Moreover, L p 10 is given by y 11 = 0 and parametrised by z 11 C Asymptotically, for z 11 0 and bounded y 11, e t, we have: (4a) y 11 1, z 11 (4b) (4c) (4d) (4e) z 11 3(1+η)z 11 θ 1 e ty 11z 11, J 11 θ 1 e t z 11, J 11 = η +θ 1 e t +O(z 11 ) = η+θ 1 e t +O(J 11 ), J 11 J 11 = 1+ (1+η)θ 1e t J 10 y 11 Integrating (4d) from τ to t, we get J 11 (t) = J 11 (τ)e (+η)(t τ) e θ 1(e t e τ) (1+o(1)) Because of (4c), z 11 is approximately equal to a small constant, and from (4a) follows that: y 11 y 11 (τ)+ t τ z 11 Thus, if t runs over an approximate disk D centered at τ with radius z 11 R, then y 11 fills and approximate disk centered at y 11 (τ) with radius R Therefore, if z 11 (τ) τ, the solution has the following properties for t D: z 11 (t) z 11 (τ) 1, and y 11 is a complex analytic diffeomorphism from D onto an approximate disk with centre y 11 (τ) and radius R If R is sufficiently large, we will have 0 y 11 (D), ie the solution of the Painlevé equation will have a pole at a unique point in D Now, it is possible to take τ to be the pole point For t τ τ, we have: d(t) d(τ) 1, ie z 11 (t) 1 d(τ) θ 1 e J11(t) t d(τ) y 11 (t) t τ t τ z 11 d(τ) θ 1e t 1 θ 1 e t, Let R 9 be a large positive real number Then the equation y 11 (t) = R 9 corresponds to (t τ)θ 1 e t d(τ) R 3, which is still small compared to τ if d(τ) is sufficiently small

16 16 NALINI JOSHI AND MILENA RADNOVIĆ DenotebyD 9 theconnectedcomponentofthesetofallt Csuchthat{t y 11 (t) R 9 } is an approximate disk with centre τ and radius d(τ) R 3 More precisely, y 11 is a complex analytic diffeomorphism from D 9 onto {y C y R 9 }, and From (4e), we have: that is, when Since R 9 1, we have J 11 J 10 1 when 1 d(t) d(τ) 1 for all t D 9 (1+η)θ 1 e t y 11 (t) t τ d(τ) t τ d(τ) R 9 d(τ) (1+η)d(τ) t τ, Thus J 11 1 for the annular disk t D 9 \ D 9 J, where D 9 is a disk centered at τ with 10 small radius compared to the radius of D 9 Lemma 46 (Behaviour near L p 8 \ Lp 1 ) For large finite R 8 > 0, consider the set of all t C, such that the solution at complex time t is close to L p 8\L p 1, with y 10 (t) R 8, but not close to L p 9 Then this set is the complement of D 9 in an approximate disk D 8 with centre τ and radius d(τ)r 8 More precisely, t y 10 defines a covering from the annular domain D 8 \D 9 onto the complement in {z C z R 8 } of an approximate disk with centre at the origin and small radius d(τ) R 9, where y 10 (t) d(τ) 1/ (t τ) Proof Set L p 8 \ Lp 1 is visible in the chart (y 10,z 10 ), where it is given by the equation z 10 = 0 and parametrized by y 10 C, see Section A13 In that chart, the line L p 9 (without one point) is given by the equation y 10 = 0 and parametrized by z 10 C The line L p 4 (without one point) is given by the equation y 10 = θ 1 e t and also parametrized by z 10 C For z 10 0, bounded e t, and y 10 bounded and bounded away from θ 1 e t, we have: z (43a) θ y 10 θ 1 e t 1 e t y 10, +y 10 (43b) (43c) (43d) (43e) From (43d): y 10, z 10 J 10 = y 10 (θ 1 e t +y 10 )z10, J 10 J 10 (1+η)θ 1e t y 10 +θ 1 e t +y 10, EJ 10 1 θ 1e t y 10 log J 10(t 1 ) J 10 (t 0 ) θ 1(e t 1 e t 0 )+(t 1 t 0 )K, K = (1+η)θ 1e τ y 10 ( τ) +y 10 ( τ),

17 ASYMPTOTICS OF P V 17 where τ is on the integration path Therefore J 10 (t 1 )/J 10 (t 0 ) 1, if forall tonthesegment fromt 0 tot 1 we have t t 0 t 0 and θ 1 e t y 10 (t) 1 t 0, y 10(t) 1 t 0 We choose t 0 on the boundary of D 9 from Lemma 45 Then we have which implies that d(τ) d(t 0 ) J 11(τ) J 10 (t 0 ) 1 and y 11(t 0 ) = R 9, z 10 = 1 y 11 +(1+η)θ 1 e t 1 1 R 9 Since D 9 is an approximate disk with centre τ and small radius d(τ) R 9, and R 9 τ 1, we have that y 11 (t) R 9 1 hence: and t t 0 t 0 y 10 1 if t = τ +r(t 0 τ), r 1, = (r 1) 1 τ 1 if r 1 1 t 0 1 τ t 0 Then equation (43c) and J 10 d(τ) yield ( ) z10 1 y10 (θ 1 e t 1/ +y 10 ) d(τ) which in combination with (43a) leads to Hence and therefore dy 10 dt d(τ) 1/ ( ) y 1/ 10, d(τ) y 10 (t) y 10 (t 0 )+d(τ) 1/ (t t 0 ), y 10 (t) d(τ) 1/ (t t 0 ) if t t 0 y 10 (t 0 ) For large finite R 8 > 0, the equation y 10 = R 8 corresponds to t t 0 d(τ)r 8, which is still small compared to t 0 τ, therefore t τ t t 0 + t 0 τ τ This proves the statement of the lemma Lemma 47 (Behaviour near L p 4) For large finite R 4 > 0, consider the set of all t C, such that the solution at complex time t is close to L p 4, with z 91 (t) R 4, but not close to L p 8 Then this set is the complement of D 8 in an approximate disk D 4 with centre τ and radius R 4 /d(τ) More precisely, t z 91 defines a covering from the annular domain D 4 \D 8 onto the complement in {z C z R 4 } of an approximate disk with centre at the origin and small radius d(τ) /R 8, where z 91 (t) d(τ)/(t τ) Proof L p 4 without one point is visible in the chart (y 91,z 91 ), where it is given by the equation y 91 = 0 and parametrized by z 91 C, see Section A1 In that chart, the line L 1 is not visible, while L 8 is given by the equation z 91 = 0

18 18 NALINI JOSHI AND MILENA RADNOVIĆ For y 91 0 and bounded z 91 and e t, we have: (44a) (44b) y 91 θ 1e t z 91, z 91 θ 1e t y 91, (44c) J 10 = y 91z91, θ 1 e t y 91 (44d) From (44d): J 10 J 10 1 η θ 1 e t log J 10(t 1 ) J 10 (t 0 ) (1+η)(t 1 t 0 ) θ 1 (e t 1 e t 0 ) ThereforeJ 10 (t 1 )/J 10 (t 0 ) 1if foralltonthesegment fromt 0 tot 1 wehave t t 0 t 0 We choose t 0 on the boundary of D 8 from Lemma 46 Then we have d(τ) d(t 0 ) J 10(τ) J 10 (t 0 ) 1 and y 10(t 0 ) = R 8, which implies that θ 1 e t = y 91 y 10 R 8 1 Hence: y 91 1 if θ 1 e t R 8 Then equation (44d) and J 10 d(τ) yield then, using (44b) we get: It follows that y 1 91 z 91 θ 1 e t d(τ), d(z91 1 ) 1 dt d(τ) z 1 91 z 91(t 0 ) 1 + t t 0 d(τ), and therefore z 91 d(τ) if t t 0 z 91 (t 0 ) 1 t t 0 For large finite R 4 > 0, the equation z 91 = R 4 corresponds to t t 0 d(τ)/r 4, which is small compared to t 0 τ, and therefore t τ t t 0 + t 0 τ τ This proves the statement of the lemma Theorem 48 Let ǫ 1, ǫ, ǫ 3 be given such that ǫ 1 > 0, 0 < ǫ < θ 1, 0 < ǫ 3 < 1 Then there exists δ 1 > 0 such that if e t 0 < ǫ 1 and d(t 0 ) < δ 1, then: satisfies: ρ = inf{r < e t 0 such that d(t) < δ 1 whenever e t 0 e t r} (i) δ 1 d(t 0 ) (ρ 1 e t 0 ) θ 1 ǫ (1 ǫ 3 ); (ii) if e t 0 e t ρ then d(t) = d(t 0 )e θ 1(t t 0 )+ǫ (t) (1+ǫ 3 (t)); (iii) if e t ρ then d(t) δ 1 (1 ǫ 3 )

19 ASYMPTOTICS OF P V 19 Proof Suppose a solution of the system () is close to the infinity set at times t 0 and t 1 It follows from Lemmas that for every solution close to I, the set of complex time t such that the solution is not close to L p 1 L p is the union of approximate disks of radius d Hence if the solution is near I for all complex times t such that e t 0 e t e t 1, then there exists a path P from t 0 to t 1, such that the solution is close to L p 1 L p for all t P and P is C 1 -close to the path: s t s 1t 1 s 0, s [0,1] Then Lemmas 41 and 4 imply that: Therefore and log E(t) E(t 0 ) = θ 1(t t 0 ) 1 0 dt+o(1) E(t) = E(t 0 )e θ 1(t t 0 )+o(1) (1+o(1)), (45) d(t) = d(t 0 )e θ 1(t t 0 )+o(1) (1+o(1)) From Lemmas 44 47, we then have that, as long as the solution is close to I, the ratio of d remains close to 1 For the first statement of the theorem, we have: and so δ 1 > d(t) d(t 0 ) e θ 1(t t 0 ) ǫ (1 ǫ 3 ) δ 1 sup d(t 0 ) e θ 1 (t t 0 ) ǫ (1 ǫ3 ) {t d(t) <δ 1 } The second statement follows from (45), while the third follows by the assumption on t In the following corollary, we summarise the results obtained in this section Corollary 49 No solution of () intersects I A solution that approached I will stay in its vicinity for a limited range of the independent variable t Moreover, if a solution is sufficiently close to I at a point t, then it will have a pole in a neighbourhood of t Proof The first two statements follow from Theorem 48, and the last one from Lemmas The limit set In this section we consider properties of the limit set of thesolutions, when Ret, ie x 0 First, wedefinethelimit set, generalising theconcept oflimitsets indynamical systems Definition 51 Let (y(t),z(t)) be a solution of () The limit set Ω (y,z) of (y(t),z(t)) is the set of all S F \I such that there exists a sequence t n C satisfying: lim Ret n = and lim(y(t n ),z(t n )) = S n n Theorem 5 There exists a compact subset K of F \I, such that the limit set Ω (y,z) of any solution (y,z) is contained in K Moreover, Ω (y,z) is a non-empty, compact and connected set, which is invariant under the flow of the autonomous system (5)

20 0 NALINI JOSHI AND MILENA RADNOVIĆ Proof For any positive numbers δ 1, r, let K δ1,r denote the set of all s F(t) such that e t r andd(s) δ 1 Since F(t)isacomplex analyticfamilyover Cofcompact surfaces, K δ1,r is also compact Furthermore K δ1,r is disjoint from the union of the infinity sets I(t), t C, and therefore K δ1,r is a compact susbset of the Okamoto space O\F 0 When r approches zero, the sets K δ1,r shrink to the set K δ1,0 = {s F(0) d(s) δ 1 } F \I, which is compact It follows from Theorem 48 that there exists δ 1 0 such that for every solution (y,z) there exists r 0 > 0 with the following property: (y(t),z(t)) K δ,r0 for every t such that e t r 0 In the sequel, we take r r 0, when it follows that (y(t),z(t)) K δ1,r whenever e t r Let T r = {t C e t r} and let Ω (y,z),r denote the closure of (y,z)(t r ) in O Since T r is connected and (y,z) continuous, Ω (y,z),r is also connected Since (y,z)(t r ) is contained in the compact subset K δ1,r, its closure Ω (y,z),r is also contained there and therefore Ω (y,z),r is a non-empty compact and connected subset of O\F(0) The intersection of a decreasing sequence of non-empty compact and connected sets is non-empty, compact and connected: therefore, as Ω (y,z),r decrease to Ω (y,z) when r tends to zero, it follows that Ω (y,z) is a nonempty, compact and connected set of O Since Ω (y,z),r K δ1,r for all r r 0, and the sets K δ1,r shrink to the compact subset K δ1,0 of F \ I as r tends to zero, it follows that Ω (y,z) K δ1,0 This proves the first statement of the theorem with K = K δ1,0 Since Ω (y,z) is the intersection of the decreasing family of compact sets Ω (y,z),r, there exists for every neighbourhood A of Ω (y,z) in O and r > 0 such that Ω (y,z),r A, hence (y(t),z(t)) A for every t C such that e t r If t j is any sequence in C such that Ret j, then the compactness of K δ1,r, in combination with (y,z)t r K δ1,r, implies that there is a subsequence j = j(k) as k, such that: (y(t j(k) ),z(t j(k) )) s as k Then it follows that s Ω (y,z) Next we prove that Ω (y,z) is invariant under the flow Φ τ of the autonomous Hamiltonian system (5) Let s Ω (y,z) and t j be a sequence in C such that Ret j and (y(t j ),z(t j )) s Since the t-dependent vector field of the system () converges in C 1 to the vector field of the autonomous system (5) as Ret, it follows from the continuous dependence on initial data and parameters, that the distance between (y(t j + τ),z(t j + τ)) and Φ τ (y(t j ),z(t j )) converges to zero as j Since Φ τ (y(t j ),z(t j )) Φ τ (s) and Ret j as j, it follows that (y(t j + τ),z(t j + τ)) Φ τ (s) and t j +t as j, hence Φ τ (s) Ω (y,z) Proposition 53 If y is a solution of (11) with essential singularity at x = 0, than the flow (y,z) of the vectory field () meets each of the pole lines L 5, L 6, L 7 infinitely many times Proof First, suppose that a solution (y(t),z(t)) intersects the union L 5 L 6 L 7 only finitely many times According to Theorem 5, the limit set Ω (y,z) is a compact set in F \ I If Ω (y,z) intersects one of the pole lines L 5, L 6, L 7 at a point p, then there exists t with arbitrarily large negative real part such that (u(t), v(t)) is arbitrarily close to p, when the transversality of the vector field to the pole line implies thet (y(τ),z(τ)) L 5 L 6 L 7 for a

21 ASYMPTOTICS OF P V 1 unique τ near t As this would imply that (y(t),z(t)) intersects L 5 L 6 L 7 infinitely many times, it follows that Ω (y,z) is a compact subset of F \ (I L 5 L 6 L 7 ) It follows that Ω (y,z) is a compact subset contained in the first affine chart, which implies that y and z remain bounded for large negative Ret Thus y, z are holomorphic functions of x = e t in a neighbourhood of x = 0, which implies that there are complex numbers y( ), z( ) which are the limit points of y(t), z(t) as Ret That means that y is analytic at x = 0, which contradicts the assumption that it has there an essential singularity Since the limit set Ω (y,z) is invariant under the autonomous flow, it means that it will contain the whole irreducible component of a curve from the pencil h c (y,z) = 0 given by (31), for some constant c It is shown in Section 31 that this pencil of curves is birationally equivalent to a pencil of conics We identified in Section 31 the three singular conics in the pencil and found the special solutions corresponding to them In all other cases, all three base points b 5, b 6, b 7 will be contained in the limit set, which are projections of the pole lines L 5 ( ), L 6 ( ), L 7 ( ) respectively For a general solution (y,z), the base point b 4 will not be contained in the limit set, because that point is not a base point of the autonomous system (4) Remark 54 If the limit set Ω (y,z) contains only one point, that point must be a fixed point of the autonomous system (5) As we obtained in Section, there are four such points One of the points has y-coordinate equal to unity and it corresponds to the rational solutions of the form κ κ+x and x+κ κ x Theorem 55 Every solution of (11) with essential singularity at x = 0 has infinitely many poles and infinitely many zeroes in each neighbourhood of that singular point Proof Applying results from Section 4 and Proposition 53, we get that each solution has a simple pole at the intersections with L 5 and L 6 and a simple zero at the intersection with L 7 6 Limit x For studying the limit x, it is convenient to represent the P V as the following system: (61) y = x( 1 ) y(y 1) z (θ 0 +η)y +(θ 0 +η θ 1 x)y θ 0, z = 1 ( (y 1)(3y 1)z ( (θ 0 +η)y θ 0 η +θ 1 x ) z + 1 ) x ǫ(θ 0 +η θ ), where θ = α, θ 0 = β, θ 1 = δ (θ 1 0), η = γ θ 1 1, ǫ = 1 (θ 0 +θ +η) Remark 61 Using the change of the independent variable t = logx, equation (61) will give P V in the form (1) The resolution of the singularities of (61) will lead to the same space of the initial values as described in Section, and shown in Figure 3 In the limit x, the fifth Painlevé equation (11) becomes: (6) y = ( 1 y + 1 y 1 ) y + δy(y +1) y 1,

22 NALINI JOSHI AND MILENA RADNOVIĆ which has a first integral: y y(y 1) + δy (y 1) The solutions of (6) are elliptic functions satisfying: y = y ( C(y 1) δy ), where C is an arbitrary constant Notice that y 0 is the only solution of that equation taking the value 0 Thus, in contrast to the case when x 0, a 3 is not a base point for the autonomous equation, while a 4 will be a base point for that equation Now, analysing the system (61) in the similar way as shown in Sections 4 and 5, it follows in the limit x, that a general solution of that system has a compact limit set which is invariant with the respect to the autonomous flow This implies, as in the proof of Proposition 53, that the flow (y,z) of the vector field (61) meets each of the pole lines L 5, L 6, L 10 infinitely many times Thus every solution of (11) with essential singularity at x = has infinitely many poles and take the value 1 infinitely many times in each neighbourhood of that singular point Appendix A Resolution of the system In this section, we give the explicit construction of the space of initial conditions for () This constructions consists of eleven successive blow-ups of points in CP Weusethefollowingnotation ThecoordinatesinthreeaffinechartsofCP aredenoted by (y 01,z 01 ), (y 0,z 0 ), and (y 03,z 03 ) The exceptional line obtained in the n-th blow-up is covered by two coordinate charts, denoted by (y n1,z n1 ) and (y n,z n ) In each of these charts, we write the system () in the corresponding coordinates and look for base points the points contained by infinitely many solutions We calculate the coordinates of the base points in local coordinates in the following way In each chart (y nj,z nj ), the system can be written in the form: y nj = P(y nj,z nj,e t ) Q(y nj,z nj,e t ), z nj = R(y nj,z nj,e t ) S(y nj,z nj,e t ), for some polynomial expressions P, Q, R, S The uniqueness property for the given initial conditions is broken whenever P = Q = 0 or R = S = 0, so solving these equations yields to base points We note that, after blowing up, a new base points can appear only on the exceptional line We remark that a base point in algebraic geometry is a joint point of all curves from a given pencil The solutions of the autonomous system (5) are algebraic curves from the pencil (31), hence the notions of base points of a system of differential equations and base points of a pencil of curves coincide in the autonomous case A1 Affine chart (y 01,z 01 ) The first affine chart is defined by the original coordinates: y 01 = y, z 01 = z The energy (6) is: E = y(y 1) z (θ 0 +η)y z +(θ 0 +η)yz θ 0 z + 1 ǫ(θ 0 +η θ )y

23 ASYMPTOTICS OF P V 3 A Affine chart (y 0,z 0 ) The second affine chart is given by the coordinates: y 0 = 1 y, z 0 = z y, y = 1 y 0, z = z 0 y 0 The line at infinity is L : y 0 = 0 The Painlevé vector field () is: y 0 = z 0 y 0 + 4z 0 y 0 +η+θ 0 (η +θ 0 )y 0 +θ 0 y 0 z 0 +θ 1 e t y 0, z 0 = 5z 0 + 8z 0 +y y0 3 y0 0 (θ 0 z 0 F)+ 3z 0(η +θ 0 z 0 ) (η +θ 0 )z 0 +θ 1 e t z 0 y 0 In this chart, there is one visible base point: a 0 (y 0 = 0,z 0 = 0) A3 Affine chart (y 03,z 03 ) The coordinates: y 03 = y z, z 03 = 1 z, y = y 03, z = 1 z 03 z 03 The line at infinity is L : z 03 = 0 The vector field () is: y 03 = 5y3 03 z y 03 z y 03(1 (η+θ 0 )y 03 ) z 03 +(Fy 03 θ 0 )z 03 +(η +θ 0 )y 03 θ 1 e t y 03, z 03 = 3y 03 4y 03 +(η +θ z03 0 )z 03 +Fz03 z +1 (η +θ 0)y 03 θ 1 e t z In this chart, there is one base point: A4 Resolution at a 0 First chart: a 1 (y 03 = 0,z 03 = 0) y 11 = y 0 z 0 = 1 z, z 11 = z 0 = z y, y = 1 y 11 z 11, z = 1 y 11 The exceptional line is L 0 : z 11 = 0, while the proper preimage of the line at infinity L is given by y 11 = 0 The vector field () in this chart is: y 11 = 3 (η +θ 0) 4 +1+(η+θ y11 0 )y 11 +Fy z z 11 y 11 z θ 1e t y 11, 11 z 11 = (η +θ 0) 3z 11 (η +θ y11z 3 11 y11 0 )z 11 Fy 11 z 11 +θ 0 y 11 z11 y 11 y +θ 1e t z In this chart, there are no base points on the exceptional line L 0

24 4 NALINI JOSHI AND MILENA RADNOVIĆ Second chart: y 1 = y 0 = 1 y, z 1 = z 0 y 0 = z, y = 1 y 1, z = z 1 The exceptional line is L 0 : y 1 = 0, while the proper preimage of the line at infinity L is not visible in this chart The vector field () is: y 1 = z 1 y 1 +η +θ 0 (η +θ 0 )y 1 +θ 0 y 1 +4z 1 y 1 z 1 +θ 1 e t y 1, z 1 = 3z 1 + 4z 1 + (η +θ 0)z 1 F (η +θ y1 0 )z 1 z1 y 1 y +θ 1e t z 1 1 In this chart, there is one base point on the exceptional line L 0 : A5 Resolution at a 1 First chart: a (y 1 = 0,z 1 = 0) y 1 = y 03 z 03 = y, z 1 = z 03 = 1 z, y = y 1, z = 1 z 1 The exceptional line is L 1 : z 1 = 0, while the proper preimage of the line at infinity L is not visible in this chart The vector field () is: y 1 = (y 1 1) y 1 z 1 θ 0 +(η +θ 0 θ 1 e t )y 1 (η+θ 0 )y 1, z 1 = 1 4y 1 +3y 1 +(η +θ 0 θ 1 e t )z 1 (η +θ 0 )y 1 z 1 +Fz 1 In this chart, there are two base points on the exceptional line L 1 : a 3 (y 1 = 0,z 1 = 0), a 4 (y 1 = 1,z 1 = 0) The energy (6) and its rate of change are: Second chart: E = y 1(y 1 1) z 1 E = θ 1 e t ( y1 y 3 1 z 1 (y 1 1) ( ) (η+θ 0 )y 1 θ 0 +Fy 1, z 1 + (η +θ ) 0)y1 θ 0 Fy 1 z 1 y = y 03 = y z, z = z 03 y 03 = 1 y, y = 1 z, z = 1 y z The exceptional line is L 1 : y = 0, while the proper preimage of the line at infinity L is z = 0

25 ASYMPTOTICS OF P V 5 The vector field () is: y = (1 (η+θ 0)y ) θ z 3 z 0 y z +Fy z z +(η+θ 0 θ 1 e t )y, z = (z 1) +η +θ y z 0 (η +θ 0 θ 1 e t )z +θ 0 z In this chart, the only visible base point on the exceptional line L 1 is (y = 0,z = 1), which is a 4 The energy (6) and its rate of change are: E = 1 y z5 E = θ 1e t yz 5 y z4 + 1 y z3 η +θ 0 y z 3 + η +θ 0 y z θ 0 y z + F z, ( Fy z 4 θ 0y z 4 +(η +θ 0)y z +z 1) A6 Resolution at a First chart: y 31 = y 1 z 1 = 1 yz, z 31 = z 1 = z, y = 1 y 31 z 31, z = z 31 The exceptional line is L : z 31 = 0, while the proper preimage of the line L 0 is y 31 = 0 The vector field () is: y 31 = 1 (η +θ 0)y 31 +Fy 31 y 31 z 31 y 31 z 31 +θ 0 y 31 z 31, z 31 = 3 + (η +θ 0) + 4z 31 (η +θ y31 0 θ 1 e t )z 31 z31 F y 31 y 31 Therearetwobasepointsontheexceptional linel ; theiry 31 coordinatesarethesolutions of the quadratic equation: 1 (η +θ 0 )y 31 +Fy 31 = 0 Thus, the base points are ) a 5 (y 31 =,z 31 = 0, a 6 (y 31 = θ 0 +η+θ The energy (6) and its rate of change are: E = 1 y31 3 z η +θ 0 31 y31 z + F 31 y31 E = θ 1 e t ( 1 y 3 31z 31 η+θ 0 y 31z 31 + F y 31 z 31 z 31 y 31 +θ 0 z 31 ),z 31 = 0 θ 0 +η θ + θ 0 +η +z 31 θ 0 z 31, y 31 z 31 y 31 )