b b b b b b b b b Geometric Singularities of Algebraic - - PowerPoint PPT Presentation

b
SMART_READER_LITE
LIVE PREVIEW

b b b b b b b b b Geometric Singularities of Algebraic - - PowerPoint PPT Presentation

b b b b b b b b b b b b b b b b b b b Geometric Singularities of Algebraic Differential Equations Werner M. Seiler ur Mathematik, Universit Institut f at Kassel (joint work with Markus Lange-Hegermann , RWTH Aachen) W.M.


slide-1
SLIDE 1

b

b

b

b

b

b

b

b

b

b

b

b

b

b

b b

b

b

b

W.M. Seiler: Geometric Singularities of Algebraic Differential Equations – 1

Geometric Singularities of Algebraic Differential Equations

Werner M. Seiler Institut f¨ ur Mathematik, Universit¨ at Kassel (joint work with Markus Lange-Hegermann, RWTH Aachen)

slide-2
SLIDE 2

b

b

b

b

b b

b

b

b

b

Introduction

Introduction Algebraic Differential Equations Vessiot Distribution and Generalised Solutions Regular Differential Equations Geometric Singularities Thomas Decomposition Detection of Singularities

W.M. Seiler: Geometric Singularities of Algebraic Differential Equations – 2

singularities of differential equations

=

singularities of solutions of differential equations

  • related, but different topics
  • no discussion of shocks, blow-ups, etc.
slide-3
SLIDE 3

b

b

b

b

b b

b

b

b

b

Introduction

Introduction Algebraic Differential Equations Vessiot Distribution and Generalised Solutions Regular Differential Equations Geometric Singularities Thomas Decomposition Detection of Singularities

W.M. Seiler: Geometric Singularities of Algebraic Differential Equations – 2

“Interpolation” between three domains:

  • differential algebra
  • differential topology
  • differential algebraic equations

together with techniques from (differential) algebraic geometry Current goal: detect all singularities of given system of (ordinary or partial) differential equations

slide-4
SLIDE 4

b

b

b

b

b b

b

b

b

b

Introduction

Introduction Algebraic Differential Equations Vessiot Distribution and Generalised Solutions Regular Differential Equations Geometric Singularities Thomas Decomposition Detection of Singularities

W.M. Seiler: Geometric Singularities of Algebraic Differential Equations – 2

Differential Algebra

  • here mainly differential ideal theory
  • covers automatically systems and all orders
  • founded by Ritt in early 20th century
  • central goal: understanding singular integrals
  • ldest example:

Taylor (1715)

  • best known example:

Clairaut equation (1734)

u = xu′ + f(u′)

mit f ′′(z) = 0 ∀z general integral:

u(x) = cx + f(c)

singular integral:

x(τ) = −f ′(τ), u(τ) = −τf ′(τ) + f(τ)

slide-5
SLIDE 5

b

b

b

b

b b

b

b

b

b

Introduction

Introduction Algebraic Differential Equations Vessiot Distribution and Generalised Solutions Regular Differential Equations Geometric Singularities Thomas Decomposition Detection of Singularities

W.M. Seiler: Geometric Singularities of Algebraic Differential Equations – 2

Differential Algebra

f(z) = −1

4z2

slide-6
SLIDE 6

b

b

b

b

b b

b

b

b

b

Introduction

Introduction Algebraic Differential Equations Vessiot Distribution and Generalised Solutions Regular Differential Equations Geometric Singularities Thomas Decomposition Detection of Singularities

W.M. Seiler: Geometric Singularities of Algebraic Differential Equations – 2

Differential Topology

  • singularities of smooth maps between manifolds
  • submanifolds of jet bundles provide geometric model for differential

equations

  • natural projections between jet bundles of different order
  • critical points = geometric singularities
  • distinction regular and irregular singularities
  • complete classifications of singularities of scalar ordinary differential

equations of first or second order

  • hardly any works on (general) systems
slide-7
SLIDE 7

b

b

b

b

b b

b

b

b

b

