Difference between revisions of "Cg fast.m"
(Created page with "Fast Clebsch-Gordan coefficient function with limited accuracy. Use clebsch_gordan.m for machine-precision results. Inputs: indices - CG coefficient indices, or...") |
|||
| Line 1: | Line 1: | ||
| − | + | {{DISPLAYTITLE:cg_fast.m}} __NOTOC__ | |
| − | + | Clebsch-Gordan coefficient: the coefficient in front of Y(L,M) spherical harmonic in the expansion of the product of Y(L1,M1) and Y(L2,M2) spherical harmonics. In the more general sense, the coefficient refers to the expansion coefficient of |L,M> angular momentum or spin state in the product basis of |L1,M1>|L2,M2> states. | |
| − | + | ==Syntax== | |
| − | + | ||
| − | + | cg=clebsch_gordan(L,M,L1,M1,L2,M2) | |
| − | + | ||
| − | + | ==Arguments== | |
| + | |||
| + | L,M,L1,M1,L2,M2 - integer or half-integer indices of | ||
| + | the angular momentum or spin states | ||
| + | |||
| + | ==Outputs== | ||
| + | |||
| + | cg - floating-point (double precision) | ||
| + | Clebsch-Gordan coefficient | ||
| + | |||
| + | ==Notes== | ||
| + | # Only some combinations of L,M,L1,M1,L2,M2 are allowed by the properties of spherical harmonics and spin states. If inadmissible indices are supplied, zero is returned. | ||
| + | # CG coefficient calculation in double-precision arithmetic is not a trivial matter for high ranks. This function produces fast answers with an accuracy of about 1e-3 up to about L=20. A slower machine precision implementation for higher ranks is available in [[clebsch_gordan.m]] function. | ||
| + | |||
| + | ==See also== | ||
| + | [[clebsch_gordan.m]], [[wigner.m]], [[wigner_3j.m]], [[wigner_6j.m]] | ||
| + | |||
| + | |||
| + | ''Version 2.3, authors: [[Ilya Kuprov]]'' | ||
Revision as of 12:35, 27 December 2018
Clebsch-Gordan coefficient: the coefficient in front of Y(L,M) spherical harmonic in the expansion of the product of Y(L1,M1) and Y(L2,M2) spherical harmonics. In the more general sense, the coefficient refers to the expansion coefficient of |L,M> angular momentum or spin state in the product basis of |L1,M1>|L2,M2> states.
Syntax
cg=clebsch_gordan(L,M,L1,M1,L2,M2)
Arguments
L,M,L1,M1,L2,M2 - integer or half-integer indices of
the angular momentum or spin states
Outputs
cg - floating-point (double precision)
Clebsch-Gordan coefficient
Notes
- Only some combinations of L,M,L1,M1,L2,M2 are allowed by the properties of spherical harmonics and spin states. If inadmissible indices are supplied, zero is returned.
- CG coefficient calculation in double-precision arithmetic is not a trivial matter for high ranks. This function produces fast answers with an accuracy of about 1e-3 up to about L=20. A slower machine precision implementation for higher ranks is available in clebsch_gordan.m function.
See also
clebsch_gordan.m, wigner.m, wigner_3j.m, wigner_6j.m
Version 2.3, authors: Ilya Kuprov
