13.07.2015 Views

Algebraic development of many-body perturbation theory in ...

Algebraic development of many-body perturbation theory in ...

Algebraic development of many-body perturbation theory in ...

SHOW MORE
SHOW LESS
  • No tags were found...

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

D Symbolic computations with NCoperators 901. Second Quantisation Representation (SQR)2. Angular Momentum Theory (AMT)3. Rayleigh–Schröd<strong>in</strong>ger Perturbation Theory (RSPT)4. Unstructured External Programm<strong>in</strong>g (UEP)The packages store Mathematica codes that are loaded <strong>in</strong>to a Mathematica session by us<strong>in</strong>gthe function «“NCoperators.m” or Get[“NCoperators.m”]. The NCoperatorsis designed so that each block can be brought <strong>in</strong>to action separately. This is easily done withGet[“PackageName.m”]. All these functions (Get[]) are located <strong>in</strong> a s<strong>in</strong>gle packageNCoperators.m.D.1 SQR and AMT blocksThe SQR and AMT blocks are closely related and therefore they will be considered together.The present blocks conta<strong>in</strong> a huge number <strong>of</strong> functions, each <strong>of</strong> them be<strong>in</strong>g useful for a specialcase <strong>of</strong> <strong>in</strong>terest. In the present overview, only a few <strong>of</strong> them will be discussed.In theoretical atomic physics, a typical task is to f<strong>in</strong>d the angular coefficients that relate twodist<strong>in</strong>ct ( momenta coupl<strong>in</strong>g schemes. Consider, for example, the momenta recoupl<strong>in</strong>g coefficientj1 j 2 (j 12 )j 3 j ∣ ∣j 2 j 3 (j 23 )j 1 j ) . The expression is known, and it is{ }(−1) 2j+j 2+j 3 −j 23[j 12 , j 23 ] 1/2 j1 j 2 j 12.j 3 j j 23Mak<strong>in</strong>g use <strong>of</strong> NCoperators, the algorithm to obta<strong>in</strong> the latter formula is displayed <strong>in</strong> Fig. 1.Fig. 1: A computation <strong>of</strong> recoupl<strong>in</strong>g coefficient ( j 1 j 2 (j 12 )j 3 j ∣ ∣ j2 j 3 (j 23 )j 1 j ) with NCoperatorsOut[25] conta<strong>in</strong>s two Clebsch–Gordan coefficients labelled by CGcoeff[] – the standardMathematica function ClebschGordan[] is <strong>in</strong>sufficient for symbolic manipulations.The function Recoupl<strong>in</strong>g[] <strong>in</strong>itiates the momenta recoupl<strong>in</strong>g followed by Refs. [10–12].SymbolSum[] plays a role <strong>of</strong> a standard Mathematica function Sum[] adapted to symboliccomputations. The function Make[] <strong>in</strong>itiates the momenta recoupl<strong>in</strong>g with<strong>in</strong> a given Clebsch–Gordan coefficient. The simplification <strong>of</strong> considered sum is performed by SumSimplify.Out[31] is pr<strong>in</strong>ted <strong>in</strong> a standard Unix output form specified by the font family «cmmi10». Obta<strong>in</strong>edformula is easy to convert to a L A TEX language by mak<strong>in</strong>g use <strong>of</strong> Mathematica function

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!