Integrability of P Es
|
|
- Judith Clark
- 5 years ago
- Views:
Transcription
1 Lecture 3 Integrability of P Es Motivated by the lattice structure emerging from the permutability/superposition properties of the Bäcklund transformations of the previous Lecture, we will now consider the integrability properties of these viewed as partial difference equations (P Es) on the two-dimensional space-time lattice. What we will discover is that the integrability can be given a precise, even algorithmic, meaning. The presence of parameters (namely the Bäcklund parameters λ and µ, which we will now reinterpret as lattice parameters) will play a crucial role in the development of the theory. Furthermore, as we will see in subsequent Lectures, the parameters render the P Es very rich: since they can be seen to represent the widths of the underlying lattice grid they allow us to recover, through continuum limits, a great wealth of other equations, semi-continuous (i.e. differential-difference type) as well as fully continuous (i.e. partial differential type) ones. The interplay between the discrete and continuous structures will prove to be one of the emerging features of the integrable systems that we study. The history of integrable difference equations goes back to seminal papers by Ablowitz & Ladik and by Hirota in the 1970 s, cf. [1,2]. The first was motivated by the search for integrable numerical algorithms through finite-difference approximations. This fits well into the general problem of the analysis of finite-difference P Es arising from numerical studies of PDEs, cf. e.g. the monograph by P.R. Garabedian, Ch. 13. More recently, systematic methods for the construction of integrable nonlinear finite-difference P Es were found, e.g. through the representation theory of infinite-dimensional Lie algebras, [3], or through singular linear integral equations and connections with Bäcklund transformations, cf. [4,5]. 3.1 Quadrilateral P Es We will investigate here partial difference equations (P E s) of the following canonical form (which we will call quadrilateral P Es) f(u, ũ, û, ũ) = 0, (3.1.1) 44
2 3.1. QUADRILATERAL P ES 45 where we adopt the canonical notation of vertices surrounding an elementary plaquette on a rectangular lattice: u := u n,m, ũ = u n+1,m û := u n,m+1, ũ = un+1,m+1 Schematically, this configuration of points is given by: u ũ û ũ The notation is inspired by the one for Bäcklund transformations, which as we have seen in the previous Lecture, give rise to purely algebraic equation as a consequence of their permutability property. However, here we will forget, at first instance, about this connection, and consider lattice equations of the form (3.1.1) in their own right as a partial difference equation (P E) on a two-dimensional lattice. Thus, n and m, play the role of independent (discrete) variables, very much like x and t being the continuous variables for an equation like the KdV equation Even though the form (3.1.1) seems very restrictive, it will turn out that it is in a sense the most elementary form as a model type of equations, and on the other hand remarkably rich. In fact, as we shall see later, P Es of the form (3.1.1) are rich enough to be candidates for discretisations of PDEs of arbitrary order in the spatial and temporal variables. A classification of P Es in the same way as of PDE s, does not yet exist. Nonetheless, we may consider a P E of the form (3.1.1) to have features reminiscent of hyperbolic PDEs. In fact, if the equation f = 0 can be solved uniquely for each elementary quadrilateral, then one may pose initial value problems (IVPs) in ways very similar to hyperbolic type of equations such as the KdV equation itself (i.e. as a discrete nonlinear evolution equation). This can be done as follows. The naive approach would be to consider IVP where we would assign values of the dependenet variable u along horizontal array of vertices in the lattice, i.e. values for u n,0, for all n, considering the variable m to be the temporal discrete variable. It is easy to see, however, that such an IVP would lead to a nonlocal problem if we want to use (3.1.1) as an iteration scheme to find all values u n,m for m > 0. In fact, to calculate any value u n,1 say for given n, we would need to involve all initial values u n,0 with n < n, and furthermore have to assume limiting behaviour as n. This would be a complicated procedure. However, there is nothing that tells us that we should identify the n- and m axes as the spatial and temporal axes respectively. The lattice picture allows us to play other and more
3 46 LECTURE 3. INTEGRABILITY OF P ES natural games, and it is the equation itself that gives the lead in provding is with the natural IVP that we could impose. Thus, changing the perspective slightly, we may tilt the lattice and rather consider initial value data to be imposed on configurations like a sawtooth (or a ladder ), such as: The analogy with hyperbolic PDEs can be further seen from the consideration of the memory each point in the lattice has of the initial values that are involved in its determination upon iteration. In fact, any given point has a backward shadow (the analogue of the so-called lightcone ) of points the values of u on which determine the value at that point, as indicated by the picture: the point at the bottom of the cone being fully determined by the initial values at the top of the cone and only and exclusively by those values! This is very much reminiscent again of what happens in the case of hyperbolic PDEs. Multilinearity In order to have a unique iteration scheme arising from an equation of the form (3.1.1) the equation should be linear in each of the variables around the quadrilateral, i.e. the equation should multilinear. The most general form of a quadrilateral lattice equation which is linear in each dependent variable around the quadrilateral, and which in addition respects reversal symmetry with respect to the shifts - and on the lattice, takes the following form: ) ( ) k 0 uûũ ũ + k 1 (uûũ + uû ũ + uũ ũ + ûũ ũ + k 2 ûũ + u ũ ) ( ) ( +k 3 (uũ + û ũ + k 4 uû + ũ ũ + k 5 u + ũ + û + ũ ) + k 6 = 0, (3.1.2) where k 1,..., k 6 are coefficients (which may depend on additional parameters). Exercise Show that, starting from a multilinear quadrilateral lattice equation (3.1.1) with 16 general coefficients, we arrive at the form (3.1.2) by assuming that the equation remains unchanged when we reverse the - or shifts, i.e. when we replace ũ by ũ (by which we mean the backward shift related to the -shift, see picture), or when we replace û by ṷ (meaning the backward shift related to the -shift, see picture).
4 3.2. INTEGRABILITY AND CONSISTENCY-AROUND-THE-CUBE 47 ṷ ũ u ũ u û It is easy to see that the lattice equation (2.2.58) derived in Lecture 2 from the Bäcklund transform of KdV is precisely of the form given above. However, not all equations of the form (3.1.2) will be of interest to us. We will be interested in those equations (for specific choices of the coefficients k 1,...,k 6 ) which have the additional property of being integrable. What we mean by this will be explored in the next section. 3.2 Integrability and Consistency-around-the-Cube We will now consider a class of quadrilateral P Es (3.1.1) in which, apart from the independent discrete variables n, m on which the variable u n,m depends, there are parameters which we associate with these independent variables. We can think of these parameters as being the parameters which measure the width of the grid in the directions associated with n and m, and we refer to them as lattice parameters. Thus, denoting these parameters by p, q, the equations under question will take the form: f(u, û, ũ ũ; p, q) = 0, (3.2.1) and if we demand that the values for u at each vertex can be solved uniquely, implying multilinearity in each of these variables around the quadrilateral, then we are led again to eq. (3.1.2) with the coefficients k 0,..., k 6 depending on the lattice parameters p and q in a specific way. The question is what criteria to use in order to chose that dependence! It is here that we will restrict ourselves to quadrilateral P Es which we regard to be integrable. The question of what is the proper definition integrability, and to answer that question in general is very difficult: a one-fits-all definition (for all the possible type of systems that we would like to regard as being integrable) is possibly not possible to give in a precise mathematical sense. However, if we restrict ourselves here to P Es, and in particular quadrilateral P Es, then we can aspire to be a bit more precise. We will explore the definition of an integrable quadrilateral P E by means of the following example, namely the example of the P E arising from the BTs for the KdV equation, eq. (2.2.58). First, we remark that the presence of the lattice parameters p, q is crucial: whereas normally we would like to consider the parameters to be chosen once and for all and then remain fixed (thus specifying a specifi equation) when solving the equation on the lattice, here we will argue that we should look at (2.2.58) as defining a whole parameter-family of equations, and that it makes sense to look at them altogether with p and q variable. In
5 48 LECTURE 3. INTEGRABILITY OF P ES doing this, however, we must attach each parameter to a specific diecrete variable such as p being associated with the variable n, and q with m. This is also natural from the way in which we derived (2.2.58) from the BT construction: each BT is attached to a parameter (λ, µ,... ), and with each parameter we can build a direction in an infinite-dimensional lattice of BTs. The main point we want to make now is the following: Statement: the infinite parameter-family of P Es represented by the lattice equation, such as (2.2.58) is compatible, i.e. in each quadrilateral sublattice of the infinite-dimensional lattice we can consistently impose a copy of the P E in terms of the relevant discrete variable and associated with a corresponding lattice parameter. Let us illustrate this by means of the example of (2.2.58), in which we will identify (for later convenience) the parameters 4λ := p 2 and 4µ := q 2. What the statement above suggests is that we can embed the equation in a multi-dimensional lattice by conbsidering the dependent variable w, not only to depend on n and m (with associated lattice paameters p and q espectively), but that w may in fact depend on an infinity of lattice variables each associated with its parameter, as follows: w = w n,m,h,... = w(n, m, h,... ; p, q, r,...) and with each of these variables we have a corresponding elementary shift on the lattice: w := w n+1,m,h, ŵ := w n,m+1,h, w := w n,m,h+1... Rewriting the equation (2.2.58) in the form (2.2.50) with the substitutions for λ and µ investigate what happens if we impose a copy of the same equation in all three lattice directions. This would lead to the system of equations: (ŵ w)(w w) = p 2 q 2, (3.2.2a) (w w)(w w) = p 2 r 2, (3.2.2b) (ŵ w)(w ŵ) = r 2 q 2. (3.2.2c) These equations are consistent if the evaluations along the cube are independent of the way of calculating the final point. Imposing initial values: w := a, w =: b, ŵ =: c, w =: d
6 3.2. INTEGRABILITY AND CONSISTENCY-AROUND-THE-CUBE 49 w w w w ŵ w ŵ w In fact: w = (p2 q 2 ) wŵ + (q 2 r 2 )ŵw + (r 2 p 2 )w w (r 2 q 2 ) w + (q 2 p 2 )w + (p 2 r 2 )ŵ. (3.2.3) independent of the way in which he value at this vertex is calculated! This property, implying that under relevant initial value problems on the 3-dimensional lattice the iteration of the solution can be performed in an unambiguous way (namely independent of the way by which we perform the calculation of the iterates) will be referred to as the consistency-around-thecube (CAC) property. It is this property that we shall consider to be the main hallmarfk of the integrability of the equation. Remark: We note also, for later reference, that in the above case the formula for w is independent of the value w at the opposite end of the main diagonal across the cube. This property we will refer to as the tetrahedron property. Other examples of integrable quadrilateral P Es: There are many examples of integrable quadrilateral P Es that have been discovered over the years. A number of these belong to the lattice KdV-family of equations, they comprise the following cases: Lattice potential KdV equation: (p q + û ũ)(p + q + u ũ) = p 2 q 2, (3.2.4) and this is, in fact, equivalent to the equation (2.2.58) by the change of (dependent) variable w = u np mq c (c a constant with respect to n and m). We will show in Lecture 4 that this equation reduces to the potential KdV equation after a double continuum limit, so that we rightly regard it as a discretisation of the latter equation. Lattice potential MKdV equation: p(v v ṽ ṽ) = q(vṽ v ṽ) (3.2.5) Solutions of this equation are related to the solutions of the previous one (3.2.4) via the relations p q + û ũ = pṽ q v v, p + q + u ũ = pv + q ṽ ṽ, (3.2.6)
7 50 LECTURE 3. INTEGRABILITY OF P ES which constitute the analogues of the Miura transformation (2.2.42). Lattice SKdV equation: (z ẑ)( z z) p2 = (z z)(ẑ z) q 2 (3.2.7) which is the Schwarzian KdV equation. This equation is invariant under Möbius transformations: z Z = αz + β γz + δ. Solutions of eq. (3.2.7) are related to the ones of (3.2.5) via the relations: p(z z) = vṽ, q(z ẑ) = v v. (3.2.8) Exercise Show by explicit computation that the P Es given by (3.2.4)-(3.2.7) possess the consistency-around-the-cube property. Exercise Show that by using the relations (3.2.6) one can derive (3.2.4) from (3.2.5) and vice versa. 3.3 Lax pair for lattice potential KdV We shall next demonstrate how the CAC property explained in the previous section gives rise to the existence of a Lax pair, i.e. an overdetenrmined linear system of difference equations, the compatibility of which is verified iff the nonlinear lattice equation is satisfied. The idea is the following: Having verified the consistency of the equation around the cube, this tells us that we can add any lattice direction to the original lattice and impose simultaneously the equation in the three two-dimensional quadrilateral sublattices. The main idea is now to consider the additional lattice variable h Z associated with lattice parameter k as an auxiliary virtual variable, whilst acknowledging only the shifts in the orignal lattice variables n and m as the operations of interest. This then suggests that the shift in the third direction: w w should not appear in any of the equations, implying that wherever w appears we should treat it as a new dependent variable w := W. Proceeding in this way, we rewrite (3.2.2b) and (3.2.2c) as follows: (W w)( W w) = k 2 p 2 W = ww + (k2 p 2 ww) W w (W ŵ)(ŵ w) = k2 q 2 W = ww + (k2 p 2 ww) W w, (3.3.1a). (3.3.1b) Noting that eqs. (3.3.1) are both fractional linear in W, we can linearise these equations by the substitution: W = F G,
8 3.3. LAX PAIR FOR LATTICE POTENTIAL KDV 51 leading to F G = wf + (k2 p 2 w w)g, F wg F Ĝ = wf + (k2 q 2 wŵ)g, F ŵg and since at least one of the two functions F or G can be chosen freely, this allows us to split in each of these equations the numerator and denominator to give: { ( F = γ wf + [k 2 p 2 w w]g ) {, ( F = γ wf + [k respectively 2 q 2 wŵ]g ), G = γ (F wg), Ĝ = γ (F ŵg), (3.3.2) in which γ and γ are to be specified later. What happens next is obvious: we introduce the two-component vector ( ) F φ =, G and write (3.3.2) as a system of two 2 2 matrix equations: φ = Lφ, φ = Mφ, (3.3.3a) with the matrices ( w k L = γ 2 p 2 w w 1 w ), M = γ ( w k 2 q 2 wŵ 1 ŵ ) (3.3.3b) How does this linear system work? The consistency relation of the linear problem (3.3.3a) is obtained from the condition that φ = φ, which is the condition that simply expresses that φ must be a proper function of the lattice variables n and m. Calculating the left hand-side of this condition we get on the one hand φ = (Lφ) = L φ = LMφ, whereas the right hand-side would be calculated as φ = (Mφ) = M φ = MLφ. Equating both sides we see that a sufficient condition for the consistency is the matricial equation: LM = ML, (3.3.4) which we will loosely refer to as the Lax equation, but it is sometimes also referred to as a discrete zero-curvature condition (the reason for this terminology will be explained in (3.3.1)). Pictorially the Lax equation is illustrated by the following diagram
9 52 LECTURE 3. INTEGRABILITY OF P ES φ L φ M M φ L φ where the vectors φ are located at the vertices of the quadrilateral and in which the matrices L and M are attached to the edges linking the vertices. Using the explicit form (3.3.3b) of the matrices L and M, and working out the condition (3.3.4) we find ( ) ( ) γγ ŵ k 2 p 2 ŵ w w k 2 q 2 wŵ = 1 w 1 ŵ ( ) ( ) = γ w k γ 2 q 2 w w w k 2 p 2 w w 1 w 1 w and we will chose the γ and γ (which were unspecified so far) such that the relation for the determinants of this equation: γγ det( L) det(m) = γ γ det( M)det(L) (3.3.5) is trivially satisified. Since in this example the determinants of L and M are given by det(l) = p 2 k 2, det(m) = q 2 k 2, respectively, it is clear that the condition (3.3.5) is satisified by simply taking γ = γ = 1. But we will encounter other examples where a nontrivial choice of γ and γ is needed in order to satisfy the determinantal condition (3.3.5). Working out all the entries on both sides of the above matrix products it is straightforward to see that the (1,1)- and (2,2) entry of the matrix equation both yield the same condition on w, namely the equation (2.2.50)! Moreover, the (2,1) entry is trivially, and the (1,2) entry we don t even need to calculate, because having checked three of the entries to be satisfied the final entry must also be satisfied by virtue of the fact that we have aranged the determinantal condition (3.3.5) to be satisfied. In conclusion, we see thus that from the Lax equation (3.3.4) we recover the lattice equation (2.2.58). We observe furthermore, that albeit both Lax matrices L and M depend on the auxiliary variable k the final nonlinear equation for w does not depend on k! Exercise Suppose that the Lax matrices L and M can be expanded in a power series in a small parameter δ and ǫ respectively as follows: L = 1 + δl 1 +, M = 1 + ǫm 1 +,
10 3.4. CLASSIFICATION OF QUADRILATERAL P ES 53 and we expand the shifted variable by Taylor expansion L = L + ǫ t L +, M = L + δ x M +, by expansion we obtain in dominant order (namely terms proportional to ǫδ) the following matricial equation: t L 1 x M 1 + [ L 1, M 1 ] = 0, (3.3.6) where [ L, M ] = LM ML denotes the usual matrix commutator bracket. Eq. (3.3.6) arises also in differential geometry and it is from there that it has the interpretation as a zero-curvature condition of differential manifolds (with the relevant interpretation of the matrices L and M). 3.4 Classification of Quadrilateral P Es In a beautiful paper by V. Adler, A. Bobenko and Yu. Suris, cf. ref. [7], the problem of classifying all quadrilateral lattices which are integrable in the sense of the CAC property discussed in section 3.2 was considered. Since this is in our view an important result, it is useful to reproduce the whole list of resulting equations here. Theorem ([7]): Consider quadrilateral P Es of the general form: Q(u, ũ, û, ũ; p, q) = 0, using the notation indicated in the following diagram (renaming the lattice parameters in order to avoid confusion with previously used notation), u p ũ q q û p ũ subject to the following restrictions: a) Linearity: Q is multilinear in its arguments, i.e., it is linear in each vertex-variable u,ũ, û, ũ; b) Symmetry: Q is invariant under the group D 4 of symmetries of the square, generated by the interchangements: Q(u, ũ, û, ũ; p, q) = ±Q(u, û, ũ, ũ; q, p) = ±Q(ũ, u, ũ, û; p, q) ; (3.4.1) c) Tetrahedral Condition: in the consistency check, the evaluation of the point on the cube given by ũ is independent of u.
11 54 LECTURE 3. INTEGRABILITY OF P ES Q-list: Q 1 : p(u û)(ũ ũ) q(u ũ)(û ũ) = δ 2 pq(q p) (3.4.2a) Q 2 : p(u û)(ũ ũ) q(u ũ)(û ũ) + pq(p q)(u + ũ + û + ũ) = = pq(p q)(p 2 pq + q 2 ) (3.4.2b) Q 3 : p(1 q 2 )(uû + ũ ũ) q(1 p 2 )(uũ + û ũ) = ( = (p 2 q 2 ) (ûũ + u ũ) + δ 2 (1 p2 )(1 q 2 ) ) 4pq Q 4 : p(uũ + û ũ) q(uû + ũ ũ) = pq qp ( ) = 1 p 2 q 2 (ûũ + u ũ) pq(1 + uũû ũ) where P 2 = p 4 γp 2 + 1, Q 2 = q 4 γq H-list: (3.4.2c) (3.4.2d) H 1 : (u ũ)(û ũ) = p 2 q 2 (3.4.3a) H 2 : (u ũ)(ũ û) = (p q)(u + ũ + û + ũ) + p 2 q 2 (3.4.3b) H 3 : p(uũ + û ũ) q(uû + ũ ũ) = δ 2 (p 2 q 2 ) (3.4.3c) A-list: A 1 : p (u + û)(ũ + ũ) q(u + ũ)(û + ũ) = δ 2 pq(p 2 q 2 ) (3.4.4a) A 2 : p(1 q 2 )(uû + ũ ũ) q(1 p 2 )(uũ + û ũ) + (p 2 q 2 )(1 + uũû ũ) = 0 (3.4.4b) The Q-equations are related through the following coalescence diagram, i.e. displaying how the equations relate to eachother through certain limits on the parameters: Q 3 (Q 3 ) δ=0 Q 4 Q 1 Q 2 (Q 1 ) δ=0 whilst the H-equations and A-equations are related through Miura type Bäcklund transforms, (J. Atkinson, 2008).
12 3.5. LATTICE KDV EQUATION 55 Remarks: It is in itself remarkable that the list of equations found is so short. In fact, all equations in the list Q are in fact special subcases of the last equation, Q4, which was discovered by V. Adler in The latter equation can also be expressed as: [ ] A[(u b)(û b) (a b)(c b)] +B [(u a)(ũ a) (b a)(c a)] (ũ b)( ũ b) (a b)(c b) [ (û a)( ũ a) (b a)(c a) + ] = = ABC(a b) (3.4.5) cf. [9]. Here the extra parameters c, C are related through C 2 = 4c 3 g 2 c g 3 and to (a, A) and (b, B) through the relations: A(c b) + B(c a) = C(a b), a + b + c = 1 ( ) 2 A + B. 4 a b A third alternative form of Adler s equation reads as follows: ( ) ( ) ( ) sn(α) vṽ + v ṽ sn(β) v v + ṽ ṽ sn(α β) ṽ v + v ṽ ( ) +ksn(α)sn(β)sn(α β) 1 + vṽ v ṽ = 0, (3.4.6) (J. Hietarinta, 2005) which is somewhat simpler. In (3.4.6) the sn denote the Jacobi elliptic functions, cf. Appendix A, of modulus k, i.e. sn(α) = sn(α; k), etc. It should be noted that all previous examples presented in section 3.2 can be recognised as special subcases of the equations in the lists Q,H and A, and which can be obtained by degeneration of Adler s equation in either of its forms. All these equations possess Lax pairs, but it was only recently that explicit solutions were found for most of the new cases including of Adler s equation, cf. e.g. [12]. Adler s equation is the discrete analogue of a famous soliton equation discovered in 1980 by I. Krichever and S. Novikov, which reads: u t = u xxx 3 u 2 xx (4u3 g 2 u g 3 ), (3.4.7) 2 u x and which generalises the Schwarzian KdV equation. 3.5 Lattice KdV Equation We finish this Lecture by deriving the lattice analogue of the KdV equation, in contrast to the lattice potential KdV equation (2.2.58) (or equivalently (3.2.4)): (w n,m w n+1,m+1 )(w n,m+1 w n+1,m ) = p 2 q 2. (3.5.1) We recall that the difference between the continuous KdV equation (2.2.38) and its potential version (2.2.46), as discussed in Lecture 2, was simply that its solutions were connected
13 56 LECTURE 3. INTEGRABILITY OF P ES through taking a derivative, u = w x. However, on the lattice there are many ways in which can do the analogue of taking the derivative, for instance we can replace it by taking a difference involving two neighbouring vertices, like n w n,m w n+1,m w n,m or m w n,m w n,m+1 w n,m, but these are not the only choices. Alternatively we can take a difference between vertices farther away or across diagonals, such as: w n,m = w n+1,m+1 w n,m or w n,m = w n,m+1 w n+1,m, or a host of other choices. As long as a well-chosen continuum limit reduces these to a derivative they can be justifiably be considered to be the discrete analogues of the operation of taking a partial derivative. The most sensible way to make a choice is to look at the equation at hand and then decide what would be for the given equation the most natural choice to apply in that case. In the case of the lattive potential KdV equation (3.5.1) the natural choice of the discrete analogue of the derivative seems to be either one of the choices given above, namely a difference across the diagonal. This suggest that as lattice KdV variables we would take either one of the two choices: Q = w w or R = ŵ w. (3.5.2) It is then straightforward from (3.2.4) (using the notation in terms of - and shifts) to derive the following equation for Q: Q Q = a Q a Q R R = a R ã R, (3.5.3) where a p 2 q 2. The equations for Q and R are simply related by the fact that QR = a. A Lax pair for either of the equations (3.5.3) can be easily derived, by starting from the Lax pair for the potential KdV equation. In fact, rather than writing it as a first order matricial system, one may derive from the Lax pair a three-point linear equation in terms of one of the components of the vector φ on which the Lax matrices act. In this way we obtain the following scalar Lax pair: remembering that QR = a. ϕ = Q ϕ + Λ ϕ, (3.5.4a) ϕ = ϕ + R ϕ, (3.5.4b) Exercise Derive the Lax pair (3.5.4) from the set of equations (3.3.2) (with γ = γ = 1) by setting ϕ = G and identifying the spectral variable by Λ = k 2 q 2. Exercise Show that the consistency condition of the Lax pair (3.5.4) viewed as an overdetermined discrete linear system, leads to the lattice KdV equation (3.5.3). 3.6 Singularity confinement In the seminal paper of B. Grammaticos et al. another point of view on integrability was developed than the one discussed in section 3.2, namely based on the issue of well-posedness
14 3.6. SINGULARITY CONFINEMENT 57 of initial value problems on the lattice. It is based on a phenomenon happening in integrable discrete systems, which nowadays is called singularity confinement and which describes what happens with singularities of the solutions of discrete system in the space of initial values. Proposition: In an integrable lattice equation of KdV type singularities induced by initial data do not propagate. As an example let us consider the lattice KdV equation (3.5.3), (setting for convenience the parameter a = 1 w.l.o.g.) R n+1,m+1 R n,m = 1 1. R n+1,m R n,m+1 We will now study how the the initial data progresses, when one hits a singularity. The initial data is given at the solid line on Figure 3.1, a, b, 0, c, d. When one proceeds from these initial values one obtains infinity at two places, then one 0 at the next level and finally two ambiguities of the type. A more detailed analysis with the initial value 0 1 = ε (small) yields the following values at the subsequent iterations at the first stage, and on the next 1 = b + 1 ε 1 a, 2 = c + 1 d 1 ε, s = a f, t = d + 1 g 1 2, 0 2 = ε ( = ε + b c a 1 ) ε d b c a 0 1 d f 1 2 g s 0 2? 1? 2 t Figure 3.1: Propagation of singularities in a 2D map. Here n grows in the SE direction and m in the SW direction.
15 58 LECTURE 3. INTEGRABILITY OF P ES Then at the next step we can resolve the ambiguities:? 1 = s = c + 1 d 1 a 1/f + O(ε)? 2 = t = b 1 a + 1 d + 1/g + O(ε) Thus the singularity is confined. Note that if the lattice equation is deformed, e.g. by taking: R n+1,m+1 R n,m = 1 R n+1,m λ R n,m+1 with λ 1, the above fine cancellation would no longer happen, and singularities would again occur at? 1 and? 2, and would persist throughout! In that case the singularities are no longer confined to a finite number of iteration steps, and we conclude that the corresponding map is not integrable. Remark: Another point of view on 2D singularity confinement is provided by the requirement of ultra-local singularity confinement (R. Sahadevan and H. Capel, 2003): One only assumes a singularity at one point and requires regularity at all other points. Referring to Figure 3.2 we have initial values are at black disks, the values at open circles are determined from them. The initial values can be chosen so that u 11 =, and the requirement is that u 12, u 21 are finite and ambiguity at u 22 can be resolved using ǫ-analysis. For example for the lattice potential KdV in the form the initial value u 01 = u 10 implies u n+1,m+1 = u n,m 1 u n,m+1 u n+1,m, u 11 =,, u 12 = u 01, u 21 = u 10,, u 22 =? u 02 u 12 u 22 u 01 u 11 u 21 u 00 u 10 u 20 Figure 3.2: Setting for ultra-local singularity confinement.
16 3.6. SINGULARITY CONFINEMENT 59 and a detailed analysis reveals that if u 10 = u 01 + ǫ then u 11 = 1 ǫ + u 00, u 12 = u 01 + ǫ + ǫ 2 (u 02 u 00 ) + O(ǫ 3 ), u 21 = u ǫ 2 ( u 20 + u 00 ) + O(ǫ 3 ), u 22 = 1 u 11 u 12 u 21 = 1 ǫ + u 1 00 ǫ + ǫ 2 (u 02 + u 02 2u 00 ) + O(ǫ 3 ) = u 02 + u 20 u 00 + O(ǫ) This resolves the singularity, and recovers initial value u 00 which was temporarily submerged to order ǫ 2. Literature 1. M.J. Ablowitz and F.J. Ladik, A nonlinear difference scheme and inverse scattering, Stud. Appl. Math. 55 (1976) ; On the solution of a class of nonlinear partial difference equations, ibid. 57 (1977) R. Hirota, Nonlinear partial difference equations I-III, J. Phys. Soc. Japan 43 (1977), , E. Date, M. Jimbo and T. Miwa, Method for generating discrete soliton equations I-V, J. Phys. Soc. Japan 51 (1982) , 52 (1983) , B. Grammaticos, A. Ramani and V. Papageorgiou, Do integrable mappings have the Painlevé property?, Phys. Rev. Lett. bf 67 (1991) F.W. Nijhoff, G.R.W. Quispel and H.W. Capel, Direct linearization of nonlinear differencedifference equations, Phys. Lett. 97A (1983) ; G.R.W. Quispel, F.W. Nijhoff, H.W. Capel and J. van der Linden, Linear integral equations and nonliner difference-difference equations, Physica 125A (1984) F.W. Nijhoff and H.W. Capel, The discrete Korteweg-de Vries equation, Acta Applicandae Mathematicae 39 (1995) F.W. Nijhoff and A.J. Walker, The discrete and continuous Painlevé VI hierarchy and the Garnier systems, Glasgow Math. J. 43A (2001) V.E. Adler, A.I. Bobenko and Yu.B. Suris, Classification of integrable equations on quadgraphs, Commun. Math. Phys. 233 (2003) F.W. Nijhoff, Lax pair for the Adler (lattice Krichever-Novikov) system, Phys. Lett. 297A (2002), V.E. Adler and Yu.B. Suris, Q4, Intl. Math. Res. Notes (2004). 11. J. Hietarinta, Searching for CAC-maps, J. Nonlin. Math. Phys. 12 Suppl. 2 (2005) J. Atkinson, J. Hietarinta and F.W. Nijhoff, Seed and soliton solutions for Adler s lattice equation, J. Phys A: Math. Theor. (40 (2007) F1-8.
arxiv:nlin/ v1 [nlin.si] 17 Sep 2006
Seed and soliton solutions for Adler s lattice equation arxiv:nlin/0609044v1 [nlin.si] 17 Sep 2006 1. Introduction James Atkinson 1, Jarmo Hietarinta 2 and Frank Nijhoff 1 1 Department of Applied Mathematics,
More informationarxiv: v1 [nlin.si] 25 Mar 2009
Linear quadrilateral lattice equations and multidimensional consistency arxiv:0903.4428v1 [nlin.si] 25 Mar 2009 1. Introduction James Atkinson Department of Mathematics and Statistics, La Trobe University,
More informationON THE SYMMETRIES OF INTEGRABLE PARTIAL DIFFERENCE EQUATIONS
Proceedings of the International Conference on Difference Equations, Special Functions and Orthogonal Polynomials, World Scientific (2007 ON THE SYMMETRIES OF INTEGRABLE PARTIAL DIFFERENCE EQUATIONS ANASTASIOS
More informationElliptic integrable lattice systems
Elliptic integrable lattice systems Workshop on Elliptic integrable systems, isomonodromy problems, and hypergeometric functions Hausdorff Center for Mathematics, Bonn, July 21-25, 2008 Frank Nijhoff University
More informationSoliton solutions to the ABS list
to the ABS list Department of Physics, University of Turku, FIN-20014 Turku, Finland in collaboration with James Atkinson, Frank Nijhoff and Da-jun Zhang DIS-INI, February 2009 The setting CAC The setting
More informationSymmetries and Group Invariant Reductions of Integrable Partial Difference Equations
Proceedings of 0th International Conference in MOdern GRoup ANalysis 2005, 222 230 Symmetries and Group Invariant Reductions of Integrable Partial Difference Equations A. TONGAS, D. TSOUBELIS and V. PAPAGEORGIOU
More informationSymmetry Reductions of Integrable Lattice Equations
Isaac Newton Institute for Mathematical Sciences Discrete Integrable Systems Symmetry Reductions of Integrable Lattice Equations Pavlos Xenitidis University of Patras Greece March 11, 2009 Pavlos Xenitidis
More informationgroup: A new connection
Schwarzian integrable systems and the Möbius group: A new connection March 4, 2009 History The Schwarzian KdV equation, SKdV (u) := u t u x [ uxxx u x 3 2 uxx 2 ] ux 2 = 0. History The Schwarzian KdV equation,
More informationGEOMETRY OF DISCRETE INTEGRABILITY. THE CONSISTENCY APPROACH
GEOMETRY OF DISCRETE INTEGRABILITY. THE CONSISTENCY APPROACH Alexander I. Bobenko Institut für Mathematik, Fakultät 2, Technische Universität Berlin, Strasse des 17. Juni 136, 10623 Berlin, Germany 1 ORIGIN
More informationFrom discrete differential geometry to the classification of discrete integrable systems
From discrete differential geometry to the classification of discrete integrable systems Vsevolod Adler,, Yuri Suris Technische Universität Berlin Quantum Integrable Discrete Systems, Newton Institute,
More informationarxiv: v1 [nlin.si] 29 Jun 2009
Lagrangian multiform structure for the lattice KP system arxiv:0906.5282v1 [nlin.si] 29 Jun 2009 1. Introduction S.B. Lobb 1, F.W. Nijhoff 1 and G.R.W. Quispel 2 1 Department of Applied Mathematics, University
More informationarxiv: v1 [nlin.si] 2 May 2017
On decomposition of the ABS lattice equations and related Bäcklund transformations arxiv:705.00843v [nlin.si] 2 May 207 Danda Zhang, Da-jun Zhang Department of Mathematics, Shanghai University, Shanghai
More informationCanonicity of Bäcklund transformation: r-matrix approach. I. arxiv:solv-int/ v1 25 Mar 1999
LPENSL-Th 05/99 solv-int/990306 Canonicity of Bäcklund transformation: r-matrix approach. I. arxiv:solv-int/990306v 25 Mar 999 E K Sklyanin Laboratoire de Physique 2, Groupe de Physique Théorique, ENS
More informationBäcklund transformation and special solutions for Drinfeld Sokolov Satsuma Hirota system of coupled equations
arxiv:nlin/0102001v1 [nlin.si] 1 Feb 2001 Bäcklund transformation and special solutions for Drinfeld Sokolov Satsuma Hirota system of coupled equations Ayşe Karasu (Kalkanli) and S Yu Sakovich Department
More informationSymbolic Computation of Lax Pairs of Nonlinear Systems of Partial Difference Equations using Multidimensional Consistency
Symbolic Computation of Lax Pairs of Nonlinear Systems of Partial Difference Equations using Multidimensional Consistency Willy Hereman Department of Applied Mathematics and Statistics Colorado School
More informationCoupled KdV Equations of Hirota-Satsuma Type
Journal of Nonlinear Mathematical Physics 1999, V.6, N 3, 255 262. Letter Coupled KdV Equations of Hirota-Satsuma Type S.Yu. SAKOVICH Institute of Physics, National Academy of Sciences, P.O. 72, Minsk,
More informationIntroduction to the Hirota bilinear method
Introduction to the Hirota bilinear method arxiv:solv-int/9708006v1 14 Aug 1997 J. Hietarinta Department of Physics, University of Turku FIN-20014 Turku, Finland e-mail: hietarin@utu.fi Abstract We give
More informationWeak Lax pairs for lattice equations
Weak Lax pairs for lattice equations Jarmo Hietarinta 1,2 and Claude Viallet 1 1 LPTHE / CNRS / UPMC, 4 place Jussieu 75252 Paris CEDEX 05, France 2 Department of Physics and Astronomy, University of Turku,
More informationarxiv: v2 [math-ph] 24 Feb 2016
ON THE CLASSIFICATION OF MULTIDIMENSIONALLY CONSISTENT 3D MAPS MATTEO PETRERA AND YURI B. SURIS Institut für Mathemat MA 7-2 Technische Universität Berlin Str. des 17. Juni 136 10623 Berlin Germany arxiv:1509.03129v2
More informationSingularities, Algebraic entropy and Integrability of discrete Systems
Singularities, Algebraic entropy and Integrability of discrete Systems K.M. Tamizhmani Pondicherry University, India. Indo-French Program for Mathematics The Institute of Mathematical Sciences, Chennai-2016
More informationOn universality of critical behaviour in Hamiltonian PDEs
Riemann - Hilbert Problems, Integrability and Asymptotics Trieste, September 23, 2005 On universality of critical behaviour in Hamiltonian PDEs Boris DUBROVIN SISSA (Trieste) 1 Main subject: Hamiltonian
More informationTHE LAX PAIR FOR THE MKDV HIERARCHY. Peter A. Clarkson, Nalini Joshi & Marta Mazzocco
Séminaires & Congrès 14, 006, p. 53 64 THE LAX PAIR FOR THE MKDV HIERARCHY by Peter A. Clarkson, Nalini Joshi & Marta Mazzocco Abstract. In this paper we give an algorithmic method of deriving the Lax
More informationA time-discretized version of the Calogero-Moser model
A time-discretized version of the Calogero-Moser model arxiv:hep-th/9403052v 8 Mar 994 Frank W. Nijhoff and Gen-Di Pang FB7 Mathematik-Informatik, Universität-Gesamthochschule Paderborn, D-33098 Paderborn,
More informationFrom Yang-Baxter maps to integrable quad maps and recurrences
From Yang-Baxter maps to integrable quad maps and recurrences arxiv:1206.1217v1 [nlin.si] 6 Jun 2012 B. Grammaticos IMNC, Universités Paris VII & XI, CNRS UMR 8165 Bât. 440, 91406 Orsay, France A. Ramani
More informationComputers and Mathematics with Applications
Computers and Mathematics with Applications 60 (00) 3088 3097 Contents lists available at ScienceDirect Computers and Mathematics with Applications journal homepage: www.elsevier.com/locate/camwa Symmetry
More informationOn Miura Transformations and Volterra-Type Equations Associated with the Adler Bobenko Suris Equations
Symmetry, Integrability and Geometry: Methods and Applications On Miura Transformations and Volterra-Type Equations Associated with the Adler Bobenko Suris Equations SIGMA 4 (2008), 077, 14 pages Decio
More informationA MULTI-COMPONENT LAX INTEGRABLE HIERARCHY WITH HAMILTONIAN STRUCTURE
Pacific Journal of Applied Mathematics Volume 1, Number 2, pp. 69 75 ISSN PJAM c 2008 Nova Science Publishers, Inc. A MULTI-COMPONENT LAX INTEGRABLE HIERARCHY WITH HAMILTONIAN STRUCTURE Wen-Xiu Ma Department
More informationFreedom in the Expansion and Obstacles to Integrability in Multiple-Soliton Solutions of the Perturbed KdV Equation
Freedom in the Expansion and Obstacles to Integrability in Multiple-Soliton Solutions of the Perturbed KdV Equation Alex Veksler 1 and Yair Zarmi 1, Ben-Gurion University of the Negev, Israel 1 Department
More informationSYMBOLIC COMPUTATION OF LAX PAIRS OF NONLINEAR PARTIAL DIFFERENCE EQUATIONS
SYMBOLIC COMPUTATION OF LAX PAIRS OF NONLINEAR PARTIAL DIFFERENCE EQUATIONS by Terry J. Bridgman c Copyright by Terry J. Bridgman, 2018 All Rights Reserved A thesis submitted to the Faculty and the Board
More informationSymmetry Reductions of (2+1) dimensional Equal Width. Wave Equation
Authors: Symmetry Reductions of (2+1) dimensional Equal Width 1. Dr. S. Padmasekaran Wave Equation Asst. Professor, Department of Mathematics Periyar University, Salem 2. M.G. RANI Periyar University,
More informationarxiv:solv-int/ v1 20 Oct 1993
Casorati Determinant Solutions for the Discrete Painlevé-II Equation Kenji Kajiwara, Yasuhiro Ohta, Junkichi Satsuma, Basil Grammaticos and Alfred Ramani arxiv:solv-int/9310002v1 20 Oct 1993 Department
More informationarxiv: v1 [nlin.si] 2 Dec 2011
DIRECT LINEARIZATION OF EXTENDED LATTICE BSQ SYSTEMS DA-JUN ZHANG & SONG-LIN ZHAO AND FRANK W NIJHOFF arxiv:1112.0525v1 [nlin.si] 2 Dec 2011 Abstract. The direct linearization structure is presented of
More informationMACSYMA PROGRAM FOR THE PAINLEVÉ TEST FOR NONLINEAR ORDINARY AND PARTIAL DIFFERENTIAL EQUATIONS
. MACSYMA PROGRAM FOR THE PAINLEVÉ TEST FOR NONLINEAR ORDINARY AND PARTIAL DIFFERENTIAL EQUATIONS Willy Hereman Mathematics Department and Center for the Mathematical Sciences University of Wisconsin at
More informationThe Construction of Alternative Modified KdV Equation in (2 + 1) Dimensions
Proceedings of Institute of Mathematics of NAS of Ukraine 00, Vol. 3, Part 1, 377 383 The Construction of Alternative Modified KdV Equation in + 1) Dimensions Kouichi TODA Department of Physics, Keio University,
More informationSymbolic Computation of Lax Pairs of Systems of Partial Difference Equations Using Consistency Around the Cube
Symbolic Computation of Lax Pairs of Systems of Partial Difference Equations Using Consistency Around the Cube Willy Hereman Department of Applied Mathematics and Statistics Colorado School of Mines Golden,
More informationarxiv:nlin/ v2 [nlin.si] 9 Oct 2002
Journal of Nonlinear Mathematical Physics Volume 9, Number 1 2002), 21 25 Letter On Integrability of Differential Constraints Arising from the Singularity Analysis S Yu SAKOVICH Institute of Physics, National
More informationHamiltonian partial differential equations and Painlevé transcendents
The 6th TIMS-OCAMI-WASEDA Joint International Workshop on Integrable Systems and Mathematical Physics March 22-26, 2014 Hamiltonian partial differential equations and Painlevé transcendents Boris DUBROVIN
More informationarxiv:math/ v1 [math.qa] 8 May 2006
Yang-Baxter maps and symmetries of integrable equations on quad-graphs Vassilios G. Papageorgiou 1,2, Anastasios G. Tongas 3 and Alexander P. Veselov 4,5 arxiv:math/0605206v1 [math.qa] 8 May 2006 1 Department
More informationBuilding Generalized Lax Integrable KdV and MKdV Equations with Spatiotemporally Varying Coefficients
Journal of Physics: Conference Series OPEN ACCESS Building Generalized Lax Integrable KdV and MKdV Equations with Spatiotemporally Varying Coefficients To cite this article: M Russo and S R Choudhury 2014
More informationContinuous limits and integrability for a semidiscrete system Zuo-nong Zhu Department of Mathematics, Shanghai Jiao Tong University, P R China
Continuous limits and integrability for a semidiscrete system Zuo-nong Zhu Department of Mathematics, Shanghai Jiao Tong University, P R China the 3th GCOE International Symposium, Tohoku University, 17-19
More informationChapter 7. Linear Algebra: Matrices, Vectors,
Chapter 7. Linear Algebra: Matrices, Vectors, Determinants. Linear Systems Linear algebra includes the theory and application of linear systems of equations, linear transformations, and eigenvalue problems.
More informationMultisolitonic solutions from a Bäcklund transformation for a parametric coupled Korteweg-de Vries system
arxiv:407.7743v3 [math-ph] 3 Jan 205 Multisolitonic solutions from a Bäcklund transformation for a parametric coupled Korteweg-de Vries system L. Cortés Vega*, A. Restuccia**, A. Sotomayor* January 5,
More informationBoussineq-Type Equations and Switching Solitons
Proceedings of Institute of Mathematics of NAS of Ukraine, Vol. 3, Part 1, 3 351 Boussineq-Type Equations and Switching Solitons Allen PARKER and John Michael DYE Department of Engineering Mathematics,
More informationMultiple-Soliton Solutions for Extended Shallow Water Wave Equations
Studies in Mathematical Sciences Vol. 1, No. 1, 2010, pp. 21-29 www.cscanada.org ISSN 1923-8444 [Print] ISSN 1923-8452 [Online] www.cscanada.net Multiple-Soliton Solutions for Extended Shallow Water Wave
More informationarxiv:nlin/ v1 [nlin.si] 13 Apr 2005
Diophantine Integrability R.G. Halburd Department of Mathematical Sciences, Loughborough University Loughborough, Leicestershire, LE11 3TU, UK (Dated: received January 1, revised March 7, 005) The heights
More informationLS.1 Review of Linear Algebra
LS. LINEAR SYSTEMS LS.1 Review of Linear Algebra In these notes, we will investigate a way of handling a linear system of ODE s directly, instead of using elimination to reduce it to a single higher-order
More informationChapter 8. Rigid transformations
Chapter 8. Rigid transformations We are about to start drawing figures in 3D. There are no built-in routines for this purpose in PostScript, and we shall have to start more or less from scratch in extending
More informationExact Solutions of The Regularized Long-Wave Equation: The Hirota Direct Method Approach to Partially Integrable Equations
Thai Journal of Mathematics Volume 5(2007) Number 2 : 273 279 www.math.science.cmu.ac.th/thaijournal Exact Solutions of The Regularized Long-Wave Equation: The Hirota Direct Method Approach to Partially
More informationQM and Angular Momentum
Chapter 5 QM and Angular Momentum 5. Angular Momentum Operators In your Introductory Quantum Mechanics (QM) course you learned about the basic properties of low spin systems. Here we want to review that
More informationSymmetries and Polynomials
Symmetries and Polynomials Aaron Landesman and Apurva Nakade June 30, 2018 Introduction In this class we ll learn how to solve a cubic. We ll also sketch how to solve a quartic. We ll explore the connections
More informationHamiltonian partial differential equations and Painlevé transcendents
Winter School on PDEs St Etienne de Tinée February 2-6, 2015 Hamiltonian partial differential equations and Painlevé transcendents Boris DUBROVIN SISSA, Trieste Cauchy problem for evolutionary PDEs with
More informationarxiv:nlin/ v2 [nlin.si] 15 Sep 2004
Integrable Mappings Related to the Extended Discrete KP Hierarchy ANDREI K. SVININ Institute of System Dynamics and Control Theory, Siberian Branch of Russian Academy of Sciences, P.O. Box 1233, 664033
More informationGALOIS GROUPS OF CUBICS AND QUARTICS (NOT IN CHARACTERISTIC 2)
GALOIS GROUPS OF CUBICS AND QUARTICS (NOT IN CHARACTERISTIC 2) KEITH CONRAD We will describe a procedure for figuring out the Galois groups of separable irreducible polynomials in degrees 3 and 4 over
More informationCONSTRUCTING SOLUTIONS TO THE ULTRA-DISCRETE PAINLEVE EQUATIONS
CONSTRUCTING SOLUTIONS TO THE ULTRA-DISCRETE PAINLEVE EQUATIONS D. Takahashi Department of Applied Mathematics and Informatics Ryukoku University Seta, Ohtsu 50-1, Japan T. Tokihiro Department of Mathematical
More informationProlongation structure for nonlinear integrable couplings of a KdV soliton hierarchy
Prolongation structure for nonlinear integrable couplings of a KdV soliton hierarchy Yu Fa-Jun School of Mathematics and Systematic Sciences, Shenyang Normal University, Shenyang 110034, China Received
More informationExact Solutions of Supersymmetric KdV-a System via Bosonization Approach
Commun. Theor. Phys. 58 1 617 6 Vol. 58, No. 5, November 15, 1 Exact Solutions of Supersymmetric KdV-a System via Bosonization Approach GAO Xiao-Nan Ô é, 1 YANG Xu-Dong Êü, and LOU Sen-Yue 1,, 1 Department
More informationM. Matrices and Linear Algebra
M. Matrices and Linear Algebra. Matrix algebra. In section D we calculated the determinants of square arrays of numbers. Such arrays are important in mathematics and its applications; they are called matrices.
More informationContents. 1 Vectors, Lines and Planes 1. 2 Gaussian Elimination Matrices Vector Spaces and Subspaces 124
Matrices Math 220 Copyright 2016 Pinaki Das This document is freely redistributable under the terms of the GNU Free Documentation License For more information, visit http://wwwgnuorg/copyleft/fdlhtml Contents
More informationChapter 12. The cross ratio Klein s Erlanger program The projective line. Math 4520, Fall 2017
Chapter 12 The cross ratio Math 4520, Fall 2017 We have studied the collineations of a projective plane, the automorphisms of the underlying field, the linear functions of Affine geometry, etc. We have
More informationHamiltonian partial differential equations and Painlevé transcendents
Winter School on PDEs St Etienne de Tinée February 2-6, 2015 Hamiltonian partial differential equations and Painlevé transcendents Boris DUBROVIN SISSA, Trieste Lecture 2 Recall: the main goal is to compare
More informationTree sets. Reinhard Diestel
1 Tree sets Reinhard Diestel Abstract We study an abstract notion of tree structure which generalizes treedecompositions of graphs and matroids. Unlike tree-decompositions, which are too closely linked
More informationSecond Order Lax Pairs of Nonlinear Partial Differential Equations with Schwarzian Forms
Second Order Lax Pairs of Nonlinear Partial Differential Equations with Schwarzian Forms Sen-yue Lou a b c, Xiao-yan Tang a b, Qing-Ping Liu b d, and T. Fukuyama e f a Department of Physics, Shanghai Jiao
More informationDiophantine integrability
Loughborough University Institutional Repository Diophantine integrability This item was submitted to Loughborough University's Institutional Repository by the/an author. Additional Information: This pre-print
More informationNUMERICAL SOLITARY WAVE INTERACTION: THE ORDER OF THE INELASTIC EFFECT
ANZIAM J. 44(2002), 95 102 NUMERICAL SOLITARY WAVE INTERACTION: THE ORDER OF THE INELASTIC EFFECT T. R. MARCHANT 1 (Received 4 April, 2000) Abstract Solitary wave interaction is examined using an extended
More informationLecture17: Generalized Solitary Waves
Lecture17: Generalized Solitary Waves Lecturer: Roger Grimshaw. Write-up: Andrew Stewart and Yiping Ma June 24, 2009 We have seen that solitary waves, either with a pulse -like profile or as the envelope
More informationThe Helically Reduced Wave Equation as a Symmetric Positive System
Utah State University DigitalCommons@USU All Physics Faculty Publications Physics 2003 The Helically Reduced Wave Equation as a Symmetric Positive System Charles G. Torre Utah State University Follow this
More informationLinearization of Mirror Systems
Journal of Nonlinear Mathematical Physics 00, Volume 9, Supplement 1, 34 4 Proceedings: Hong Kong Linearization of Mirror Systems Tat Leung YEE Department of Mathematics, The Hong Kong University of Science
More informationIntroduction to Group Theory
Chapter 10 Introduction to Group Theory Since symmetries described by groups play such an important role in modern physics, we will take a little time to introduce the basic structure (as seen by a physicist)
More informationSymbolic Computation of Lax Pairs of Integrable Nonlinear Partial Difference Equations on Quad-Graphs
Symbolic Computation of Lax Pairs of Integrable Nonlinear Partial Difference Equations on Quad-Graphs Willy Hereman & Terry Bridgman Department of Mathematical and Computer Sciences Colorado School of
More informationA variational perspective on continuum limits of ABS and lattice GD equations
A variational perspective on continuum limits of ABS and lattice GD equations Mats Vermeeren Technische Universität Berlin vermeeren@mathtu-berlinde arxiv:80855v [nlinsi] 5 Nov 08 Abstract A pluri-lagrangian
More informationFrom Continuous to Discrete Equations via Transformations
Lecture 2 From Continuous to Discrete Equations via Transformations In this chapter we will discuss how difference equations arise from transformations applied to differential equations. Thus we find that
More informationIntroduction and some preliminaries
1 Partial differential equations Introduction and some preliminaries A partial differential equation (PDE) is a relationship among partial derivatives of a function (or functions) of more than one variable.
More informationSupersymmetric Sawada-Kotera Equation: Bäcklund-Darboux Transformations and. Applications
Supersymmetric Sawada-Kotera Equation: Bäcklund-Darboux Transformations and arxiv:1802.04922v1 [nlin.si] 14 Feb 2018 Applications Hui Mao, Q. P. Liu and Lingling Xue Department of Mathematics, China University
More informationSymmetry Methods for Differential and Difference Equations. Peter Hydon
Lecture 2: How to find Lie symmetries Symmetry Methods for Differential and Difference Equations Peter Hydon University of Kent Outline 1 Reduction of order for ODEs and O Es 2 The infinitesimal generator
More informationRelativistic Collisions as Yang Baxter maps
Relativistic Collisions as Yang Baxter maps Theodoros E. Kouloukas arxiv:706.0636v2 [math-ph] 7 Sep 207 School of Mathematics, Statistics & Actuarial Science, University of Kent, UK September 9, 207 Abstract
More informationBlack Holes, Integrable Systems and Soft Hair
Ricardo Troncoso Black Holes, Integrable Systems and Soft Hair based on arxiv: 1605.04490 [hep-th] In collaboration with : A. Pérez and D. Tempo Centro de Estudios Científicos (CECs) Valdivia, Chile Introduction
More informationSIGNATURES OVER FINITE FIELDS OF GROWTH PROPERTIES FOR LATTICE EQUATIONS
SIGNATURES OVER FINITE FIELDS OF GROWTH PROPERTIES FOR LATTICE EQUATIONS JOHN A. G. ROBERTS 1, DINH T. TRAN 1,2 Abstract. We study integrable lattice equations and their perturbations over finite fields.
More informationIntroduction to Partial Differential Equations
Introduction to Partial Differential Equations Philippe B. Laval KSU Current Semester Philippe B. Laval (KSU) Key Concepts Current Semester 1 / 25 Introduction The purpose of this section is to define
More informationLinear Algebra. The analysis of many models in the social sciences reduces to the study of systems of equations.
POLI 7 - Mathematical and Statistical Foundations Prof S Saiegh Fall Lecture Notes - Class 4 October 4, Linear Algebra The analysis of many models in the social sciences reduces to the study of systems
More informationThe Solitary Wave Solutions of Zoomeron Equation
Applied Mathematical Sciences, Vol. 5, 011, no. 59, 943-949 The Solitary Wave Solutions of Zoomeron Equation Reza Abazari Deparment of Mathematics, Ardabil Branch Islamic Azad University, Ardabil, Iran
More informationDeterminant Expressions for Discrete Integrable Maps
Typeset with ps2.cls Full Paper Determinant Expressions for Discrete Integrable Maps Kiyoshi Sogo Department of Physics, School of Science, Kitasato University, Kanagawa 228-8555, Japan Explicit
More informationMTH Linear Algebra. Study Guide. Dr. Tony Yee Department of Mathematics and Information Technology The Hong Kong Institute of Education
MTH 3 Linear Algebra Study Guide Dr. Tony Yee Department of Mathematics and Information Technology The Hong Kong Institute of Education June 3, ii Contents Table of Contents iii Matrix Algebra. Real Life
More informationBreaking soliton equations and negative-order breaking soliton equations of typical and higher orders
Pramana J. Phys. (2016) 87: 68 DOI 10.1007/s12043-016-1273-z c Indian Academy of Sciences Breaking soliton equations and negative-order breaking soliton equations of typical and higher orders ABDUL-MAJID
More informationReal Analysis Prof. S.H. Kulkarni Department of Mathematics Indian Institute of Technology, Madras. Lecture - 13 Conditional Convergence
Real Analysis Prof. S.H. Kulkarni Department of Mathematics Indian Institute of Technology, Madras Lecture - 13 Conditional Convergence Now, there are a few things that are remaining in the discussion
More informationp,q H (X), H (Y ) ), where the index p has the same meaning as the
There are two Eilenberg-Moore spectral sequences that we shall consider, one for homology and the other for cohomology. In contrast with the situation for the Serre spectral sequence, for the Eilenberg-Moore
More informationRepresentation Theory
Frank Porter Ph 129b February 10, 2009 Chapter 3 Representation Theory 3.1 Exercises Solutions to Problems 1. For the Poincare group L, show that any element Λ(M,z) can be written as a product of a pure
More informationMTH 2032 Semester II
MTH 232 Semester II 2-2 Linear Algebra Reference Notes Dr. Tony Yee Department of Mathematics and Information Technology The Hong Kong Institute of Education December 28, 2 ii Contents Table of Contents
More informationLinear equations in linear algebra
Linear equations in linear algebra Samy Tindel Purdue University Differential equations and linear algebra - MA 262 Taken from Differential equations and linear algebra Pearson Collections Samy T. Linear
More informationElementary Operations and Matrices
LECTURE 4 Elementary Operations and Matrices In the last lecture we developed a procedure for simplying the set of generators for a given subspace of the form S = span F (v 1,..., v k ) := {α 1 v 1 + +
More informationarxiv:nlin/ v1 [nlin.si] 25 Sep 2006
Remarks on the conserved densities of the Camassa-Holm equation Amitava Choudhuri 1, B. Talukdar 1a and S. Ghosh 1 Department of Physics, Visva-Bharati University, Santiniketan 73135, India Patha Bhavana,
More information2. Examples of Integrable Equations
Integrable equations A.V.Mikhailov and V.V.Sokolov 1. Introduction 2. Examples of Integrable Equations 3. Examples of Lax pairs 4. Structure of Lax pairs 5. Local Symmetries, conservation laws and the
More informationOn Local Time-Dependent Symmetries of Integrable Evolution Equations
Proceedings of Institute of Mathematics of NAS of Ukraine 2000, Vol. 30, Part 1, 196 203. On Local Time-Dependent Symmetries of Integrable Evolution Equations A. SERGYEYEV Institute of Mathematics of the
More informationNumerical Analysis Lecture Notes
Numerical Analysis Lecture Notes Peter J Olver 3 Review of Matrix Algebra Vectors and matrices are essential for modern analysis of systems of equations algebrai, differential, functional, etc In this
More informationNUMERICAL METHODS FOR SOLVING NONLINEAR EVOLUTION EQUATIONS
NUMERICAL METHODS FOR SOLVING NONLINEAR EVOLUTION EQUATIONS Thiab R. Taha Computer Science Department University of Georgia Athens, GA 30602, USA USA email:thiab@cs.uga.edu Italy September 21, 2007 1 Abstract
More informationA Numerical Solution of the Lax s 7 th -order KdV Equation by Pseudospectral Method and Darvishi s Preconditioning
Int. J. Contemp. Math. Sciences, Vol. 2, 2007, no. 22, 1097-1106 A Numerical Solution of the Lax s 7 th -order KdV Equation by Pseudospectral Method and Darvishi s Preconditioning M. T. Darvishi a,, S.
More informationThe Higher Dimensional Bateman Equation and Painlevé Analysis of Nonintegrable Wave Equations
Symmetry in Nonlinear Mathematical Physics 1997, V.1, 185 192. The Higher Dimensional Bateman Equation and Painlevé Analysis of Nonintegrable Wave Equations Norbert EULER, Ove LINDBLOM, Marianna EULER
More informationLecture 10: A (Brief) Introduction to Group Theory (See Chapter 3.13 in Boas, 3rd Edition)
Lecture 0: A (Brief) Introduction to Group heory (See Chapter 3.3 in Boas, 3rd Edition) Having gained some new experience with matrices, which provide us with representations of groups, and because symmetries
More informationThe Symmetric Space for SL n (R)
The Symmetric Space for SL n (R) Rich Schwartz November 27, 2013 The purpose of these notes is to discuss the symmetric space X on which SL n (R) acts. Here, as usual, SL n (R) denotes the group of n n
More informationarxiv: v1 [nlin.si] 27 Nov 2016
A noncommutative discrete potential KdV lift S. Konstantinou-Rizos 1,2 and T. E. Kouloukas 3 arxiv:1611.08923v1 [nlin.si] 27 Nov 2016 1 Institute of Mathematical Physics and Seismodynamics, Chechen State
More informationMILNOR SEMINAR: DIFFERENTIAL FORMS AND CHERN CLASSES
MILNOR SEMINAR: DIFFERENTIAL FORMS AND CHERN CLASSES NILAY KUMAR In these lectures I want to introduce the Chern-Weil approach to characteristic classes on manifolds, and in particular, the Chern classes.
More information