Message Boards Message Boards

6 Replies
0 Total Likes
View groups...
Share this post:

Sorting noncommuting operators?

I'd like to be able to sort noncommuting operators, using their commutation relations. This shows up a lot in quantum mechanics. Any suggestions?

POSTED BY: Frank Kampas
6 Replies
In[1]:= Quiet @ Needs["HighEnergyPhysics`FeynCalc`"]

Loading FeynCalc from C:\Users\Frank\AppData\Roaming\Mathematica\Applications\HighEnergyPhysics

FeynCalc 8.2.0 For help, type ?FeynCalc, open FeynCalcRef8.nb or visit

Loading FeynArts, see for documentation

FeynArts 3.7 patched for use with FeynCalc

In[2]:= DeclareNonCommutative[a, ad]

In[3]:= Commutator[a, ad] = 1

Out[3]= 1

In[4]:= DotSimplify[ Commutator[a,]]

Out[4]= 2 ad

In[5]:= DotSimplify[ Commutator[a,]]

Out[5]= 3

In[6]:= DotSimplify[ Commutator[a,]]

Out[6]= 4
POSTED BY: Frank Kampas

Looks like FeynCalc is the way to go.

Here's an example problem from an online course I'm taking:

Recalling that [a^,a^†]=1, calculate the commutator:


a is an annihilation operator and a dagger is a creation operator.

POSTED BY: Frank Kampas

Also possibly of relevance is the section "Some noncommutative algebraic manipulation" in the notebook available here.

POSTED BY: Daniel Lichtblau
Mathematica 10.0 for Linux x86 (64-bit)
Copyright 1988-2014 Wolfram Research, Inc.

In[1]:= !!fordan
(* FeynCalc can be installed by
(* load FeynCalc *)

(* FI just sets the output to FeynCalcInternalForm, kind of FullForm *)

(* remember in FeynCalc:  .  is used as non-commutative operator *)

DeclareNonCommutative[a,b,c];  (* declares a,b,c to be non commutative *)

test = a . ( 2*b - z c).a ;

(* while now we add a commutator  [a,c] = 1 *)

Commutator[a, c] =1;

Expand @ DotSimplify @ test

In[2]:= <<fordan
Loading FeynCalc from /home/rolfm/.Mathematica/Applications/HighEnergyPhysics
FeynCalc 8.2.0 Type ?FeynCalc for help or visit
$PrePrint is set to FeynCalcForm. Use FI and FC to change the display format.
Loading FeynArts, see for documentation
FeynArts 3.7 patched for use with FeynCalc

Out[2]= -(a*z) + 2*a . b . a - z*c . a . a
POSTED BY: Rolf Mertig

Some simple examples of input and desired output would (as ever) be useful.

POSTED BY: Daniel Lichtblau

One possibility is to use DotSimplify from FeynCalc. Or you look at the source and use the relevant parts.

POSTED BY: Rolf Mertig
Reply to this discussion
Community posts can be styled and formatted using the Markdown syntax.
Reply Preview
or Discard

Group Abstract Group Abstract