Introduction

Introduction Algebraic Differential Equations Vessiot Distribution and Generalised Solutions Regular Differential Equations Geometric Singularities Thomas Decomposition Detection of Singularities

W.M. Seiler: Geometric Singularities of Algebraic Differential Equations – 2

Differential Topology

(u′)2 + u2 + x2 = 1

slide-8
SLIDE 8

b

b

b

b

b b

b

b

b

b

Introduction

Introduction Algebraic Differential Equations Vessiot Distribution and Generalised Solutions Regular Differential Equations Geometric Singularities Thomas Decomposition Detection of Singularities

W.M. Seiler: Geometric Singularities of Algebraic Differential Equations – 2

Differential Algebraic Equations

  • mainly analytic theory of quasi-linear systems A(x, u)u′ = F(x, u)

(including very large systems!)

  • A not necessarily of maximal rank and rank may jump
  • impasse points already discussed in 1960s by electrical engineers
  • lead to jump phenomena in solutions
  • n one side interpreted as sign of bad model. . .
  • . . . on the other side experiments often show similar behaviour
slide-9
SLIDE 9

b

b

b

b

b

b

b

b

b

b

b

b

b

b

b

b

b

b

b

Algebraic Differential Equations

Introduction Algebraic Differential Equations Vessiot Distribution and Generalised Solutions Regular Differential Equations Geometric Singularities Thomas Decomposition Detection of Singularities

W.M. Seiler: Geometric Singularities of Algebraic Differential Equations – 3

  • consider holomorphic function f : U ⊆
Cn → Cm, u = f(z)

from now on: n, m fixed, U ignored

  • q-jet [f](q)

z

  • equivalence class of all holomorphic functions

g :

Cn → Cm with same Taylor polynomial of degree q around

z ∈ U as f

  • jet bundle Jq = Jq( Cn,
Cm)
  • set of all q-jets [f](q)

z

  • manifold of dimension dq = n + m

n+q

q

  • (may be identified with
Cdq)
  • local coordinates (z, u(q)) corresponding to expansion point z

and derivatives up to order q

  • natural projections for 0 ≤ r < q

πq

r :

  • Jq

− → Jr [f](q)

z

− → [f](r)

z

πq :

  • Jq

− →

Cn

[f](q)

z

− → z

slide-10
SLIDE 10

b

b

b

b

b

b

b

b

b

b

b

b

b

b

b

b

b

b

b

Algebraic Differential Equations

Introduction Algebraic Differential Equations Vessiot Distribution and Generalised Solutions Regular Differential Equations Geometric Singularities Thomas Decomposition Detection of Singularities

W.M. Seiler: Geometric Singularities of Algebraic Differential Equations – 3

Definition:

  • algebraic jet set of order q
  • locally Zariski closed set Rq ⊆ Jq

(i.e.: difference of two varieties)

  • algebraic differential equation of order q
  • algebraic jet set Rq ⊆ Jq such that restricted projection πq|Rq

dominant (i.e.: image Zariski dense in

Cn)

both generalisation and restriction of classical geometric definition:

  • nly polynomial non-linearities admitted
  • Rq may have algebraic singularities
  • equations and inequations admitted
  • dominance replaces surjectivity, permits “special points” in
Cn
  • πq|Rq not necessarily submersive
slide-11
SLIDE 11

b

b

b

b

b

b

b

b

b

b

b

b

b

b

b

b

b

b

b

Algebraic Differential Equations

Introduction Algebraic Differential Equations Vessiot Distribution and Generalised Solutions Regular Differential Equations Geometric Singularities Thomas Decomposition Detection of Singularities

W.M. Seiler: Geometric Singularities of Algebraic Differential Equations – 3

holomorphic function f defines section

σf :

Cn → Cn × Cm = J0, z →
  • z, f(z)
  • = [f](0)

z

(graph of f is image of σf ) consider prolonged section

jqσf :

Cn → Jq, z → [f](q)

z

Def:

