Solutions by quadratures of complex Bernoulli differential equations and their quantum deformation

It is shown that the complex Bernoulli differential equations admitting the supplementary structure of a Lie-Hamilton system related to the book algebra $\mathfrak{b}_2$ can always be solved by quadratures, providing an explicit solution of the equations. In addition, considering the quantum deformation of Bernoulli equations, their canonical form is obtained and an exact solution by quadratures is deduced as well. It is further shown that the approximations of $k^{th}$-order in the deformation parameter from the quantum deformation are also integrable by quadratures, although an explicit solution cannot be obtained in general. Finally, the multidimensional quantum deformation of the book Lie-Hamilton systems is studied, showing that, in contrast to the multidimensional analogue of the undeformed system, the resulting system is coupled in a non-trivial form.


Introduction
The obtainment of explicit solutions to (ordinary or partial) differential equations, as well as finding effective criteria that determine whether such solutions can be found at all, is certainly one of the central research problems, not only for the structural theory of differential equations, but also for their applications, notably for dynamical systems.The beginning of the systematic study and classification of differential equations goes back to the XIX century, in an analytical context, where the pioneering work of Goursat, Picard or Painlevé, among others, emerged [1][2][3], complementing the purely geometrical approach initiated by Poincaré, that has evolved to constitute one of the most important current techniques [4][5][6].A somewhat different ansatz, based on group theory, and taking into account the symmetry properties, was developed by Lie, leading to the modern symmetry method, that provides a formal explanation for several of the already known solution methods [7][8][9].Indeed, the Lie symmetry method is an effective tool for either reducing the order of an equation, to find a canonical representative of the equation, and furthermore, to determine criteria that ensure that the equation is linearizable.In this context, several procedures have been proposed to linearize and solve systems of differential equations, as well as to guarantee their integrability [10,11].One auxiliary tool, whenever an explicit solution is either not available or excessively cumbersome, is to determine the existence of a (nonlinear) superposition principle, thus reducing the problem to finding a certain number of particular solutions that, in combination with some significant constants, allows to write the general solution.A characterization of this property was obtained by Lie himself in [12] using group theoretical methods, hence establishing a technique that has been extended by several authors in different directions (see e.g.[13][14][15][16][17][18][19][20] and references therein).In some circumstances, even if it is known that an equation admits a superposition principle, it may go unnoticed that, after an appropriate local diffeomorphism, the differential equation admits an explicit solution in the new coordinates.
Consider for instance a complex function w(t) = u(t) + iv(t) and the first-order equation dw dt = f (t)w + g(t)w n (1.1) for n = 2 and arbitrary real functions f (t), g(t), which is a special case of the so-called complex Bernoulli equation [21].Separating the real and imaginary parts, we are led to the first-order system du dt = f (t)u + g(t)(u 2 − v 2 ), dv dt = f (t)v + 2g(t)uv. (1.2) Solving the second equation with respect to u we are led to and the second-order equation reads as By means of the change of the dependent variable v = ξ −2 , the latter equation reduces to This equation, which is not linear, can be shown to possess a Lie point symmetry algebra isomorphic to sl(2, R) [22] (equation (1.5) actually appears as a perturbation of linear homogeneous ODE preserving an sl(2, R)-subalgebra of Noether symmetries), and by application of the Lie symmetry method [23], it can be reduced to a linear equation by means of a non-point transformation [24], from which the general solution of (1.1) is deduced.However, this procedure cannot be applied directly, without enormous computational complication, for integer values n ≥ 3, as the real and imaginary parts are not separable: (1.6) The system, however, is endowed with an additional structure that simplifies its analysis.It can be easily seen that the above system can be written in terms of the time-dependent vector field with the t-independent vector fields given by satisfying the commutator (1.9) This implies that (1.6) possesses the structure of a Lie system (see e.g.[12,25] and references therein) with a Vessiot-Guldberg algebra isomorphic to the so-called book algebra b 2 [26][27][28][29][30]. Lie systems are characterized by the remarkable property of admitting a (nonlinear) superposition principle.As will be deduced later, the system (1.6) admits a supplementary structure, namely that of a Lie-Hamilton system (LH in short), that allows us to reduce it to a linear system, from which an exact solution by quadratures will be obtained.Actually, as shown in [27], the complex Bernoulli equation (1.1) with complex coefficients admits the structure of a Lie system, but not of an LH system for arbitrary (complex) choices of the coefficient functions f, g.Hence, if an LH is known to be solvable by quadratures for a (local) system of coordinates, then, for any (local) diffeomorphism, the transformed system will also be solvable by quadratures.This fact will enable us to obtain an exact solution of (1.6) for arbitrary values of n and real coefficients f (t), g(t).
LH systems constitute a natural generalization of Lie systems, with the salient property of being related to dynamical systems through a Hamiltonian [25,26].Precisely this fact provides a systematic procedure to compute the constants of the motion needed for finding a superposition principle, using the so-called coalgebra formalism, a technique developed specifically in the context of (quantum) Hamiltonian systems [31], hence valid for LH systems, but not applicable to generic Lie systems.For this reason, LH systems have been studied focusing mainly on their coalgebra symmetry [28,29].However, depending on the specific realization of the LH system in terms of Hamiltonian vector fields, suitable coordinate frames may be found that allow us either to linearize or to integrate directly the system, making the explicit construction of superposition rules superfluous.First results in this new direction, that amounts to analyze the equivalent realizations by Hamiltonian vector fields associated to a given Lie algebra, were obtained in [30] for some particular types of LH systems based on b 2 .This work is structured as follows.In Section 2, we review the general properties of LH systems associated to the Lie algebra b 2 and reconsider their exact integrability.As a new result, we determine the most general form of first-order systems equivalent to these b 2 -LH systems by local diffeomorphisms and, using a basis different from that considered in the classification of LH systems [26], we obtain explicit solutions.In Section 3, we focus on complex Bernoulli equations admitting the structure of an LH system, and compute their explicit solution, completing and expanding the results of [29], where no exact solutions were provided.Section 4 deals with the quantum deformation of b 2 -LH systems [28][29][30], for which a solution by quadratures is computed, comparing the results with the explicit solvability of the approximations of k th -order in the quantum deformation parameter, which have not been considered previously in the literature.The latter general results for deformed b 2 -LH systems are applied in Section 5 to the construction of the corresponding deformed complex Bernoulli differential equations and, furthermore, to the obtention of their exact solution.In Section 6 we show that, for multidimensional quantum deformations, the structure of the resulting system is radically different from the undeformed case, although even in this case, the resulting equations can still be solved by quadratures.
Finally, in Section 7 we draw some conclusions and comment on possible extensions or generalizations of our results to the other isomorphism classes of LH systems, as well as their potential use in real-world applications.

