On the Question of the B\"acklund Transformations and Jordan Generalizations of the Second Painlev\'e Equation

We demonstrate the way to derive the second Painlev\'e equation $P_2$ and its B\"acklund transformations from the deformations of the Nonlinear Schr\"odinger equation (NLS), all the while preserving the strict invariance with respect to the Schlesinger transformations. The proposed algorithm allows for a construction of Jordan algebra-based completely integrable multiple-field generalizations of $P_2$ while also producing the corresponding B\"acklund transformations. We suggest calling such models the JP-systems. For example, a Jordan algebra $J_{_{{\rm Mat}(N,N)}}$ with the Jordan product in the form of a semi-anticommutator is shown to generate an integrable matrix generalization of $P_2$, whereas the $V_{_N}$ algebra produces a different JP-system that serves as a generalization of the Sokolov's form of a vectorial NLS.

We demonstrate the way to derive the second Painlevé equation P2 and its Bäcklund transformations from the deformations of the Nonlinear Schrödinger equation (NLS), all the while preserving the strict invariance with respect to the Schlesinger transformations. The proposed algorithm allows for a construction of Jordan algebra-based completely integrable multiple-field generalizations of P2 while also producing the corresponding Bäcklund transformations. We suggest calling such models the JP-systems. For example, a Jordan algebra J Mat(N,N ) with the Jordan product in the form of a semi-anticommutator is shown to generate an integrable matrix generalization of P2, whereas the V N algebra produces a different JP-system that serves as a generalization of the Sokolov's form of a vectorial NLS.