f (resp. σf ) (classical) solution of differential equation Rq ⊆ Jq

  • im (jqσf) ⊆ Rq
slide-12
SLIDE 12

b

b

b

b

b

b

b

b

b

b

b

b Vessiot Distribution and Generalised Solutions

Introduction Algebraic Differential Equations Vessiot Distribution and Generalised Solutions Regular Differential Equations Geometric Singularities Thomas Decomposition Detection of Singularities

W.M. Seiler: Geometric Singularities of Algebraic Differential Equations – 4

What distinguishes Jq from

Cdq?
  • contact structure on Jq

Def: contact distribution Cq ⊂ TJq generated by vector fields

C(q)

i

= ∂zi +

  • α
  • 0≤|µ|<q

µ+1i∂uα

µ

1 ≤ i ≤ n Cµ

α = ∂uα

µ

1 ≤ α ≤ m, |µ| = q

Prop: section γ :

Cn → Jq of the form γ = jqσf for function f

⇐ ⇒ Tim(γ) ⊂ Cq

Proof: chain rule!

slide-13
SLIDE 13

b

b

b

b

b

b

b

b

b

b

b

b Vessiot Distribution and Generalised Solutions

Introduction Algebraic Differential Equations Vessiot Distribution and Generalised Solutions Regular Differential Equations Geometric Singularities Thomas Decomposition Detection of Singularities

W.M. Seiler: Geometric Singularities of Algebraic Differential Equations – 4

consider prolonged solution jqσf of equation Rq ⊆ Jq:

  • integral elements
  • im(jqσf)

ur ρ ∈ im(jqσf)

  • solution of Rq

= ⇒ Tρ

  • im(jqσf)
  • ⊆ TρRq
  • prolonged section

= ⇒ Tρ

  • im(jqσf)
  • ⊆ Cq|ρ
slide-14
SLIDE 14

b

b

b

b

b

b

b

b

b

b

b

b Vessiot Distribution and Generalised Solutions

Introduction Algebraic Differential Equations Vessiot Distribution and Generalised Solutions Regular Differential Equations Geometric Singularities Thomas Decomposition Detection of Singularities

W.M. Seiler: Geometric Singularities of Algebraic Differential Equations – 4

consider prolonged solution jqσf of equation Rq ⊆ Jq:

  • integral elements
  • im(jqσf)

ur ρ ∈ im(jqσf)

  • solution of Rq

= ⇒ Tρ

  • im(jqσf)
  • ⊆ TρRq
  • prolonged section

= ⇒ Tρ

  • im(jqσf)
  • ⊆ Cq|ρ

Def: Vessiot space in point ρ on algebraic jet set Rq

Vρ[Rq] = TρRq ∩ Cq|ρ

  • dim Vρ[Rq] generally depends on ρ
  • regular distribution only on Zariski open subset of Rq
  • computing Vessiot distribution V[Rq] requires only linear algebra
slide-15
SLIDE 15

b

b

b

b

b

b

b

b

b

b

b

b Vessiot Distribution and Generalised Solutions

Introduction Algebraic Differential Equations Vessiot Distribution and Generalised Solutions Regular Differential Equations Geometric Singularities Thomas Decomposition Detection of Singularities

W.M. Seiler: Geometric Singularities of Algebraic Differential Equations – 4

consider prolonged solution jqσf of equation Rq ⊆ Jq:

  • integral elements
  • im(jqσf)

ur ρ ∈ im(jqσf)

  • solution of Rq

= ⇒ Tρ

  • im(jqσf)
  • ⊆ TρRq
  • prolonged section

= ⇒ Tρ

  • im(jqσf)
  • ⊆ Cq|ρ

Def: Vessiot space in point ρ on algebraic jet set Rq

Vρ[Rq] = TρRq ∩ Cq|ρ

  • (geometric) symbol:

Nq,ρ = TρRq ∩ Vρπq

q−1 ⊆ V[Rq]

  • decompose V[Rq] = Nq ⊕ H with complement H (non-unique)
  • if dim H = n
  • H horizontal space of Vessiot connection
  • if Vessiot connection flat (i. e. differential equation integrable)
  • integral manifolds images of prolonged sections
