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