I. INTRODUCTION
The triumphant emergence of the Painlevé equations dates back to the very end of the nineteenth century, when Emile Picard posed the following question [1]: what kind of second order O.D.E.s contain no movable singularities except for poles? By 1910, Paul Panilevé and his student Bertrand Gambier [2][3][4] proved that there are 50 types of second order O.D.E.s with such a property. Six of them were shown to be irreducible to either elementary or classical special functions. These interesting solutions have been dubbed the Painlevé transcendents, and the corresponding six equations, the Painlevé equations I-VI. The Painlevé equations have been extensively studied as isomonodromic deformations of linear systems [5][6][7][8][9], and to this day, they remain one of the most important ingredients of the integrable systems theory and the one most shrouded in mystery. The equations arise in the problems that involve self-similar solutions of integrable hierarchies [10,11]; they helped establish a (to this day not sufficiently understood) relationship between the integrability and the problems of preservation of O.D.E.s' monodromy [12,13]; and they positively proliferate when one studies the dressing chains of discrete symmetries [14]. Additionally, besides the multitude of articles that are dedicated to the equations themselves, there are also many works on the various generalizations of said equations, running the gamut from discrete to multiple-fields matrix models.
However, the importance of the Painlevé equations is not limited to the field of mathematics, since many of them arise in a number of interesting physical problems. For example, the solutions of Painlevé III emerge in the studies of the spin-spin correlation function in the 2D Ising model [15] and occur in the scaling functions for two-dimensional polymers [16,17]; Painlevé IV is a feature in generalizations of odd superpotentials when one studies the exact nonsingular cosmological solutions on a 3D brain that interacts with five-dimensional gravity and the bulk scalar field [18]; Painlevé V is necessary in the descriptions of a density matrix of an impenetrable Bose gas [19]; and the conformal field theory was even shown to be intimately connected to Painlevé VI [20], as well as to V and III [21] (see also [22]).
In this article, we will concentrate on Painlevé II, which for brevity we will henceforth call P 2 : these polynomials were shown to crop up in general solutions of soliton P.D.E.s, from the Korteweg-De Vries (KdV) equation and NLS to Kadomtsev-Petviashvili (KP) equations (see also [33]), while in [29], it has been demonstrated that any rational solution of P 2 can be constructed via the Bäcklund transformations (23)-(25) (we will proceed to explicitly derive these transformations later in the article) from the solution to the homogeneous version of (1), i.e., the one with α = 0. The article [35] has introduced a powerful new way to derive the Painlevé equations from integrable evolution equations. Its gist was to utilize the deformations of NLS. First, one introduces an auxiliary linear evolution equation and one auxiliary spectral problem, in which the potential depends on the spectral parameter as a polynomial of either the first or second order. While the former polynomial leads to KdV, the latter one produces NLS, which is natural, since a quadratic dependence upon the spectral parameter is but a linear Zakharov-Shabat problem, albeit written in some special gauge [36]. In the next step, one performs a deformation of both the spectral problem and of the compatibility equation. Finally, one looks for the travelling wave solutions of the resulting equations, only to end up with the sought after Painlevé equation. It is a very attractive method, but it is also admittedly somewhat overwrought. And as we shall see below, it is quite possible to derive P 2 from the (trivially) deformed NLS system in a more direct and natural way; it appears that the key to this lies not in the auxiliary spectral problem, but in the famous symmetries of NLS known as Schlesinger transformations (a term we use following the seminal work [37]). Furthermore, it is these very transformations that end up serving as the Bäcklund transformations for P 2 . Even better, this approach is versatile enough to be easily generalizable to multiple-fields models. To be more precise, recall that all multiple-fields integrable NLS generalizations can be embedded into one general formalism, associated with a unital Jordan algebra [38]. For example, if one starts with a Jordan algebra J Mat(N,N ) with the Jordan product in the form of an anticommutator divided by two, it will lead to a system of matrix-valued NLS equations. Such systems are called the JS-systems (where JS stands for Jordan-Shrödinger), and since we shall soon demonstrate that our approach naturally leads to a matrix-valued analogue of P 2 (and to its Bäcklund transformations), we think it natural to introduce the term JP-systems for such equations (JP here stands for Jordan-Painlevé). While the matrix-valued Painlevés themselves have been a subject of many research papers, our method clearly delineates them as but the first and simplest exhibit in a larger JP-systems menagerie. Other types of Jordan algebras will instead produce their own distinct JP-systems.
Before we continue, let us briefly discuss the structure of this paper. Section II is dedicated to a discussion of relevant symmetries, known as the Schlesinger transformations for NLS, and to their relationship with the Toda chain equations. We delve deeper in Section III, where we study the NLS deformations, prove the invariance of the trivial deformations with respect to the Schlesinger transformations, and then derive the system of homogeneous split P 2 equations. This system is equivalent to a single fourth order O.D.E., but its order can be lowered, and it is exactly what we will do in Section IV. There, we derive the complete P 2 equation along with the corresponding Bäcklund transformations and show a simple way to generalize the results for the matrix-valued analogue of P 2 . While we are at it, we also demonstrate how a number of this equation's properties, originally proven only via rather difficult calculations, turn out to be essentially self-evident when viewed through the lens of our approach. Similarly, deriving the Bäcklund transformations for the matrix-valued P 2 is very easy and is accomplished at the end of that Section. After that, we move on to Section V, where we describe a new JP-system which serves as a generalization of a vectorial Sokolov model. In Section VI, we take a small step out of the confines of the discussion focused on Painlevé II and very briefly describe the first few steps in a possible application of our method to the task of Jordan generalizations of Painlevé IV (P 4 ). Finally, we conclude this paper in Section VII by briefly discussing some of the most straightforward and potentially interesting steps to further the subsequent development of this rich topic.

Consider a split Nonlinear Schrödinger equation (NLS):
and It is easy to see that Hence (3) and (4) are explicitly invertible Bäcklund auto-transformations for (2).