slide-16
SLIDE 16

b

b

b

b

b

b

b

b

b

b

b

b Vessiot Distribution and Generalised Solutions

Introduction Algebraic Differential Equations Vessiot Distribution and Generalised Solutions Regular Differential Equations Geometric Singularities Thomas Decomposition Detection of Singularities

W.M. Seiler: Geometric Singularities of Algebraic Differential Equations – 4

Def: differential equation Rq ⊆ Jq

  • generalised solution
  • n-dimensional integral manifold N ⊆ Rq
  • f Vessiot distribution V[Rq]
  • geometric solution
  • projection πq

0(N) of generalised solution

N ⊆ Rq

comparison with classical solutions:

  • geometric solution not necessarily graph of function f
  • geometric solution πq

0(N) graph of classical solution

⇐ ⇒ N everywhere transversal to πq

  • geometric solution allow for modelling of multivalued solutions

(“breaking waves”)

slide-17
SLIDE 17

b

b

b

b

b

b

b

b

b

b

b

b Vessiot Distribution and Generalised Solutions

Introduction Algebraic Differential Equations Vessiot Distribution and Generalised Solutions Regular Differential Equations Geometric Singularities Thomas Decomposition Detection of Singularities

W.M. Seiler: Geometric Singularities of Algebraic Differential Equations – 4

Vessiot distribution and generalised solutions for sphere example

(u′)2 + u2 + x2 = 1

slide-18
SLIDE 18

b

b

b

b

b

b

b

b

b

b

b

b

b

b

b b

b

b

b

b

Regular Differential Equations

Introduction Algebraic Differential Equations Vessiot Distribution and Generalised Solutions Regular Differential Equations Geometric Singularities Thomas Decomposition Detection of Singularities

W.M. Seiler: Geometric Singularities of Algebraic Differential Equations – 5

P =

C[x1, . . . , xn] polynomial ring in n variables, point ρ ∈ Jq

Φ : Jq →

C holomorphic function; ¯

q maximal order of jet variable uα

µ

actually appearing in Φ

  • principal part of Φ in ρ

ppρ Φ =

m

  • α=1
  • |µ|=¯

q

∂Φ ∂uα

µ

(ρ)xµeα ∈ Pm Rq described by equations Φ1 = 0, . . . , Φr = 0 (inequations irrelevant)

  • in every point ρ ∈ Rq polynomial module

M[ρ] = ppρ Φ1, . . . , ppρ Φr

Def: Hilbert function of Rq in ρ

  • Hilbert function of factor module Pm/M[ρ]
slide-19
SLIDE 19

b

b

b

b

b

b

b

b

b

b

b

b

b

b

b b

b

b

b

b

Regular Differential Equations

Introduction Algebraic Differential Equations Vessiot Distribution and Generalised Solutions Regular Differential Equations Geometric Singularities Thomas Decomposition Detection of Singularities

W.M. Seiler: Geometric Singularities of Algebraic Differential Equations – 5

Def: algebraic differential equation Rq ⊆ Jq regular

  • Rq smooth (i.e. manifold)
  • Hilbert functions independent of point ρ ∈ Rq

idea: uniform behaviour of all “characteristic values” over Rq, in particular

  • dim Vρ[Rq]
  • dim Nq,ρ
  • size of formal solution space in ρ
slide-20
SLIDE 20

b

b

b

b

b

b

b

Geometric Singularities

Introduction Algebraic Differential Equations Vessiot Distribution and Generalised Solutions Regular Differential Equations Geometric Singularities Thomas Decomposition Detection of Singularities

W.M. Seiler: Geometric Singularities of Algebraic Differential Equations – 6

consider algebraic differential equation Rq ⊆ Jq

  • algebraic singularities
  • singularities in the sense of algebraic

geometry

  • ignored in sequel (not much known)
  • determination classical problem in algebraic geometry

