Consider an operator
with smooth coefficients and look for a factorization
Let us write down the equations on p i {\displaystyle p_{i}} explicitly, keeping in mind the rule of left composition, i.e. that
Then in all cases
where the notation L = p 1 ∂ x + p 2 ∂ y {\displaystyle {\mathcal {L}}=p_{1}\partial _{x}+p_{2}\partial _{y}} is used.
Without loss of generality, a 20 ≠ 0 , {\displaystyle a_{20}\neq 0,} i.e. p 1 ≠ 0 , {\displaystyle p_{1}\neq 0,} and it can be taken as 1, p 1 = 1. {\displaystyle p_{1}=1.} Now solution of the system of 6 equations on the variables
can be found in three steps.
At the first step, the roots of a quadratic polynomial have to be found.
At the second step, a linear system of two algebraic equations has to be solved.
At the third step, one algebraic condition has to be checked.
Step 1. Variables
can be found from the first three equations,
The (possible) solutions are then the functions of the roots of a quadratic polynomial:
Let ω {\displaystyle \omega } be a root of the polynomial P 2 , {\displaystyle {\mathcal {P}}_{2},} then
Step 2. Substitution of the results obtained at the first step, into the next two equations
yields linear system of two algebraic equations:
In particularly, if the root ω {\displaystyle \omega } is simple, i.e.
equations have the unique solution:
At this step, for each root of the polynomial P 2 {\displaystyle {\mathcal {P}}_{2}} a corresponding set of coefficients p j {\displaystyle p_{j}} is computed.
Step 3. Check factorization condition (which is the last of the initial 6 equations)
written in the known variables p j {\displaystyle p_{j}} and ω {\displaystyle \omega } ):
If
the operator A 2 {\displaystyle {\mathcal {A}}_{2}} is factorizable and explicit form for the factorization coefficients p j {\displaystyle p_{j}} is given above.
Similar to the case of the operator A 2 , {\displaystyle {\mathcal {A}}_{2},} the conditions of factorization are described by the following system:
with L = p 1 ∂ x + p 2 ∂ y , {\displaystyle {\mathcal {L}}=p_{1}\partial _{x}+p_{2}\partial _{y},} and again a 30 ≠ 0 , {\displaystyle a_{30}\neq 0,} i.e. p 1 = 1 , {\displaystyle p_{1}=1,} and three-step procedure yields:
At the first step, the roots of a cubic polynomial
have to be found. Again ω {\displaystyle \omega } denotes a root and first four coefficients are
At the second step, a linear system of three algebraic equations has to be solved:
At the third step, two algebraic conditions have to be checked.
Definition The operators A {\displaystyle {\mathcal {A}}} , A ~ {\displaystyle {\tilde {\mathcal {A}}}} are called equivalent if there is a gauge transformation that takes one to the other:
BK-factorization is then pure algebraic procedure which allows to construct explicitly a factorization of an arbitrary order LPDO A ~ {\displaystyle {\tilde {\mathcal {A}}}} in the form
with first-order operator L = ∂ x − ω ∂ y + p {\displaystyle {\mathcal {L}}=\partial _{x}-\omega \partial _{y}+p} where ω {\displaystyle \omega } is an arbitrary simple root of the characteristic polynomial
Factorization is possible then for each simple root ω ~ {\displaystyle {\tilde {\omega }}} iff
for n = 2 → l 2 = 0 , {\displaystyle n=2\ \ \rightarrow l_{2}=0,}
for n = 3 → l 3 = 0 , l 31 = 0 , {\displaystyle n=3\ \ \rightarrow l_{3}=0,l_{31}=0,}
for n = 4 → l 4 = 0 , l 41 = 0 , l 42 = 0 , {\displaystyle n=4\ \ \rightarrow l_{4}=0,l_{41}=0,l_{42}=0,}
and so on. All functions l 2 , l 3 , l 31 , l 4 , l 41 , l 42 , . . . {\displaystyle l_{2},l_{3},l_{31},l_{4},l_{41},\ \ l_{42},...} are known functions, for instance,
and so on.
Theorem All functions
are invariants under gauge transformations.
Definition Invariants l 2 = a 00 − L ( p 6 ) + p 3 p 6 , l 3 = a 00 − L ( p 9 ) + p 3 p 9 , l 31 , . . . . . {\displaystyle l_{2}=a_{00}-{\mathcal {L}}(p_{6})+p_{3}p_{6},l_{3}=a_{00}-{\mathcal {L}}(p_{9})+p_{3}p_{9},l_{31},.....} are called generalized invariants of a bivariate operator of arbitrary order.
In particular case of the bivariate hyperbolic operator its generalized invariants coincide with Laplace invariants (see Laplace invariant).
Corollary If an operator A ~ {\displaystyle {\tilde {\mathcal {A}}}} is factorizable, then all operators equivalent to it, are also factorizable.
Equivalent operators are easy to compute:
and so on. Some example are given below:
Factorization of an operator is the first step on the way of solving corresponding equation. But for solution we need right factors and BK-factorization constructs left factors which are easy to construct. On the other hand, the existence of a certain right factor of a LPDO is equivalent to the existence of a corresponding left factor of the transpose of that operator.
Definition The transpose A t {\displaystyle {\mathcal {A}}^{t}} of an operator A = ∑ a α ∂ α , ∂ α = ∂ 1 α 1 ⋯ ∂ n α n . {\displaystyle {\mathcal {A}}=\sum a_{\alpha }\partial ^{\alpha },\qquad \partial ^{\alpha }=\partial _{1}^{\alpha _{1}}\cdots \partial _{n}^{\alpha _{n}}.} is defined as A t u = ∑ ( − 1 ) | α | ∂ α ( a α u ) . {\displaystyle {\mathcal {A}}^{t}u=\sum (-1)^{|\alpha |}\partial ^{\alpha }(a_{\alpha }u).} and the identity ∂ γ ( u v ) = ∑ ( γ α ) ∂ α u , ∂ γ − α v {\displaystyle \partial ^{\gamma }(uv)=\sum {\binom {\gamma }{\alpha }}\partial ^{\alpha }u,\partial ^{\gamma -\alpha }v} implies that A t = ∑ ( − 1 ) | α + β | ( α + β α ) ( ∂ β a α + β ) ∂ α . {\displaystyle {\mathcal {A}}^{t}=\sum (-1)^{|\alpha +\beta |}{\binom {\alpha +\beta }{\alpha }}(\partial ^{\beta }a_{\alpha +\beta })\partial ^{\alpha }.}
Now the coefficients are
A t = ∑ a ~ α ∂ α , {\displaystyle {\mathcal {A}}^{t}=\sum {\tilde {a}}_{\alpha }\partial ^{\alpha },} a ~ α = ∑ ( − 1 ) | α + β | ( α + β α ) ∂ β ( a α + β ) . {\displaystyle {\tilde {a}}_{\alpha }=\sum (-1)^{|\alpha +\beta |}{\binom {\alpha +\beta }{\alpha }}\partial ^{\beta }(a_{\alpha +\beta }).}
with a standard convention for binomial coefficients in several variables (see Binomial coefficient), e.g. in two variables
In particular, for the operator A 2 {\displaystyle {\mathcal {A}}_{2}} the coefficients are a ~ j k = a j k , j + k = 2 ; a ~ 10 = − a 10 + 2 ∂ x a 20 + ∂ y a 11 , a ~ 01 = − a 01 + ∂ x a 11 + 2 ∂ y a 02 , {\displaystyle {\tilde {a}}_{jk}=a_{jk},\quad j+k=2;{\tilde {a}}_{10}=-a_{10}+2\partial _{x}a_{20}+\partial _{y}a_{11},{\tilde {a}}_{01}=-a_{01}+\partial _{x}a_{11}+2\partial _{y}a_{02},}
For instance, the operator
is factorizable as
and its transpose A 1 t {\displaystyle {\mathcal {A}}_{1}^{t}} is factorizable then as [ . . . ] [ ∂ x − ∂ y + 1 2 ( y + x ) ] . {\displaystyle {\big [}...{\big ]}\,{\big [}\partial _{x}-\partial _{y}+{\tfrac {1}{2}}(y+x){\big ]}.}
Weiss (1986) ↩
R. Beals, E. Kartashova. Constructively factoring linear partial differential operators in two variables. Theor. Math. Phys. 145(2), pp. 1510-1523 (2005) https://doi.org/10.1007%2Fs11232-005-0178-7 ↩
E. Kartashova. A Hierarchy of Generalized Invariants for Linear Partial Differential Operators. Theor. Math. Phys. 147(3), pp. 839-846 (2006) https://doi.org/10.1007%2Fs11232-006-0079-4 ↩
E. Kartashova, O. Rudenko. Invariant Form of BK-factorization and its Applications. Proc. GIFT-2006, pp.225-241, Eds.: J. Calmet, R. W. Tucker, Karlsruhe University Press (2006); arXiv https://arxiv.org/abs/math-ph/0607040/ ↩