In mathematics, differential algebra is, broadly speaking, the area of mathematics consisting in the study of differential equations and differential operators as algebraic objects in view of deriving properties of differential equations and operators without computing the solutions, similarly as polynomial algebras are used for the study of algebraic varieties, which are solution sets of systems of polynomial equations. Weyl algebras and Lie algebras may be considered as belonging to differential algebra.
More specifically, differential algebra refers to the theory introduced by Joseph Ritt in 1950, in which differential rings, differential fields, and differential algebras are rings, fields, and algebras equipped with finitely many derivations.
A natural example of a differential field is the field of rational functions in one variable over the complex numbers,
C(t),
t.
Joseph Ritt developed differential algebra because he viewed attempts to reduce systems of differential equations to various canonical forms as an unsatisfactory approach. However, the success of algebraic elimination methods and algebraic manifold theory motivated Ritt to consider a similar approach for differential equations. His efforts led to an initial paper Manifolds Of Functions Defined By Systems Of Algebraic Differential Equations and 2 books, Differential Equations From The Algebraic Standpoint and Differential Algebra. Ellis Kolchin, Ritt's student, advanced this field and published Differential Algebra And Algebraic Groups.
A derivation on a ring is a function
\partial:R\toR
\partial(r1r2)=(\partialr1)r2+r1(\partialr2)
r1
r2
R.
A derivation is linear over the integers since these identities imply
\partial(0)=\partial(1)=0
\partial(-r)=-\partial(r).
R
r\inR.
A differential field is a differential ring that is also a field. A differential algebra
A
K
K
K
A
K.
K
K
A Witt algebra is a differential ring that contains the field
\Q
\Q,
\Q
The constants of a differential ring are the elements
r
\partialr=0
\partial.
In the following identities,
\delta
R.
r\inR
c
R
\deltac=0
r\inR
u
R,
n
r\inR
u1,\ldots,un
R,
e1,\ldots,en
A derivation operator or higher-order derivation is the composition of several derivations. As the derivations of a differential ring are supposed to commute, the order of the derivations does not matter, and a derivation operator may be written aswhere
\delta1,\ldots,\deltan
e1,\ldots,en
The sum
o=e1+ … +en
o=1
o=0
A derivative of an element
x
x,
e1 | |
\delta | |
1 |
\circ … \circ
en | |
\delta | |
n |
(x).
I
R
R
\partial
x\inI.
The radical of a differential ideal is the same as its radical as an algebraic ideal, that is, the set of the ring elements that have a power in the ideal. The radical of a differential ideal is also a differential ideal. A radical or perfect differential ideal is a differential ideal that equals its radical. A prime differential ideal is a differential ideal that is prime in the usual sense; that is, if a product belongs to the ideal, at least one of the factors belongs to the ideal. A prime differential ideal is always a radical differential ideal.
A discovery of Ritt is that, although the classical theory of algebraic ideals does not work for differential ideals, a large part of it can be extended to radical differential ideals, and this makes them fundamental in differential algebra.
The intersection of any family of differential ideals is a differential ideal, and the intersection of any family of radical differential ideals is a radical differential ideal.It follows that, given a subset
S
The algebraic ideal generated by
S
S,
(S)
\langleS\rangle.
The differential ideal generated by
S
S
[S].
S
[S]
The radical differential ideal generated by
S
\{S\}.
A differential polynomial over a differential field
K
K,
So, let
K
K(X)=K(x1,\ldots,xn)
\partiali
\partialixi=1
\partialixj=0
i ≠ j
For defining the ring of differential polynomials over
K
Y=\{y1,\ldots,yn\}
\partial1,\ldots,\partialn,
\Deltayi,
\Delta
K\{Y\}
n=1,
K\{y\}=K\left[y,\partialy,\partial2y,\partial3y,\ldots\right].
Even when
n=1,
Firstly, a finite number of differential polynomials involves together a finite number of indeterminates. Its follows that every property of polynomials that involves a finite number of polynomials remains true for differential polynomials. In particular, greatest common divisors exist, and a ring of differential polynomials is a unique factorization domain.
The second fact is that, if the field
K
K
R
R\{y\}
This Noetherian property implies that, in a ring of differential polynomials, every radical differential ideal is finitely generated as a radical differential ideal; this means that there exists a finite set of differential polynomials such that is the smallest radical differential idesl containing . This allows representing a radical differential ideal by such a finite set of generators, and computing with these ideals. However, some usual computations of the algebraic case cannot be extended. In particular no algorithm is known for testing membership of an element in a radical differential ideal or the equality of two radical differential ideals.
Another consequence of the Noetherian property is that a radical differential ideal can be uniquely expressed as the intersection of a finite number of prime differential ideals, called essential prime components of the ideal.
Elimination methods are algorithms that preferentially eliminate a specified set of derivatives from a set of differential equations, commonly done to better understand and solve sets of differential equations.
Categories of elimination methods include characteristic set methods, differential Gröbner bases methods and resultant based methods.
Common operations used in elimination algorithms include 1) ranking derivatives, polynomials, and polynomial sets, 2) identifying a polynomial's leading derivative, initial and separant, 3) polynomial reduction, and 4) creating special polynomial sets.
The ranking of derivatives is a total order and an admisible order, defined as:
Each derivative has an integer tuple, and a monomial order ranks the derivative by ranking the derivative's integer tuple. The integer tuple identifies the differential indeterminate, the derivative's multi-index and may identify the derivative's order. Types of ranking include:
\forallyi,yj\inY, \forall\theta\mu,\theta\nu\in\Theta : \operatorname{ord}(\theta\mu)\ge\operatorname{ord}(\theta\nu) ⇒ \theta\muyi\ge\theta\nuyj
\forallyi,yj\inY, \forall\theta\mu,\theta\nu\in\Theta : yi\geyj ⇒ \theta\muyi\ge\theta\nuyj
In this example, the integer tuple identifies the differential indeterminate and derivative's multi-index, and lexicographic monomial order, , determines the derivative's rank.
e1 | |
η(\delta | |
1 |
\circ … \circ
en | |
\delta | |
n |
(yj))=(j,e1,\ldots,en)
η(\theta\muyj)\gelexη(\theta\nuyk) ⇒ \theta\muyj\ge\theta\nuyk.
This is the standard polynomial form:
p=ad ⋅
d+ | |
u | |
p |
ad-1 ⋅
d-1 | |
u | |
p |
+ … +a1 ⋅ up+a0
up
ad,\ldots,a0
\deg | |
up |
(p)=d
Ip=ad
d | |
u | |
p |
Sp=
\partialp | |
\partialup |
Separant set is
SA=\{Sp\midp\inA\}
IA=\{Ip\midp\inA\}
Partially reduced (partial normal form) polynomial with respect to polynomial indicates these polynomials are non-ground field elements, , and
q
up
Partially reduced polynomial with respect to polynomial becomes reduced (normal form) polynomial with respect to if the degree of in is less than the degree of in .
An autoreduced polynomial set has every polynomial reduced with respect to every other polynomial of the set. Every autoreduced set is finite. An autoreduced set is triangular meaning each polynomial element has a distinct leading derivative.
Ritt's reduction algorithm identifies integers and transforms a differential polynomial using pseudodivision to a lower or equally ranked remainder polynomial that is reduced with respect to the autoreduced polynomial set . The algorithm's first step partially reduces the input polynomial and the algorithm's second step fully reduces the polynomial. The formula for reduction is:
fred\equiv
\prod | |
Ak\inA |
| |||||
I | |||||
Ak |
⋅
| |||||
S | |||||
Ak |
⋅ f,\pmod{[A]}with
i | |
Ak |
,
s | |
Ak |
\inN.
Set is a differential chain if the rank of the leading derivatives is and is reduced with respect to
Autoreduced sets and each contain ranked polynomial elements. This procedure ranks two autoreduced sets by comparing pairs of identically indexed polynomials from both autoreduced sets.
A1< … <Am\inA
B1< … <Bn\inB
i,j,k\inN
rankA<rankB
k\le\operatorname{minimum}(m,n)
Ai=Bi
Ak<Bk
\operatorname{rank}A<\operatorname{rank}B
n<m
Ai=Bi
1\lei\len
\operatorname{rank}A=\operatorname{rank}B
n=m
Ai=Bi
1\lei\len
A characteristic set is the lowest ranked autoreduced subset among all the ideal's autoreduced subsets whose subset polynomial separants are non-members of the ideal .
The delta polynomial applies to polynomial pair whose leaders share a common derivative, . The least common derivative operator for the polynomial pair's leading derivatives is , and the delta polynomial is:
\operatorname{\Delta-poly}(p,q)=Sq ⋅
\thetapqp | |
\thetap |
-Sp ⋅
\thetapqq | |
\thetaq |
A coherent set is a polynomial set that reduces its delta polynomial pairs to zero.
A regular system contains a autoreduced and coherent set of differential equations and a inequation set with set reduced with respect to the equation set.
Regular differential ideal and regular algebraic ideal are saturation ideals that arise from a regular system. Lazard's lemma states that the regular differential and regular algebraic ideals are radical ideals.
The Rosenfeld–Gröbner algorithm decomposes the radical differential ideal as a finite intersection of regular radical differential ideals. These regular differential radical ideals, represented by characteristic sets, are not necessarily prime ideals and the representation is not necessarily minimal.
The membership problem is to determine if a differential polynomial is a member of an ideal generated from a set of differential polynomials . The Rosenfeld–Gröbner algorithm generates sets of Gröbner bases. The algorithm determines that a polynomial is a member of the ideal if and only if the partially reduced remainder polynomial is a member of the algebraic ideal generated by the Gröbner bases.
The Rosenfeld–Gröbner algorithm facilitates creating Taylor series expansions of solutions to the differential equations.
Example 1: is the differential meromorphic function field with a single standard derivation.
Example 2: is a differential field with a linear differential operator as the derivation, for any polynomial
p(y)
Define as shift operator for polynomial .
A shift-invariant operator commutes with the shift operator: .
The Pincherle derivative, a derivation of shift-invariant operator , is .
Ring of integers is
(Z.\delta)
\delta(m+1)=\delta(m)+\delta(1)=\delta(m) ⇒ \delta(m+1)=\delta(m)
\delta(1)=0 \wedge \delta(m+1)=\delta(m) ⇒ \forall m\inZ, \delta(m)=0
Field of rational numbers is
(Q.\delta)
\forallr\inQ, \exists a\inZ, b\inZ/\{0\}, r=
a | |
b |
\delta(r)=\delta\left(
a | |
b |
\right)=
\delta(a) ⋅ b-a ⋅ \delta(b) | |
b2 |
=0
Constants form the subring of constants .
Element simply generates differential ideal in the differential ring .
Any ring with identity is a algebra. Thus a differential ring is a algebra.
If ring is a subring of the center of unital ring , then is an algebra. Thus, a differential ring is an algebra over its differential subring. This is the natural structure of an algebra over its subring.
Ring has irreducible polynomials, (normal, squarefree) and (special, ideal generator).
Ring has derivatives and
The
leading derivatives, and initials are:.
Symbolic integration uses algorithms involving polynomials and their derivatives such as Hermite reduction, Czichowski algorithm, Lazard-Rioboo-Trager algorithm, Horowitz-Ostrogradsky algorithm, squarefree factorization and splitting factorization to special and normal polynomials.
Differential algebra can determine if a set of differential polynomial equations has a solution. A total order ranking may identify algebraic constraints. An elimination ranking may determine if one or a selected group of independent variables can express the differential equations. Using triangular decomposition and elimination order, it may be possible to solve the differential equations one differential indeterminate at a time in a step-wise method. Another approach is to create a class of differential equations with a known solution form; matching a differential equation to its class identifies the equation's solution. Methods are available to facilitate the numerical integration of a differential-algebraic system of equations.
In a study of non-linear dynamical systems with chaos, researchers used differential elimination to reduce differential equations to ordinary differential equations involving a single state variable. They were successful in most cases, and this facilitated developing approximate solutions, efficiently evaluating chaos, and constructing Lyapunov functions. Researchers have applied differential elimination to understanding cellular biology, compartmental biochemical models, parameter estimation and quasi-steady state approximation (QSSA) for biochemical reactions. Using differential Gröbner bases, researchers have investigated non-classical symmetry properties of non-linear differential equations. Other applications include control theory, model theory, and algebraic geometry. Differential algebra also applies to differential-difference equations.
A vector space is a collection of vector spaces with integer degree for . A direct sum can represent this graded vector space:
V\bullet=oplusmVm
A differential graded vector space or chain complex, is a graded vector space with a differential map or boundary map with
dm\circdm+1=0
A cochain complex is a graded vector space with a differential map or coboundary map with
dm+1\circdm=0
A differential graded algebra is a graded algebra with a linear derivation with
d\circd=0
\foralla,b\inA, d(a ⋅ b)=d(a) ⋅ b+(-1)|a| ⋅ a ⋅ d(b)
|a|
a
A Lie algebra is a finite-dimensional real or complex vector space with a bilinear bracket operator with Skew symmetry and the Jacobi identity property.
[X,Y]=-[Y,X]
[X,[Y,Z]]+[Y,[Z,X]]+[Z,[X,Y]]=0
X,Y,Z\inl{g}
The adjoint operator, is a derivation of the bracket because the adjoint's effect on the binary bracket operation is analogous to the derivation's effect on the binary product operation. This is the inner derivation determined by .
\operatorname{ad}X([Y,Z])=[\operatorname{ad}X(Y),Z]+[Y,\operatorname{ad}X(Z)]
The universal enveloping algebra of Lie algebra is a maximal associative algebra with identity, generated by Lie algebra elements and containing products defined by the bracket operation. Maximal means that a linear homomorphism maps the universal algebra to any other algebra that otherwise has these properties. The adjoint operator is a derivation following the Leibniz product rule.
U(l{g})
X ⋅ Y-Y ⋅ X=[X,Y]
\operatorname{ad}X(Y ⋅ Z)=\operatorname{ad}X(Y) ⋅ Z+Y ⋅ \operatorname{ad}X(Z)
X,Y,Z\inU(l{g})
The Weyl algebra is an algebra over a ring with a specific noncommutative product:
pi ⋅ qi-qi ⋅ pi=1, : i\in\{1,...,n\}
pi ⋅ qj-qj ⋅ pi=0ifi\nej, pi ⋅ pj-pj ⋅ pi=0, qi ⋅ qj-qj ⋅ qi=0
qj(yk)=yj ⋅ yk, qj(c)=c ⋅ yjwithc\inK, pj(yj)=1, pj(yk)=0ifj\nek, pj(c)=0withc\inK
The associative, possibly noncommutative ring has derivation .
The pseudo-differential operator ring is a left containing ring elements :
ai\inA, i,imin\inN, |imin|>0 : L=
n | |
\sum | |
i\geimin |
ai ⋅ \partiali
The derivative operator is .
The binomial coefficient is
l({i\atopk}r)
Pseudo-differential operator multiplication is:
n | |
\sum | |
i\geimin |
ai ⋅ \partiali ⋅
m | |
\sum | |
j\gejmin |
bi ⋅ \partialj=\sumi,j;kl({i\atopk}r) ⋅ ai ⋅
k(b | |
d | |
j) |
⋅ \partiali+j-k
The Ritt problem asks is there an algorithm that determines if one prime differential ideal contains a second prime differential ideal when characteristic sets identify both ideals.
The Kolchin catenary conjecture states given a dimensional irreducible differential algebraic variety and an arbitrary point , a long gap chain of irreducible differential algebraic subvarieties occurs from to V.
The Jacobi bound conjecture concerns the upper bound for the order of an differential variety's irreducible component. The polynomial's orders determine a Jacobi number, and the conjecture is the Jacobi number determines this bound.