(Jacobi criterion)

  • geometric singularities
  • critical points of restricted projection

ˆ πq

q−1 : Rq −

→ πq

q−1(Rq)

(i.e. Tρˆ

πq

q−1 not surjective)

  • points where dimension of symbol jumps
slide-21
SLIDE 21

b

b

b

b

b

b

b

Geometric Singularities

Introduction Algebraic Differential Equations Vessiot Distribution and Generalised Solutions Regular Differential Equations Geometric Singularities Thomas Decomposition Detection of Singularities

W.M. Seiler: Geometric Singularities of Algebraic Differential Equations – 6

let Rq ⊆ Jq be union of algebraic jet sets; smooth point ρ ∈ Rq is

  • regular
  • ρ has open neighbourhood where V[Rq] regular and

V[Rq] = Nq ⊕ H with dim H = n

  • regular singular
  • ρ has open neighbourhood where V[Rq]

regular, but dim Hρ < n

  • irregular singular
  • V[Rq] not regular on any neighbourhood of ρ
  • purely irregular singular
  • irregular singularity with dim Hρ = n

difference to classical definitions:

  • partial differential equations require consideration of neighbourhood of

point

  • “right” dimension of V[Rq] a priori not known
  • for ordinary differential equations pointwise analysis sufficient
slide-22
SLIDE 22

b

b

b

b

b

b

b

Geometric Singularities

Introduction Algebraic Differential Equations Vessiot Distribution and Generalised Solutions Regular Differential Equations Geometric Singularities Thomas Decomposition Detection of Singularities

W.M. Seiler: Geometric Singularities of Algebraic Differential Equations – 6

formally integrable ODE:

Rq ⊆ Jq

local description:

Φ(z, u(q)) = 0 Rq of finite type

  • almost everywhere dim Vρ[Rq] = 1

Prop: point ρ ∈ Rq

  • regular

⇐ ⇒ rank

  • C(q)Φ
  • ρ = m
  • regular singular

⇐ ⇒ ρ not regular and rank

  • C(q)Φ | C(q)

transΦ

  • ρ = m
slide-23
SLIDE 23

b

b

b

b

b

b

b

Geometric Singularities

Introduction Algebraic Differential Equations Vessiot Distribution and Generalised Solutions Regular Differential Equations Geometric Singularities Thomas Decomposition Detection of Singularities

W.M. Seiler: Geometric Singularities of Algebraic Differential Equations – 6

formally integrable ODE:

Rq ⊆ Jq

local description:

Φ(z, u(q)) = 0 Rq of finite type

  • almost everywhere dim Vρ[Rq] = 1

Thm: assume Rq without irregular singularities

  • ρ ∈ Rq regular point

= ⇒

(i) unique classical solution f exists with ρ ∈ im jqσf (ii) solution f may be continued in any direction until jqσf reaches either boundary of Rq or a regular singularity

  • ρ ∈ Rq regular singularity

= ⇒

dichotomy (i) either two classical solutions f1, f2 exist with ρ ∈ im jqσfi (either both start or both end in ρ) (ii)

  • r one classical solution f exists with ρ ∈ im jqσf whose

derivative of order q + 1 in z = πq(ρ) is not defined

slide-24
SLIDE 24

b

b

b

b

b

b

b

Geometric Singularities

Introduction Algebraic Differential Equations Vessiot Distribution and Generalised Solutions Regular Differential Equations Geometric Singularities Thomas Decomposition Detection of Singularities

W.M. Seiler: Geometric Singularities of Algebraic Differential Equations – 6

formally integrable ODE:

Rq ⊆ Jq

local description:

Φ(z, u(q)) = 0 Rq of finite type

  • almost everywhere dim Vρ[Rq] = 1

Proof:

V[Rq] locally generated by vector field X ρ regular singularity = ⇒ X vertical wrt πq

dichotomy

  • does ∂z-component of X change sign in ρ?
slide-25
SLIDE 25

b

b

b

b

b

b

b

Geometric Singularities