Lie-Hamilton systems on the book algebra revisited
We briefly review the main features of LH systems based on the two-dimensional book Lie algebra b 2 (details concerning the general formalism of LH system can be found in [26,27,29,30] and references therein).Consider the basis b 2 = span{v A , v B } with Lie bracket (2.1) The generator v A can be seen as a dilation operator, while v B can be interpreted as a translation operator.As a common property to any Lie algebra [32,33], b 2 can be endowed with a trivial Hopf algebra structure through the (primitive) coproduct ∆ : This map further defines an algebra homomorphism for arbitrary Lie algebra A, satisfying the coassociativity constraint Hence, the pair (A, ∆) defines a coalgebra.
We now consider the symplectic representation D of b 2 defined by (see e.g.[29]) in terms of the canonical variables (x, y) and with respect to the canonical symplectic form The Hamiltonian functions (2.4) satisfy the following Poisson bracket with respect to ω: A realization of the b 2 -generators in terms of vector fields with Cartesian coordinates (x, y) ∈ R 2 is directly deduced from (2.4), using the inner product associated to ω: from which the vector fields result.It is straightforward to verify that the invariance condition is satisfied.The Hamiltonian functions (2.4) and the vector fields (2.8) determine a t-dependent Hamiltonian (respectively a t-dependent vector field) depending on two real arbitrary parameters b A (t) and b B (t) and given by (2.10) Either of these expressions leads to the same (linear) system of non-autonomous ODEs on R 2 , namely This system is easily solved by quadratures, with its exact solution being given by where c 1 and c 2 are integration constants determined by the initial conditions and a is an appropriate real number that ensures the existence of the integrals over the compact interval [a, t].
The relations (2.8) clearly show that (2.11) determines a Lie system [12,16,25] with associated Vessiot-Guldberg Lie algebra isomorphic to b 2 .The existence of Hamiltonian vector fields with respect to a Poisson structure implies that the system has a richer geometrical structure, namely that of an LH system (see [26,27] and references therein).
In spite of the simplicity of the equations (2.11) in Cartesian coordinates, LH systems associated to the book algebra b 2 can appear in disguised form, due to an inappropriate choice of coordinates, which eventually makes it difficult to recognize them as LH systems, as well as to find a suitable integration strategy.
In the new coordinates, we have that where, as a shorthand notation we have defined And the Hamiltonian functions (2.4) are transformed as follows: A routine computation shows that the transformed Hamiltonian vector fields are given by obviously satisfying the invariance condition (2.9) and the commutator (2.18) Independently on the particular shape of the resulting functions, the general solution of the system (2.18) can be obtained explicitly from (2.12), using that Φ is (locally) invertible.Supposed that ), the solution of (2.18) is given in terms of the integrals in (2.12) as It should be mentioned that, following the classification of planar LH systems [26], the Lie algebra b 2 admits only one equivalence class of realizations as vector fields on the real plane [34].This means that any first-order system that possesses the structure of a Lie system with Vessiot-Guldberg algebra isomorphic to b 2 is locally diffeomorphic to the system (2.11), and hence integrable by quadratures.This result can be stated as follows: Proposition 2.1.Any first-order system such that the associated t-dependent vector field admits the decomposition for some functions f (t), g(t) and vector fields X 1 , X 2 satisfying with λ, µ constants, can be solved by quadratures.
Proof.It suffices to verify that the vector fields X 1 , X 2 span a Lie algebra isomorphic to b 2 , the rest following from the previous discussion.If either λ = 0 or µ = 0, a scaling transformation shows that X 1 , X 2 satisfy the commutator (2.8), hence the Lie algebra is isomorphic to b 2 .On the contrary, if λµ = 0, considering the new vector field showing that the Lie algebra is isomorphic to b 2 .As Lie system (indeed, as LH system), (2.20) can always be rewritten as where We observe that, within the classification of LH systems on R 2 of [26], the book algebra b 2 corresponds to the class I r=1 14A ≃ R ⋉ R, although expressed in a different basis that is computationally more cumbersome than that considered in this work.Although, in general, it is a routine task to verify whether a first-order system (2.20) admits b 2 as Vessiot-Guldberg algebra, the effective determination of the local diffeomorphism is far from being trivial.It is precisely the additional structure of an LH system with the symplectic form and the invariance condition which allows us to systematically compute the appropriate change of coordinates, and hence to explicitly solve the system.
As an illustrative example, let us consider the system (1.2) from this perspective, instead of solving or transforming the second-order ODE (1.4).Clearly, the system is given by X = f (t)X 1 + g(t)X 2 , with showing that the system (1.2) is locally diffeomorphic to (2.11).The noncanonical symplectic form in the (u, v) coordinates is given by from which we deduce the Hamiltonian functions ) It is straightforward to verify that, for the parameter functions, we have the identities f (t) = b A (t) and g(t) = b B (t).Now, the comparison of the Hamiltonian and symplectic forms in both systems of coordinates allows us to compute easily the (local) change of coordinates, which is given by (2.29) Inserting (2.12) into the preceding expressions provides the exact solution of (1.2) in the coordinates u, v, without requiring further computation.This shows that, at least for the case of b 2 , the main characteristics of LH systems, namely the coalgebra formalism for the construction of constants of the motion and the obtainment of a superposition principle, are not required, as the system can be solved explicitly by quadratures.It should be observed that, as the Lie algebra b 2 has no nontrivial Casimir invariant, the superposition principles considered in [29] were deduced considering an embedding of the book algebra b 2 into the oscillator algebra h 4 .
In this framework, it is worth mentioning that (systems of) differential equations that correspond to LH systems based on b 2 have been considered in various contexts, such as the (generalised) Buchdahl equations arising in the study of relativistic fluids [35,36] or some particular Lotka-Volterra systems with t-dependent coefficients; they have been shown to be LH systems in [26,27].In addition, b 2 -LH systems that generalize time-dependent epidemic models with stochastic fluctuations have recently been developed in [30].Nevertheless, an explicit diffeomorphism has only been deduced in the latter work, as well as in [29] for the complex Bernoulli equations here studied, by deriving the corresponding change of variables.