Remark 1
The article [39] has introduced the concept of the so-called Ablowitz-Kaup-Newell-Segur (AKNS) hierarchy of nonlinear integrable P.D.E.s. The hierarchy starts with a (split) NLS, then goes to the coupled modified Kortewegde Vries (cmKdV) system, and the third entry is the Lakshmanan-Porsezian-Daniel (LPD) equations. (For a detailed description of the AKNS hierarchy, see [40].) It is possible to prove by induction that all nonlinear equations belonging to the AKNS hierarchy are invariant with respect to the Schlesinger transformations (3) and (4). However, the proof is quite cumbersome and lies out of scope of this work and, as such, will be the subject of a follow-up article.
The STs are intimately related to the Toda chain equations. If we define then the transformation leads to a welcome surprise in the shape of Toda lattice: where the prime denotes the derivative with respect to the variable x. It is safe to say that the relationship between Equations (3) and (4), on the one hand, and the Toda chain (6), on the other hand, is truly one of most peculiar properties of ST.
There are many known ways to generalize the Schlesinger's approach. For example, one might jump from the (1 + 1) NLS (2) to the (1 + 2) Davey-Stewartson (DS) equations [41]. Alternatively, one might instead utilize the Jordan generalizations of NLS, introduced by Svinolupov and Yamilov [38]; this approach is known to produce a number of non-trivial structures known as JS-systems and JT-systems. However, in this article we will concentrate on the third avenue of research that couples ST with P 2 . In particular, we will demonstrate that the famous Bäcklund transformations for P 2 are exactly the Schlesinger transformations (3) and (4). We will perform it by starting out with the NLS deformations.

III. THE NLS DEFORMATIONS
The NLS deformation can be written in the following form [35]: where D −1 x denotes an indefinite integral with respect to the variable x. As has been noticed by Shabat and Yamilov, these equations are related to four of the Painlevé equations: P 2 , P 3 , P 4 , and P 5 .
Let c 1 = 0, c 2 = 0. Then the deformations (7) and (8) are trivial and reducible to the split NLS (via the shift-and Galilean point transformations). In this case, (7) and (8) are invariant with respect to ST, coupled with the point transformations. The latter can be taken care of by setting c 3 = 0. This yields the following: (7) and (8)  Consider a travelling wave solution of (7) and (8): which abides by the following conditions: in accordance with Proposition 1 (keep in mind that c 1 = c 3 = 0). After the point transformation u → −u and assuming that c 2 = 0, systems (7) and (8) reduce to which can be called a split homogeneous P 2 equation, because when u = v = w(x, 0) ≡ w, it further reduces to a special case of (1) with α = 0. Of course, the new split P 2 is invariant with respect to ST (3) and (4), which, for this case, take the form and, naturally, log u n , log v n still satisfy the Toda chain (6). At first glance, the results we have gained are less than spectacular. First of all, it is easy to see that the new STs (10) and (11) do not preserve the reduction u = v, thus creating the impression that the equation P 2 (1) might, after all, be unrelated to both (10) and (11) and the corresponding Toda chain. Secondly, the r.h.s. of (9) bears no constant terms and thus might at best be reducible to a homogeneous P 2 with α = 0. All of this is exacerbated by the fact that, in general, the equation we end up with will be of a fourth order: However, as we shall see below, the order of (12) can be reduced. Before we do that, though, we shall point out a number of very interesting properties it possess. (12) is invariant with respect to transformation u → βu, with β = const: (12) is invariant with respect to three discrete transformations, generated by the nonlinear operators P 0 , P 1 , P 2 :

Property 2 Equation
where u is a solution of 4th order Equation (12).
Property 3 The operators P 0 , P 1 , and P 2 have the following properties: where Eu = u and Remark 2 The interweaving relations (17) are nontrivial because For example, Using these properties one can construct an infinite set of exact solutions of (12), building upon some simple seminal solution.
Example 1 Let us begin by adopting a trivial solution u 0 = 1 of (12). We obtain an infinite set of progressively more complex solutions from it. Let us list some of them: x(x 9 − 30x 6 − 1400) , and so on.