Introduction Algebraic Differential Equations Vessiot Distribution and Generalised Solutions Regular Differential Equations Geometric Singularities Thomas Decomposition Detection of Singularities

W.M. Seiler: Geometric Singularities of Algebraic Differential Equations – 6

formally integrable ODE:

Rq ⊆ Jq

local description:

Φ(z, u(q)) = 0 Rq of finite type

  • almost everywhere dim Vρ[Rq] = 1

let ρ ∈ Rq be an irregular singularity

  • consider simply connected domain U ⊂ Rq without irregular

singularities such that ρ ∈ U

  • in U Vessiot distribution V[Rq] generated by vector field X

Thm: generically every continuation of X to ρ vanishes Consequence: solution behaviour in neighbourhood of isolated irregular singularity ρ analysable with dynamical systems theory (mainly determined by eigenstructure of Jacρ X)

slide-26
SLIDE 26

b

b

b

b

b

b

b

Geometric Singularities

Introduction Algebraic Differential Equations Vessiot Distribution and Generalised Solutions Regular Differential Equations Geometric Singularities Thomas Decomposition Detection of Singularities

W.M. Seiler: Geometric Singularities of Algebraic Differential Equations – 6

Example:

(u′)3 + uu′ − x = 0

(hyperbolic gather) singularity curve (criminant):

3(u′)2 + u = 0

slide-27
SLIDE 27

b

b

b

b

b

b

b

Geometric Singularities

Introduction Algebraic Differential Equations Vessiot Distribution and Generalised Solutions Regular Differential Equations Geometric Singularities Thomas Decomposition Detection of Singularities

W.M. Seiler: Geometric Singularities of Algebraic Differential Equations – 6

Example:

(u′)3 + uu′ − x = 0

(hyperbolic gather) second derivative of solution touching “tip” of discriminant does not exist solutions “change direction” when crossing discriminant

slide-28
SLIDE 28

b

b

b

b

b

b

b

Geometric Singularities

Introduction Algebraic Differential Equations Vessiot Distribution and Generalised Solutions Regular Differential Equations Geometric Singularities Thomas Decomposition Detection of Singularities

W.M. Seiler: Geometric Singularities of Algebraic Differential Equations – 6

Example:

(u′)3 + uu′ − x = 0

(hyperbolic gather) neighbourhood of irregular singularity the two solutions tangential to eigenvectors of Jac X intersect

slide-29
SLIDE 29

b

b

b

b

b

b

b

b

b

b

b

Thomas Decomposition

Introduction Algebraic Differential Equations Vessiot Distribution and Generalised Solutions Regular Differential Equations Geometric Singularities Thomas Decomposition Detection of Singularities

W.M. Seiler: Geometric Singularities of Algebraic Differential Equations – 7

Algebraic Case polynomial ring

C[z1, . . . , zn] with total ordering on variables
  • leader ld p
  • largest variable in polynomial p
  • consider p as univariate polynomial in ld p
  • initial init p
  • leading coefficient of p
  • separant sep p
  • ∂p/∂(ld p)

algebraic system

  • finite set of polynomial equations and inequations

S =

  • p1 = 0, . . . , ps = 0, q1 = 0, . . . , qt = 0
  • solution set

(locally closed wrt Zariski topology)

Sol S =

  • z ∈
Cn | pi(z) = 0, qj(z) = 0
slide-30
SLIDE 30

b

b

b

b

b

b

b

b

b

b

b

Thomas Decomposition

Introduction Algebraic Differential Equations Vessiot Distribution and Generalised Solutions Regular Differential Equations Geometric Singularities Thomas Decomposition Detection of Singularities

W.M. Seiler: Geometric Singularities of Algebraic Differential Equations – 7

Algebraic Case Def: simple algebraic system

  • triangular:
  • {ld pi, ld qj} \ {1}
  • = s + t
  • non-vanishing initials:

no equation init pi = 0 or init qj = 0 has solution in Sol S

  • square-free:

dito for separants Def: Thomas decomposition of algebraic system S

  • finitely many

