Wigner3j.m
6 Downloads
Updated 21 Jun 2004
No License
Wigner3j(j1,j2,j,m1,m2,m) returns the Wigner 3j-symbol, where j1, j2, j, m1, m2, and m are half-integers. Physically, the Wigner 3j-symbol is closely related to the Clebsch-Gordon coefficient <j1,j2,m1,m2|j1,j2,j,m>, the square of which is the probability that a system of two particles with angular momentum j1 and j2 respectively and z-component of angular momentum m1 and m2 respectively has total angular momentum j and z-component of total angular momentum m.
I haven't tested this function thoroughly yet, so there may be some bugs. Please let me know if you find any.
Cite As
David Terr (2021). Wigner3j.m (https://www.mathworks.com/matlabcentral/fileexchange/5275-wigner3j-m), MATLAB Central File Exchange. Retrieved .
Comments and Ratings (10)
MATLAB Release Compatibility
Platform Compatibility
Windows macOS LinuxCategories
Tags
Acknowledgements
Inspired: Wigner3j symbol
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!Discover Live Editor
Create scripts with code, output, and formatted text in a single executable document.
Select a Web Site
Choose a web site to get translated content where available and see local events and offers. Based on your location, we recommend that you select: .
Select web siteYou can also select a web site from the following list:
Americas
- América Latina (Español)
- Canada (English)
- United States (English)
Europe
- Belgium (English)
- Denmark (English)
- Deutschland (Deutsch)
- España (Español)
- Finland (English)
- France (Français)
- Ireland (English)
- Italia (Italiano)
- Luxembourg (English)
- Netherlands (English)
- Norway (English)
- Österreich (Deutsch)
- Portugal (English)
- Sweden (English)
- Switzerland
- United Kingdom (English)
The file is related to the author's ClebschGordan.m file. Adopt Kobi Kraus's suggestion and replace factorial(p) as gamma(p+1) will be much more efficient.
Still, the file does not check if the condition m1 + m2 = m3 is fulfilled. This should still be added!! In the current version it's very problematic to use the code! False results e.g. for angular momentum couplings are inevitable with this code.
Ops, did not read the table caption, it was 3jsymbol squared in table. now it match better for a few test cases
Made quick test against alpha=beta=gamma=0 3jsymbols in Table 2, in D.M. Brink and R.G. Satchler "Angular Momentum" Oxford Science Publications and can not get a match
It doesn't respect wigner 3j symetries
Wigner3j(1,1,1,1,-1,-1) is not 0
The file does not check wether the condition m1 + m2 = m3 is fulfilled. This should be added.
The loop in rows 59-64 can be replaced by a simple direct sum:
t = tmin : tmax;
wigner = sum( (-1).^t ./ ( factorial(t) .* factorial(t-t1) .* factorial(t-t2) ...
.* factorial(t3-t) .* factorial(t4-t) .* factorial(t5-t) ) );