Exact solutions of complex Bernoulli equations
As mentioned in the introduction, our aim is to apply Proposition 2.1 to the complex Bernoulli differential equations with t-dependent real coefficients, in order to obtain an explicit solution for any value of n.
Let us consider the differential equation (1.1) where w is a complex function, n is a real number and a 1 (t), a 2 (t) are arbitrary real t-dependent parameters.We introduce polar coordinates (r, θ) in the form w = re iθ , thus finding that the differential equation (3.1) leads to a first-order coupled system: This system was already shown in [27,29] to be an LH system, although its exact solution was not provided.According to the results in the preceding section, the system (3.2) is linearizable, i.e., it can be reduced via a change of variables to a linear system, although it is not immediately obvious how to find the appropriate new variables.It is precisely the LH formalism which provides a direct and systematic procedure to find such a reduction and hence, to solve the system explicitly.To this extent, we first reformulate the equations in terms of the t-dependent vector field with and commutator showing that Y determines a Lie system with Vessiot-Guldberg Lie algebra isomorphic to b 2 .A compatible symplectic form ω = f (r, θ)dr ∧ dθ is obtained imposing the relation (2.9), from which we deduce that The associated Hamiltonian functions h 1 , h 2 (see relation (2.7)) of the LH system are given by while the corresponding Poisson bracket with respect to the symplectic form (3.6) reads As both LH systems (2.11) and (3.2) are based on the same Lie algebra b 2 , there must exist a diffeomorphism between the expressions (2.4) with the canonical symplectic form (2.5) and variables (x, y) and (3.7) with the non-canonical symplectic form (3.6) and the variables (r, θ).This change of variables is now easily found comparing the expressions for the Hamiltonian functions and symplectic form in both systems of coordinates, leading to the explicit expression (3.9) The relationship between the corresponding vector fields, Hamiltonian functions and parameters is easily seen to be From these relations, an explicit solution of the Bernoulli equation (3.2) is obtained by merely introducing (2.12) into (3.9),where a 1 (t), a 2 (t) are arbitrary functions (and n / ∈ {0, 1}): It can be routinely checked that equations (3.2) are satisfied, observing that with the choice of the negative square root following from the relation according to (3.9).
Using the same argumentation as above, it is not difficult to verify that the differential equation with w complex, can also be solved by quadratures.With the same change of coordinates w = re iθ , the associated first-order system corresponding to the real and imaginary parts is given by In this case, the corresponding Hamiltonian vector fields are with commutator [Z 1 , Z 2 ] = (n − 1)Z 2 .Indeed, the vector fields Y 2 and Z 2 are related through the orthogonal transformation from which the assertion follows at once.

