NT /
3nfJuelich
Back up to 3b
Channel definitions used by the Juelich group
The following is a sample output from using the 3nf_read.py script:
- The script reads a 3bME file from Epelbaum and prints some information.
- Note in particular the channel definitions and numbering.
- This particular example is the J=1/2, T=1/2, pi=+ channel (v1p5.chiral-lec-c-laminf-p16n5-1-v4-1+P.T=1_2.form)
- Note that there is a two-body cutoff j_{NN} <=5
% python 3nf_read.py
Reading file v3n.form
Momentum mesh:
np12 = 16 pp = [0.001, 0.08096, 0.32084, 0.72064, 1.28036, 2.0, 2.57143, 3.14286, 3.71429, 4.28571, 4.85714, 5.42857, 6.0, 7.33333, 8.66667, 10.0]
np3 = 16 qp = [0.001, 0.08096, 0.32084, 0.72064, 1.28036, 2.0, 3.42857, 4.85714, 6.28571, 7.71429, 9.14286, 10.5714, 12.0, 14.6667, 17.3333, 20.0]
NN Channels:
j cutoff: j_NN <= 5
# l s j t
--- -- -- -- --
0 0 0 0 1
1 1 1 0 1
2 1 0 1 0
3 0 1 1 0
4 2 1 1 0
5 1 1 1 1
6 2 0 2 1
7 1 1 2 1
8 3 1 2 1
9 2 1 2 0
10 3 0 3 0
11 2 1 3 0
12 4 1 3 0
13 3 1 3 1
14 4 0 4 1
15 3 1 4 1
16 5 1 4 1
17 4 1 4 0
18 5 0 5 0
19 4 1 5 0
20 6 1 5 0
21 5 1 5 1
NNN channels:
3-body quantum numbers: 2J = 1 Pi = 1 2T = 1
# l_12 s_12 j_12 t_12 l_3 2I_3 2J 2T
--- --- --- --- --- --- --- --- ---
0 0 0 0 1 0 1 1 1
1 1 1 0 1 1 1 1 1
2 1 0 1 0 1 1 1 1
3 1 0 1 0 1 3 1 1
4 0 1 1 0 0 1 1 1
5 0 1 1 0 2 3 1 1
6 2 1 1 0 0 1 1 1
7 2 1 1 0 2 3 1 1
8 1 1 1 1 1 1 1 1
9 1 1 1 1 1 3 1 1
10 2 0 2 1 2 3 1 1
11 2 0 2 1 2 5 1 1
12 1 1 2 1 1 3 1 1
13 1 1 2 1 3 5 1 1
14 3 1 2 1 1 3 1 1
15 3 1 2 1 3 5 1 1
16 2 1 2 0 2 3 1 1
17 2 1 2 0 2 5 1 1
18 3 0 3 0 3 5 1 1
19 3 0 3 0 3 7 1 1
20 2 1 3 0 2 5 1 1
21 2 1 3 0 4 7 1 1
22 4 1 3 0 2 5 1 1
23 4 1 3 0 4 7 1 1
24 3 1 3 1 3 5 1 1
25 3 1 3 1 3 7 1 1
26 4 0 4 1 4 7 1 1
27 4 0 4 1 4 9 1 1
28 3 1 4 1 3 7 1 1
29 3 1 4 1 5 9 1 1
30 5 1 4 1 3 7 1 1
31 5 1 4 1 5 9 1 1
32 4 1 4 0 4 7 1 1
33 4 1 4 0 4 9 1 1
34 5 0 5 0 5 9 1 1
35 5 0 5 0 5 11 1 1
36 4 1 5 0 4 9 1 1
37 4 1 5 0 6 11 1 1
38 6 1 5 0 4 9 1 1
39 6 1 5 0 6 11 1 1
40 5 1 5 1 5 9 1 1
41 5 1 5 1 5 11 1 1
Number of channels: 22 (NN) 42 (NNN)
Number of MEs per alpha,alpha' channel: 65536
Reading matrix elements!
- Only channel (alpha',alpha)= (0, 0) will be read.
- The rest of the file will be scanned.
Storing momentum grid data (q3',p12', q3, p12) in file: pq_grid.dat