simple systems S1, . . . , Sk such that Sol S disjoint union of all Sol Si

  • exists always
  • depends on ordering of variables
  • decomposes according to fibre cardinality for coordinate projections
  • can be determined algorithmically

(subresultants, case distinctions

  • expensive!)
slide-31
SLIDE 31

b

b

b

b

b

b

b

b

b

b

b

Thomas Decomposition

Introduction Algebraic Differential Equations Vessiot Distribution and Generalised Solutions Regular Differential Equations Geometric Singularities Thomas Decomposition Detection of Singularities

W.M. Seiler: Geometric Singularities of Algebraic Differential Equations – 7

Algebraic Case consider V

  • y3 + (3x + 1)y2 + (3x2 + 2x)y + x3

x y

Thomas decomposition

  • S1 =
  • y3 +(3x+1)y2 +(3x2 +2x)y +x3 = 0, 27x3 −4x = 0
  • S2 =
  • 6y2 −(27x2 −12x−6)y −3x2 +2x = 0, 27x3 −4x = 0
slide-32
SLIDE 32

b

b

b

b

b

b

b

b

b

b

b

Thomas Decomposition

Introduction Algebraic Differential Equations Vessiot Distribution and Generalised Solutions Regular Differential Equations Geometric Singularities Thomas Decomposition Detection of Singularities

W.M. Seiler: Geometric Singularities of Algebraic Differential Equations – 7

differential case ring of differential polynomials

  • K =
C(z1, . . . , zn)

rational functions

  • derivations

δi = ∂/∂zi

  • differential unknowns:

U = {u1, . . . , um}

  • jet variables uα

µ = δµuα

  • K{U} =
K

µ | 1 ≤ α ≤ m, µ ∈

Nn
  • (polynomial ring in infinitely many variables)

derivations can be extended:

δiuα

µ = uα µ+1i

  • distinguish:
  • algebraic ideal:

p1, . . . , ps

  • differential ideal:

p1, . . . , ps∆

  • set D =
C
  • zi, uα

µ

K{U}
  • Dq =
C
  • zi, uα

µ | |µ| ≤ q

  • coordinate ring of Jq
slide-33
SLIDE 33

b

b

b

b

b

b

b

b

b

b

b

Thomas Decomposition

Introduction Algebraic Differential Equations Vessiot Distribution and Generalised Solutions Regular Differential Equations Geometric Singularities Thomas Decomposition Detection of Singularities

W.M. Seiler: Geometric Singularities of Algebraic Differential Equations – 7

differential case ranking on

K{U}
  • total ordering ≺ of jet variables
  • uα ≺ δiuα

µ ≺ uβ ν

= ⇒ δiuα

µ ≺ δiuβ ν

extend concepts like leader, initial or separant differential system

  • finite set of differential polynomial equations and

inequations

S =

  • p1 = 0, . . . , ps = 0, q1 = 0, . . . , qt = 0
  • solution set
  • consider formal solutions

(different function spaces possible)

slide-34
SLIDE 34

b

b

b

b

b

b

b

b

b

b

b

Thomas Decomposition

Introduction Algebraic Differential Equations Vessiot Distribution and Generalised Solutions Regular Differential Equations Geometric Singularities Thomas Decomposition Detection of Singularities

W.M. Seiler: Geometric Singularities of Algebraic Differential Equations – 7

differential case Def: simple differential system

  • simple as algebraic system in the finitely many occuring jet variables
  • involutive for Janet division
  • no leader of inequation derivative of leader of equation

Def: Thomas decomposition of differential system S

  • finitely many

simple systems S1, . . . , Sk such that Sol S disjoint union of all Sol Si

  • exists always
  • algorithmically computable via combination of algebraic Thomas

decomposition and Janet-Riquier theory

slide-35
SLIDE 35

b

b

b

b

b

b

b

b

b

Detection of Singularities

Introduction Algebraic Differential Equations Vessiot Distribution and Generalised Solutions Regular Differential Equations Geometric Singularities Thomas Decomposition Detection of Singularities

