Beispielanwendungen - Fortsetzung

Dirac - Gleichung

Die freien Felder, die die Feldgleichungen lösen, werden wiederum vorgegeben:

Ψ := p r = 1 2 ( Sqrt [ m / ( V p 0 ) ] c ^ p 1 , p 2 , p 3 r u r Exp [ exp ] ) + p r = 1 2 ( Sqrt [ m / ( V p 0 ) ] d ^ p 1 , p 2 , p 3 r , v r Exp [ + exp ] ) Ψ_ ¯ := ( Hermitian [ Ψ ] γ 0 ) // QASimplify [ ]

  1. Der Hamiltonoperator wird definiert und ausgewertet

    H := ( Ψ ¯ // QASimplify[] ) i = 1 3 ( γ i ( D x i ( Ψ // Evaluate ) /. p p 1 ) ) + m ( Ψ ¯ // QASimplify[] ) ( Ψ /. p p 1 ) H := V H 3 x H // QAEvaluate2 // Timing { 154.141 , 2 1 ^ p p 0 + p c ^ p 1 , p 2 , p 3 1 , c ^ p 1 , p 2 , p 3 1 p 0 + p c ^ p 1 , p 2 , p 3 2 , c ^ p 1 , p 2 , p 3 2 p 0 + p d ^ p 1 , p 2 , p 3 1 , d ^ p 1 , p 2 , p 3 1 p 0 + p d ^ p 1 , p 2 , p 3 2 , d ^ p 1 , p 2 , p 3 2 p 0 }

  2. Der 3er - Impuls wird definiert und ausgewertet

    P i_ := ( Ψ ¯ // QASimplify [ ] ) ( γ 0 ( D x i ( Ψ // Evaluate ) /. p p1 ) ) P i_ := V P i 3 x Map [ ( P i // QAEvaluate2 ) & , { 1 , 2 , 3 } ] // Timing { 116.688 , { p c ^ p 1 , p 2 , p 3 1 , c ^ p 1 , p 2 , p 3 1 p 1 p c ^ p 1 , p 2 , p 3 2 , c ^ p 1 , p 2 , p 3 2 p 1 p d ^ p 1 , p 2 , p 3 1 , d ^ p 1 , p 2 , p 3 1 p 1 p d ^ p 1 , p 2 , p 3 2 , d ^ p 1 , p 2 , p 3 2 p 1 , p c ^ p 1 , p 2 , p 3 1 , c ^ p 1 , p 2 , p 3 1 p 2 p c ^ p 1 , p 2 , p 3 2 , c ^ p 1 , p 2 , p 3 2 p 2 p d ^ p 1 , p 2 , p 3 1 , d ^ p 1 , p 2 , p 3 1 p 2 p d ^ p 1 , p 2 , p 3 2 , d ^ p 1 , p 2 , p 3 2 p 2 , p c ^ p 1 , p 2 , p 3 1 , c ^ p 1 , p 2 , p 3 1 p 3 p c ^ p 1 , p 2 , p 3 2 , c ^ p 1 , p 2 , p 3 2 p 3 p d ^ p 1 , p 2 , p 3 1 , d ^ p 1 , p 2 , p 3 1 p 3 p d ^ p 1 , p 2 , p 3 2 , d ^ p 1 , p 2 , p 3 2 p 3 } }

  3. Der 4er - Strom wird definiert und die Ladung ausgewertet (ohne multiplikative Ladungskonstante)

    J i : = ( Ψ QASimplify [ ] ) ( γ i ( Ψ /. p p 1 ) ) J i : = v J i 3 x QAEvaluate2 [ J 0 ] // Timing { 41.5313 , 2 1 ^ p 1 + p c ^ p 1 , p 2 , p 3 1 , c ^ p 1 , p 2 , p 3 1 + p c ^ p 1 , p 2 , p 3 2 , c ^ p 1 , p 2 , p 3 2 p d ^ p 1 , p 2 , p 3 1 , d ^ p 1 , p 2 , p 3 1 p d ^ p 1 , p 2 , p 3 2 , d ^ p 1 , p 2 , p 3 2 }