IV. MATRIX P2 RELOADED
Now, let us explain how to lower the power of differential Equation (12). We begin by noticing that the first integral of (9) would be its Wrońskian, i.e., After we remove function v from (20), we end up with a new third order equation: and it is straightforward to show that any solution of (21) will serve as such for (12) as well.
We are not finished. The 3rd order equation is homogeneous with respect to combination {u, u , u , u }, which implies that the Cole-Hopf transformation q = (log u) will further reduce it to a second order O.D.E.: where α = C − 1/2. It is this equation which we will subsequently call P 2 , and for a good reason: it is reducible to the canonic form (1) via two simple substitutions: However, in order to prevent an unnecessary cluttering of our formulas by throwaway coefficients, for now we will continue working with (22) under the assumption q = q(x, α).
It is both important and interesting to note that the P i symmetries are in fact the famous Bäcklund transformations for P 2 . Since we have chosen to embrace a non-canonical form of P 2 , it is instructive to show these transformations. Here they are: The transformations (23)-(25) are well-known and well-studied, although they are usually derived from the Bäcklund transformations for KdV equations with the aid of a self-similar change of variables. We can see now that (23)-(25) serve as discrete symmetries for P 2 , i.e., the explicitly invertible Bäclund transformations for (9), which are, in turn, nothing but ST for the trivial NLS deformation. Additionally, their connection with the Toda chain becomes very straightforward.
The (24) and (25) are well-studied in the literature, so we can now safely leave them and move on to the Jordan generalizations. As a particular example, consider the following matrix equation P 2 : where W (x, α) and α = const are N × N matrices. Our goal is to derive (26) in the same manner as (22) before that.
To this end, we write a matrix equivalent of (9): where U and V are, again, N × N square matrices. The sought after analogue of ST Equations (10) and (11) (the analogous symmetries for the Davey-Stewartson equation have been studied in [41] by one of the authors) will be Note that the symmetries (28) and (29) are mutually inverse: By removing V , we end up with a matrix analogue of (12): whose order is reducible via the first integral Remark 4 One can also pick an alternative integral: U V − U V = B. However, since the system (27) is invariant with respect to inversion U ↔ V , this option simply means that B = −A.
Now, we define the matrix-valued function q as such that it satisfies the condition so after a few simple calculations (31), we obtain Of course, if A is proportional to a unitary matrix A = cE, then (33) will be identical to (22), except that q will be an N × N matrix, whereas 2α = 2c − 1 will be equal to a product of E and a number. The equation then turns into (26) via the same point transformation as the one used in the scalar case.
It has been previously pointed out by Balandin in [42] that the diagonality of A might be a necessary condition for the integrability of matrix P 2 . We now see that the non-diagonality actually leads to the integrability (for our purposes, the integrability of a nonlinear equation means the existence of an explicitly invertible Bäcklund transformation-an analogue of the Schlesinger transformations) of some non-local integro-differential matrix analogue of P 2 . However, this equation might be rewritten in a completely local form via a different substitution: which results in a different, and rather elegant, matrix O.D.E.: In a scalar case, the commutator vanishes, delivering us P 2 .
Remark 5 This work was reaching its conclusion when we learned of a paper by Adler and Sokolov [43], published in 2021. There, the authors introduced three matrix-valued integrable generalizations of P 2 , denoted by Adler and Sokolov as P 0 2 , P 1 2 , and P 1 2 . Out of these three, the first one was (34) exactly, whereas the second generalization can be derived from (34) via a specific matrix shift of an independent variable. Interestingly, the authors of [43] have introduced their versions of equations mostly voluntarily (essentially acting on the principle of "in a scalar case a term with the commutator vanishes; let's add it to the mix!"), while in our approach, (34) is a necessary outcome of previous calculations. Nevertheless, we admit that we cannot claim the honour of discovering (34); it rightfully belongs to Adler and Sokolov.
To conclude our discussion of Jordan generalizations, let us jot down the Bäcklund transform for (34). Let Then, the direct Bäcklund transformation will have the form while its inverse where