Deformed Lie-Hamilton systems from the quantum book algebra
Among the most interesting properties of LH systems, that distinguish them clearly from classical Lie systems, we enumerate the possibility of combining them with quantum groups, in order to obtain more general systems of differential equations.Even if such deformed systems are no longer described in terms of Lie algebras, we can still derive "deformed superposition principles", using the formalism of Poisson-Hopf deformations of LH systems (see [28,29] and references therein).In this section, the general solution by quadratures of the quantum deformations of b 2 -systems is obtained, with the results being particularized to the deformed complex Bernoulli differential equations, for which a deformed superposition principle was previously studied in [29].
We introduce the coboundary quantum deformation of the book Lie algebra b 2 (2.1) coming from the classical r-matrix which is a solution of the classical Yang-Baxter equation, and where z is the quantum deformation parameter such that q = e z .A structure of Lie bialgebra is then determined through the cocommutator map δ : b 2 → b 2 ∧ b 2 obtained from the classical r-matrix (4.1) as (see e.g.[33]) This element is just the skew-symmetric part of the first-order term ∆ 1 in the quantum deformation parameter z of the full coproduct ∆ z , that is, where ∆ 0 (v i ) = v i ⊗1+1⊗v i is the primitive (non-deformed) coproduct (2.2) and σ is the flip operator: For the particular case of b 2 , the quantum book algebra U z (b 2 ) ≡ b z,2 is defined by the following deformed coproduct, fulfilling the coassociativity property (2.3): The compatible deformed commutation relation is thus given by A deformed symplectic representation D z of b z,2 in terms of the canonical variables (x, y) used in Section 2 and possessing the same symplectic form (2.5) reads as closing on the following deformed Poisson bracket with respect to ω: Using the relation (2.7), we easily obtain the corresponding vector fields as (4.9) These vector fields span a smooth distribution, in the sense of Stefan-Sussmann [37,38], through the commutator [X z,A , X z,B ] = e zx X z,B , ( and satisfy the invariance condition under the Lie derivative (2.9) and the symplectic form ω (2.5) .
This construction leads to a deformed t-dependent Hamiltonian and a t-dependent vector field depending on two real arbitrary parameters b A (t) and b B (t) as follows (compare with (2.10)): giving rise to the following first-order non-autonomous system on R 2 : This system generalizes (2.11), in the sense that the expressions (4.7)-(4.12)reduce to (2.4)-(2.11)for lim z → 0. Indeed, the deformed equation can be transformed into another independent on the deformation parameter z by means of the change of variables leading to These equations can be seen as the general "canonical" system of differential equations for the set of deformed LH systems based in the quantum book algebra b z,2 , as given by (4.5).As the first of the equations is separable, the system can also be solved by quadratures.After a short computation, we arrive at the solution where and γ(t) is the same as given in (2.12).Equation (4.15) encompasses the solution of the quantum deformation for all LH systems based on the Lie algebra b 2 , provided that a proper diffeomorphism (change of variables) is found.We hence conclude that any quantum deformation of a b 2 -LH system can always be solved by quadratures, and is moreover explicitly integrable.
In the original coordinates x, y, the exact solution is then given by Observe that the non-deformed limit z → 0 is well defined, reducing to the exact solution (2.12) of the system (2.11).

On the approximations of order k in the deformation parameter
In a certain sense, the presence of the quantum deformation parameter z can be regarded as the introduction of a perturbation into the initial LH system (2.11), in such a manner that a nonlinear interaction or coupling between the variables (x, y) in the deformed LH system (4.12)arises through the term e zx y.This fact can be clearly appreciated by considering a power series expansion in z and truncating it at a certain order.The resulting ODEs can be seen as an approximation to the deformed system.
For the first-order, we obtain the system which holds for a small value of z.In this approximation, we find that z introduces a quadratic term x 2 into the first equation, becoming a (real) Bernoulli equation, while in the second, we get a nonlinear interaction term xy (this does not alter the fact that the equation is linear in y).The change of variables x = 1 zu (4.19) reduces the Bernoulli equation to the linear equation from which the solution by quadratures of the system follows at once.After a short algebraic manipulation, we arrive at the general solution for linear approximation in z, in the coordinates x, y:  In the quadratic approximation, we get the system The first equation is clearly of Abel type [39], and can be easily seen to admit the integrating factor from which the first integral is obtained.In general, even if the last integral can be explicitly solved, the solution of the ODE (hence, the system) cannot be written in explicit form.The same pattern holds for approximations at the order k ≥ 2, leading to the separable equation As the polynomial is irreducible over the rationals, having at most one real non-rational root, the application of the method of partial fractions or the Hermite method is not of great use in simplifying the integral dx xP k (x) , without numerical approximations or integration by series.In any case, as the genus of the algebraic curve y − xP k (x) = 0 is zero for any k ≥ 1, the integral is always rational, hence an elementary function, even if the solution cannot be presented in closed explicit form [40,41].
We conclude that, independently on the integral t a b A (s)ds, the system obtained from the k th -order approximation cannot be solved explicitly, although it is integrable by quadratures.

Deformed complex Bernoulli differential equations
In order to construct the quantum deformation of the complex Bernoulli differential equations along the lines described in Section 3, we keep the same symplectic form (3.6), apply the change of variables (3.9) to the deformed Hamiltonian functions (4.7) and vector fields (4.9) in the variables (x, y), and maintain the relations (3.10).The deformed Hamiltonian functions in variables (r, θ) turn out to be such that with respect to the non-canonical symplectic from (3.6).The corresponding deformed vector fields Y z,i , related with (5.1) by means of (2.7), are given by that satisfy the deformed commutator again corresponding to a distribution in the Stefan-Sussmann sense.We next consider the deformed t-dependent Hamiltonian and vector field (4.11) expressed in terms of the polar variables (r, θ) given in (3.9), thus finding that the deformed Bernoulli system of differential equations is given by (see e.g.[29]) Taking into account the general solution (4.17) for the deformed LH systems with a quantum b 2symmetry in (4.12), the change of variables (3.9), along with the relations (3.10), provide the exact solution for the deformed complex Bernoulli equations (5.5) for arbitrary values of the parameters a 1 (t), a 2 (t) and coefficients n / ∈ {0, 1}: t a a 1 (s)ds. ( It can be seen that the system (5.5) is properly satisfied.For the corresponding calculations observe that (again, only one of the square roots is admissible) which is a consequence of the relation (3.13).
Finally, we note that, as expected, all the above expressions reduce to those given in Section 3 under the non-deformed limit z → 0.

Higher dimensional quantum deformations of b 2 -LH systems
So far, we have revisited the exact solution for generic b 2 -LH systems previously given in [30], and applied this result to obtain the general solution for the complex Bernoulli equation [27,29].In addition, we have also derived the general solution for the quantum deformation of b 2 -LH systems and, as a byproduct, deduced the exact solution for the deformed complex Bernoulli equation considered in [29].These results correspond to systems of differential equations in two variables, with either the Cartesian (x, y) or polar reference (r, θ).One of the remarkable features of the underlying (deformed) coalgebra structure of LH systems, determined by the coproduct map, is that it allows us to generalize them to higher dimensions [27,29].In this section, we explicitly present the two-dimensional counterpart of the previous results and analyze the strong differences between classical and quantum deformed models.The very same features will appear in any higher dimension.
Let us consider the primitive (non-deformed) coproduct ∆ given by (2.2) on the tensor product R 2 ⊗ R 2 .We construct Hamiltonian functions h (2) i on (R 2 ) 2 with Cartesian coordinates (x 1 , y 1 , x 2 , y 2 ) using the symplectic realization D in (2.4) (see e.g.[29] for details) by means of the prescription h (2) which obviously satisfies the Lie bracket (2.6) with respect to the symplectic form The associated Poisson bracket is then given by h A , h B , ( providing the two-dimensional Hamiltonian (compare with (2.10)) A + b B (t)h (2) From the point of view of differential equations, the system resulting from the latter Hamiltonian merely consists of two copies of the non-autonomous ODEs (2.11) on (R 2 ) 2 , with the canonical variables (x 1 , y 1 ) and (x 2 , y 2 ).This property is common to all Lie systems, independently on the dimension of the symplectic representation, and implies that the general solution of such multidimensional systems is immediately obtained from the solution (2.12) given for the one-dimensional case.
The introduction of a quantum deformation, on the contrary, changes qualitatively the structure of the deformed system.The quantum deformation of b 2 , with a deformed coproduct ∆ z defined in (4.5), as well as the realization D z given in (4.7), leads to a "two-particle" symplectic realization with the following Hamiltonian functions (see [29] for details): (6.5) With respect to the symplectic form (6.2), we obtain the bracket z,B ω = − 1 − e −zh (2)   z,B z (6.6) and deformed Hamitonian (compare with Eq. (4.11)) (6.7) The corresponding equations of the motion are thus given by (6.8) Clearly, the limit for z → 0 leads to the undeformed (and uncoupled) 2-dimensional system.The main feature of the quantum deformed system is that it does not longer consist of two copies of the initial one-dimensional deformed system (4.12), but adds an additional nontrivial coupling of the variables.In this sense, the resulting differential equations (6.8) provide a different system that must be solved independently.This fact is a direct consequence of the non-trivial coproduct ∆ z in (4.5), which conveys interacting terms from its structure, being a general property for quantum deformations of LH systems (see [29] for further properties).
As expected, by means of the change of variables x k (t) = − 1 z ln u k (t), the system (6.8) can be transformed into a linear system that can again be easily solved by quadratures.We observe, in particular, that (6.9) no longer depends on the deformation parameter z.In a certain sense, we can be consider this system as the prototype of a "canonical form" for the quantum deformation in z.Notice, in particular, that the first and third equations (second and fourth equations, respectively) are not copies, but differ in a genuine deformation term.
We conclude that, if one looks for higher dimensional Bernoulli equations in the non-deformed case, one will merely find a system consisting of copies of the undeformed ODE (see Section 3), while the deformed Bernoulli equations considered in Section 5 lead to new systems for each of the considered dimensions.This implies, in particular, that if an LH system can be solved explicitly, its multidimensional version will have the same property.Whether this holds for the quantum deformed systems, is still an open question.

Conclusions and final remarks
In this work we have analyzed the exact integrability of LH systems, enlarging previous work begun in [30].Although LH systems (and their quantum deformations) are primarily characterized by the possibility of obtaining a nonlinear (deformed) superposition principle [27,29], it has been shown that, for the case of the Lie algebra b 2 , the geometric formalism can be used to obtain exact solutions of firstorder systems, both at the classical and quantum deformed level, whenever an appropriate system of coordinates is found.Using a suitable local diffeomorphism, a realization by Hamiltonian vector fields is found that enables us to reduce the differential equations to a linear system, from which an explicit solution is obtained by quadratures.The generic form of first-order systems (locally) diffeomorphic to b 2 -LH systems is found using the equivalent reformulation of the system in terms of Hamiltonian functions and the compatibility condition with the canonical symplectic form.In general, it is a routine task to check whether a given system of first-order ordinary differential equations defines a Lie system, and whether it admits a compatible symplectic structure that transforms it into an LH one.If such a system is known to admit a (local) system of coordinates that allows an integration by quadratures of the resulting differential equations, the approach proposed in this work may be more adequate than a direct analysis of the system by other methods, at least computationally, as the solution in arbitrary coordinates can be deduced directly applying the inverse of a local diffeomorphism.The argument remains valid for quantum deformations, even if formally they do not define an LH system any more, and where the deformation parameter can be seen as a perturbation term of the initial system that makes it possible to determine the approximation of k th -order in the deformation parameter.The procedure has been applied to obtain exact solutions of the so-called complex Bernoulli equations, whenever they form an LH system.Besides the quantum deformation, the approximation of k th -order in the deformation parameter has been studied, showing that it can be integrated by quadratures, although an explicit solution is no longer possible for orders k ≥ 3.These results may be of interest for solving numerically these equations, as the integrals can always be approximated with standard methods.A question that arises naturally in this context is whether the resulting numerical approximation improves or is computationally more efficient than well-established techniques, such as the predictor-corrector or the Heun methods [42].A comparison of both approaches is worthy to be analyzed in detail.
Concerning the multidimensional deformation, it differs considerably from the multidimensional undeformed analogue, in the sense that genuine deformation terms are obtained and prevent that the resulting system merely consists of copies of the one-dimensional equation.In this situation, illustrated for the two-dimensional case, it has also been shown that the deformation is solvable by quadratures.
Besides the complex Bernoulli equation, using Proposition 2.1, the results on b 2 -LH systems can also be used to determine exact solutions of other first-order systems strongly related to these equations.Consider for example the change of variables x = u p v −r , y = −v q u −m , Ξ = Λ(u p−m−1 v q−r−1 ), Λ = mr − pq = 0. ( With these new coordinates, the system (2.11) is transformed into (see Eq. (2.18)) where the variables are non-trivially coupled.We observe that each of these equations, when considered separately from the other, corresponds formally to a real Bernoulli equation, so that system (7.2) can be interpreted as a coupled Bernoulli system.Independently of the functions chosen for the coefficients b i (t), the system can be solved by quadratures using (2.12).Considering now the composition of the change of variables leading to the "canonical" form (4.14) of the quantum deformation of (2.11) with the change of variables defined above, we are led to the system du dt = b A (t) (r + q)v(t) r − q u(t) p Λ u(t) which no longer corresponds to a coupling of Bernoulli equations for arbitrary values of p, q, r, m, although its general solution can always be obtained using Eq.(4.17).For some particular values, the resulting system also corresponds to the non-trivial coupling of Bernoulli equations.For example, for the values p = −1, m = 2 and r = −q, the system (7.3)reduces to du dt = b A (t)u(t) + b B (t) v(t) q u(t) 3 , where both equations are still of Bernoulli type.
Potential applications or extensions of this work concern the analysis and exact solutions of secondorder ordinary differential equations admitting b 2 as Vessiot-Guldberg algebra [43], as well as their quantum deformation and multidimensional counterparts.Further, considering b 2 as a subalgebra of more general Lie algebras associated to LH systems, it is conceivable to extend the study of exact solutions to the quantum deformed systems of LH systems related to the oscillator algebra h 4 [26,27,29], either from the formal point of view, or taking into account specific applications [30].
In a general context, a problem not yet studied in detail is the systematic analysis of the solvability by quadratures, or the obtainment of criteria that ensure the existence of explicit solutions of LH systems and their (multidimensional) generalizations and deformations, as well as the generic description of systems that are (locally) equivalent to an LH system associated to a given Lie algebra.Another question to be settled is the precise identification of LH structures with systems of ordinary differential equations appearing in real-world applications, specifically with those systems possessing time-dependent coefficients, such as e.g.some specific types of Lotka-Volterra systems (see [26,27] and references therein), where the procedure may provide an alternative approach to obtain either exact or approximate solutions, once the parameters have been appropriately adjusted.Work in these various directions is currently in progress.