W.M. Seiler: Geometric Singularities of Algebraic Differential Equations – 8

starting point: differential system S goal: all geometric singularities in given order q differential computation:

  • differential Thomas decomposition (other methods also possible)
  • simple differential systems Si
  • ne Si corresponds to general integral

all others yield singular integrals

  • all other kind of singularities eliminated
slide-36
SLIDE 36

b

b

b

b

b

b

b

b

b

Detection of Singularities

Introduction Algebraic Differential Equations Vessiot Distribution and Generalised Solutions Regular Differential Equations Geometric Singularities Thomas Decomposition Detection of Singularities

W.M. Seiler: Geometric Singularities of Algebraic Differential Equations – 8

starting point: differential system S goal: all geometric singularities in given order q algebraic analysis of simple differential system S

  • introduce suitable ideals
  • I(S) = p1, . . . , ps∆ : h∞ ⊆ D with

h =

i sep (pi) init (pi)

  • Iq(S) = I(S) ∩ Dq,

Kq = qj | ord(qj) ≤ qDq

  • algebraic jet set

Rq = Sol

  • Iq(S)
  • \ Sol
  • Kq(S)
  • ⊆ Jq
slide-37
SLIDE 37

b

b

b

b

b

b

b

b

b

Detection of Singularities

Introduction Algebraic Differential Equations Vessiot Distribution and Generalised Solutions Regular Differential Equations Geometric Singularities Thomas Decomposition Detection of Singularities

W.M. Seiler: Geometric Singularities of Algebraic Differential Equations – 8

starting point: differential system S goal: all geometric singularities in given order q algebraic analysis of simple differential system S

  • ansatz for Vessiot distribution of Rq

X =

  • i

aiC(q)

i

+

  • α,µ

µCµ α

extended polynomial ring DV

q = Dq[a, b] with b ≻ a ≻ u ≻ z

  • compute algebraic Thomas decomposition of system over DV

q

consisting of generators of Iq(S) plus equations for Vessiot distribution (linear in a, b)

  • solve “parametric linear system”
  • simple systems SV

i

and Si = SV

i ∩ Dq

slide-38
SLIDE 38

b

b

b

b

b

b

b

b

b

Detection of Singularities

Introduction Algebraic Differential Equations Vessiot Distribution and Generalised Solutions Regular Differential Equations Geometric Singularities Thomas Decomposition Detection of Singularities

W.M. Seiler: Geometric Singularities of Algebraic Differential Equations – 8

starting point: differential system S goal: all geometric singularities in given order q Def: regularity decomposition in order q of simple differential system

  • write Sol Iq(S) ⊆ Jq as disjoint union of finitely many regular

algebraic jet sets R(i)

q

⊂ Jq (components in order q)

  • singular closure R(i)

q

  • f component R(i)

q

  • union with all

components R(j)

q

lying in Zariski closure of R(i)

q

  • constituent in order q
  • algebraically simple system S′ such that

Sol (S′) ⊆ Sol

  • Iq(S)
  • and set of leaders of equations in S′ equal to

set of jet variables in ld p1, . . . , ld ps∆ ∩ Dq

slide-39
SLIDE 39

b

b

b

b

b

b

b

b

b

Detection of Singularities

Introduction Algebraic Differential Equations Vessiot Distribution and Generalised Solutions Regular Differential Equations Geometric Singularities Thomas Decomposition Detection of Singularities

W.M. Seiler: Geometric Singularities of Algebraic Differential Equations – 8

starting point: differential system S goal: all geometric singularities in given order q final analysis:

  • Thm:

algorithm yields regularity decomposition with all constituents regular algebraic differential equations

  • Prop:

union of solution sets of constituents Zariski dense in

Sol

  • Iq(S)
  • consider for each constituent singular closure

(component may lie in closure of several constituents!)

  • taxonomy of singularities via leaders of systems SV

i

(do variables a appear as leader?) and comparison with constituent