V. THE VECTORIAL PAINLEVÉ EQUATION
Now, let us take one more step and attempt to tackle a vectorial generalization of P 2 . In order to do this, we will use the model of vector NLS originally introduced by Sokolov (and studied in [38]). Let U = {u 1 (x), u 2 (x), .., u N (x)} T and V = {v 1 (x), v 2 (x), .., v N (x)} T be two N -dimensional vectors with a standard Euclidean scalar product What kind of equation are we going to call an N -component vector P 2 ? It will be a system of 2N equations, with the following vector representation: One might at first hesitate calling (37) the P 2 equation, but it is in fact a literal generalization of P 2 . For example, in the case of N = 2, (37) has the form and for N = 3, Consider (38). Solve the first two equations for v 1 , v 2 and substitute them into the remaining equation. The result will be a system of two equations of fourth order for u 1 and u 2 that we will omit here owing to its cumbersomeness. However, if we assume that u 1 = u 2 = u, both of these equations will be reducible to (12), which is just another form of P 2 . In a similar vein, the system (39) can be coaxed into producing a system of three equations of fourth order with a total of 92 terms in each; they are tamed by the conditions u 1 = u 2 = u 3 = u, contracting down to a single Equation (12).
The vector form of P 2 permits an explicit self-transformation of the kind which allows us to construct infinitely diverse families of exact solutions for (37) based on the initial trivial offering. For example, if we take an arbitrary constant N -vector w = {w 1 , w 2 , w 3 , ..., w N } T , then the simplest solutions of (37) will be of the form Applying to them our transformation (40), after the n-th iteration, we will gather the following solution: where the scalar functions v (n) = v (n) (x), u (n) = u (n) (x) are derived with the aid of the Schlesinger transformation: Here are the fruits of our labours. Applying (44) for (42) yields , etc. Thus, we obtain a chain of exact solutions u (k) , k ∈ N, which we shall recognize straight away, because we have seen them before. In fact, they are almost identical to the u k solutions (18) we procured in Section III. A closer examination reveals the following relationship between the old solutions u k and the new one: Another interesting starting point for generating new solutions would be the null-vector V ≡ 0. This choice not only converts the system (37) into a linear equation, but it also implies that every component u k of vector U satisfies the Airy equation: Two linearly independent solutions of (46) are the so-called Airy functions Ai(x) and Bi(x), which for x ∈ R can be written as [44] Ai One can therefore write general solutions of Equation (46) as a linear combination of Ai(x) and Bi(x): and applying the transformations (40) and (41) to (48) (keeping in mind that v k ≡ 0), we will end up with a rather different family of exact solutions that not only contains 2N arbitrary parameters {a k } N k=1 , {b k } N k=1 , but is also completely entirely determined via the Airy functions and their derivatives. (This is but one of many examples where the solutions of the integrable nonlinear differential equations can be written in terms of the Airy functions; for instance, in [45], the Airy functions are shown to be an important ingredient in the solution of the Cauchy problem for the Novikov-Veselov equation.) Before the conclusion, let us emphasize the fact that (37) is a JP-system, whose Jordan algebra V N is determined by the Jordan product (e, x)y + (e, y)x − (x, y)e, where e, x, and y are the elements of N -dimensional vector space; ( . , . ) is the scalar product; and e = (1, 0, 0, ..) T is a unit in the V N algebra.

VI. ON A WAY TO THE PAINLEVÉ IV
In the previous sections, we have succeeded in developing a technique for the construction of solutions for the P 2 equation, thus completing the main goal of this article. Nevertheless, the success motivates us to at least take a short glance at the possibility of extending our technique for the remaining five Panilevé equations. Thus, in this section we will briefly explain how to do that for Painlevé IV. First of all, we have to point our that NLS, whose deformations acted as a sort of a launchpad for our investigation, is tightly related to yet another famous integrable equation known as the Kadomtsev-Petviashvili equation (KP). This relationship, discovered and used in [46,47] to generate new rogue wave soliton solutions (i.e., the solutions that are localized in both space and time; originally discovered by Howell Peregrine in 1983 for NLS [48], the rogue wave solutions have been steadily cropping up in almost every field of mathematical physics, from the collapse of intrathermocline eddies in the ocean [49] to a magnetic "impacton" arising during a collision of two positon solutions in ferromagnetic nanowires [50]), can only be properly understood in the framework of the AKNS hierarchy (see Remark 1 and article [40]). For our purposes, it is important to single out the existence of an explicit relationship between the solutions of (split) NLS and KP. The latter equation can be written as a compatibility condition of two different linear equations; their shapes determine two types of dressing chains of discrete symmetries. Say we choose a first one [51]: with α 2 = ±1. Let us add a periodic boundary condition where f = f n , s = f n+1 . For a period N = 3, choose c = −2y/α and introduce new field variables g n , n = 1, 2, 3 as follows:

VII. CONCLUSIONS
In this article, we have demonstrated how the coupling of NLS deformations with the Schlesinger transformation naturally produces the Second Painlevé equation P 2 and its Bäcklund transformations. The resulting scheme is versatile enough to help generalize the JS-system formalism for the JP-systems, by which we understand the multiplefields integrable generalization of P 2 , based upon the unital Jordan algebras. This opens up a number of very promising avenues of research, including the ultimate goal of classification and description of JP-systems. Another interesting possibility lies in applying the new results to the task of constructing new discrete analogues of P 2 (cf. [53]), associated with the Jordan algebras. The novel approach developed herein is based on the invertible Bäcklund transformations, which serve as a crucial part of the JP-systems theory. This further reinforces our opinion that the process of the discretization of the JP-system will see no substantial difficulties. Now let us summarize the main results.
1. We have demonstrated how the trivial deformations of NLS, invariant with respect to the Schlesinger transformations, can be reduced to a system of two split O.D.E.s of second order that inherits this symmetry.
2. The resulting system is shown to be equivalent to a single fourth order O.D.E.; its order can be reduced, producing a familiar equation: P 2 .
3. Inherent Schlesinger symmetries appear to be nothing more but well-known explicitly invertible Bäcklund transformations for P 2 . Thus, we obtain a very simple and concise method for the derivation of both P 2 and its Bäcklund autotransformations.
4. The simplicity of the method paves the way for its generalization for multiple fields models-a problem which is traditionally considered a difficult one. We demonstrate how the new approach dispels the difficulties in establishing multicomponent P 2 generalizations associated with arbitrary unital Jordan algebras. This implies that the multicomponent integrable (in the sense of having an analogue of the Schlesinger transformation) P 2 generalizations can actually be classified, for example, by a set of structural constants of a corresponding Jordan algebra. We have called such models JP-systems.

5.
Another interesting observation that needs to be emphasized is that existence of explicitly invertible Bäcklund transformations opens a way to generalize the proposed method to the discrete Painlevé equations. Taking into account the rapid growth of interest in the discretization of Painlevé equations (see [53][54][55][56]), it would not be unreasonable to expect the said method to be useful to the research.
One last question must be addressed before we wrap up. In this article, we have developed a method of generalization for the Second Painlevé equation. What about the remaining five? In the last Section, we briefly touched upon Painlevé IV and how it can be naturally introduced into our considerations via the fascinating relationship between NLS and KP equations. As for the rest, at this stage we can only surmise that it might be possible to construct a similar scheme for at last some of them, since these equations also posses the Schlesinger symmetries and the chains of discrete symmetries [52]. This, however, is of course a matter for another time and another article. For now, let us simply conclude by stating out humble hope that the approach developed in this article might be of use in both the theory of Painlevé equations and in mathematical physics in general.