Garnet correlation coefficients

Comparing to correlation coefficients from William D. Carlson and Clark R. Wilson, Phys Chem Minerals 43, 363-369 (2016) doi:10.1007/s00269-016-0800-2

Garnet structure includes pyrope, which we use as our example structure, with space group 230 (Ia3d) with stoichiometry Mg3Al2Si3O12. The occupied Wyckoff positions for this are (lattice constant \(a_0\)=1.1459 nm):

Wykcoff site

chemistry

position

24c

Mg

1/8 0 1/4

16a

Al

0 0 0

24d

Si

3/8 0 1/4

96h

O

.03284 .05014 .65330

Data from G. V. Gibbs and J. V. Smith, “Refinement of the crystal structure of synthetic pyrope.” American Mineralogist 50 2023-2039 (1965), PDF.

[1]:
import sys
sys.path.extend(['../'])
import numpy as np
import onsager.crystal as crystal
import onsager.OnsagerCalc as onsager

Create garnet crystal (lattice constant in nm). Wyckoff positions cut and pasted from Bilbao crystallographic server.

[2]:
# a0 = 1.1459
# alatt = a0*np.eye(3)
a0 = 1.
alatt = a0*np.array([[-0.5,0.5,0.5],[0.5,-0.5,0.5],[0.5,0.5,-0.5]])
invlatt = np.array([[0,1,1],[1,0,1],[1,1,0]])
x,y,z = (.03284,.05014,.65330)
uMg = ((1/8,0,1/4),(3/8,0,3/4),(1/4,1/8,0),(3/4,3/8,0),
       (0,1/4,1/8),(0,3/4,3/8),(7/8,0,3/4),(5/8,0,1/4),
       (3/4,7/8,0),(1/4,5/8,0),(0,3/4,7/8),(0,1/4,5/8))
uAl = ((0,0,0),(1/2,0,1/2),(0,1/2,1/2),(1/2,1/2,0),
       (3/4,1/4,1/4),(3/4,3/4,3/4),(1/4,1/4,3/4),(1/4,3/4,1/4))
uSi = ((3/8,0,1/4),(1/8,0,3/4),(1/4,3/8,0),(3/4,1/8,0),
       (0,1/4,3/8),(0,3/4,1/8),(3/4,5/8,0),(3/4,3/8,1/2),
       (1/8,1/2,1/4),(7/8,0,1/4),(0,1/4,7/8),(1/2,1/4,1/8))
uO = ((x,y,z),(-x+1/2,-y,z+1/2),(-x,y+1/2,-z+1/2),(x+1/2,-y+1/2,-z),
      (z,x,y),(z+1/2,-x+1/2,-y),(-z+1/2,-x,y+1/2),(-z,x+1/2,-y+1/2),
      (y,z,x),(-y,z+1/2,-x+1/2),(y+1/2,-z+1/2,-x),(-y+1/2,-z,x+1/2),
      (y+3/4,x+1/4,-z+1/4),(-y+3/4,-x+3/4,-z+3/4),(y+1/4,-x+1/4,z+3/4),(-y+1/4,x+3/4,z+1/4),
      (x+3/4,z+1/4,-y+1/4),(-x+1/4,z+3/4,y+1/4),(-x+3/4,-z+3/4,-y+3/4),(x+1/4,-z+1/4,y+3/4),
      (z+3/4,y+1/4,-x+1/4),(z+1/4,-y+1/4,x+3/4),(-z+1/4,y+3/4,x+1/4),(-z+3/4,-y+3/4,-x+3/4),
      (-x,-y,-z),(x+1/2,y,-z+1/2),(x,-y+1/2,z+1/2),(-x+1/2,y+1/2,z),
      (-z,-x,-y),(-z+1/2,x+1/2,y),(z+1/2,x,-y+1/2),(z,-x+1/2,y+1/2),
      (-y,-z,-x),(y,-z+1/2,x+1/2),(-y+1/2,z+1/2,x),(y+1/2,z,-x+1/2),
      (-y+1/4,-x+3/4,z+3/4),(y+1/4,x+1/4,z+1/4),(-y+3/4,x+3/4,-z+1/4),(y+3/4,-x+1/4,-z+3/4),
      (-x+1/4,-z+3/4,y+3/4),(x+3/4,-z+1/4,-y+3/4),(x+1/4,z+1/4,y+1/4),(-x+3/4,z+3/4,-y+1/4),
      (-z+1/4,-y+3/4,x+3/4),(-z+3/4,y+3/4,-x+1/4),(z+3/4,-y+1/4,-x+3/4),(z+1/4,y+1/4,x+1/4))
# tovec = lambda x: np.array(x)
# tovec2 = lambda x: np.array((x[0]+1/2,x[1]+1/2,x[2]+1/2))
tovec = lambda x: np.dot(invlatt, x)
pyrope = crystal.Crystal(alatt, [[vec(w) for w in ulist for vec in (tovec,)]
                                 for ulist in (uMg, uAl, uSi, uO)],
                         ['Mg','Al','Si','O'])
# print(pyrope)

Next, we construct a diffuser based on vacancies for our Mg ion. We need to create a sitelist (which will be the Wyckoff positions) and a jumpnetwork for the transitions between the sites. There are tags that correspond to the unique states and transitions in the diffuser. The first cutoff is \(\sim 0.31a_0\), but that connects half of the Mg cation sites to each other; increasing the cutoff to \(\sim 0.51a_0\) introduces a second network that completes the connections.

[3]:
chem = 0  # 0 is the index corresponding to our Mg atom in the crystal
cutoff = 0.31*a0  # had been 0.51*a0
sitelist = pyrope.sitelist(chem)
jumpnetwork = pyrope.jumpnetwork(chem, cutoff)
Mgdiffuser = onsager.VacancyMediated(pyrope, chem, sitelist, jumpnetwork, 1)
print(Mgdiffuser)
Diffuser for atom 0 (Mg), Nthermo=1
#Lattice:
  a1 = [-0.5  0.5  0.5]
  a2 = [ 0.5 -0.5  0.5]
  a3 = [ 0.5  0.5 -0.5]
#Basis:
  (Mg) 0.0 = [ 0.25   0.375  0.125]
  (Mg) 0.1 = [ 0.75   0.125  0.375]
  (Mg) 0.2 = [ 0.125  0.25   0.375]
  (Mg) 0.3 = [ 0.375  0.75   0.125]
  (Mg) 0.4 = [ 0.375  0.125  0.25 ]
  (Mg) 0.5 = [ 0.125  0.375  0.75 ]
  (Mg) 0.6 = [ 0.75   0.625  0.875]
  (Mg) 0.7 = [ 0.25   0.875  0.625]
  (Mg) 0.8 = [ 0.875  0.75   0.625]
  (Mg) 0.9 = [ 0.625  0.25   0.875]
  (Mg) 0.10 = [ 0.625  0.875  0.75 ]
  (Mg) 0.11 = [ 0.875  0.625  0.25 ]
  (Al) 1.0 = [ 0.  0.  0.]
  (Al) 1.1 = [ 0.5  0.   0.5]
  (Al) 1.2 = [ 0.   0.5  0.5]
  (Al) 1.3 = [ 0.5  0.5  0. ]
  (Al) 1.4 = [ 0.5  0.   0. ]
  (Al) 1.5 = [ 0.5  0.5  0.5]
  (Al) 1.6 = [ 0.   0.   0.5]
  (Al) 1.7 = [ 0.   0.5  0. ]
  (Si) 2.0 = [ 0.25   0.625  0.375]
  (Si) 2.1 = [ 0.75   0.875  0.125]
  (Si) 2.2 = [ 0.375  0.25   0.625]
  (Si) 2.3 = [ 0.125  0.75   0.875]
  (Si) 2.4 = [ 0.625  0.375  0.25 ]
  (Si) 2.5 = [ 0.875  0.125  0.75 ]
  (Si) 2.6 = [ 0.625  0.75   0.375]
  (Si) 2.7 = [ 0.875  0.25   0.125]
  (Si) 2.8 = [ 0.75   0.375  0.625]
  (Si) 2.9 = [ 0.25   0.125  0.875]
  (Si) 2.10 = [ 0.125  0.875  0.25 ]
  (Si) 2.11 = [ 0.375  0.625  0.75 ]
  (O) 3.0 = [ 0.70344  0.68614  0.08298]
  (O) 3.1 = [ 0.10316  0.62046  0.41702]
  (O) 3.2 = [ 0.39684  0.81386  0.5173 ]
  (O) 3.3 = [ 0.79656  0.87954  0.9827 ]
  (O) 3.4 = [ 0.08298  0.70344  0.68614]
  (O) 3.5 = [ 0.41702  0.10316  0.62046]
  (O) 3.6 = [ 0.5173   0.39684  0.81386]
  (O) 3.7 = [ 0.9827   0.79656  0.87954]
  (O) 3.8 = [ 0.68614  0.08298  0.70344]
  (O) 3.9 = [ 0.62046  0.41702  0.10316]
  (O) 3.10 = [ 0.81386  0.5173   0.39684]
  (O) 3.11 = [ 0.87954  0.9827   0.79656]
  (O) 3.12 = [ 0.87954  0.39684  0.08298]
  (O) 3.13 = [ 0.81386  0.79656  0.41702]
  (O) 3.14 = [ 0.62046  0.70344  0.5173 ]
  (O) 3.15 = [ 0.68614  0.10316  0.9827 ]
  (O) 3.16 = [ 0.10316  0.9827   0.68614]
  (O) 3.17 = [ 0.70344  0.5173   0.62046]
  (O) 3.18 = [ 0.79656  0.41702  0.81386]
  (O) 3.19 = [ 0.39684  0.08298  0.87954]
  (O) 3.20 = [ 0.5173   0.62046  0.70344]
  (O) 3.21 = [ 0.9827   0.68614  0.10316]
  (O) 3.22 = [ 0.08298  0.87954  0.39684]
  (O) 3.23 = [ 0.41702  0.81386  0.79656]
  (O) 3.24 = [ 0.29656  0.31386  0.91702]
  (O) 3.25 = [ 0.89684  0.37954  0.58298]
  (O) 3.26 = [ 0.60316  0.18614  0.4827 ]
  (O) 3.27 = [ 0.20344  0.12046  0.0173 ]
  (O) 3.28 = [ 0.91702  0.29656  0.31386]
  (O) 3.29 = [ 0.58298  0.89684  0.37954]
  (O) 3.30 = [ 0.4827   0.60316  0.18614]
  (O) 3.31 = [ 0.0173   0.20344  0.12046]
  (O) 3.32 = [ 0.31386  0.91702  0.29656]
  (O) 3.33 = [ 0.37954  0.58298  0.89684]
  (O) 3.34 = [ 0.18614  0.4827   0.60316]
  (O) 3.35 = [ 0.12046  0.0173   0.20344]
  (O) 3.36 = [ 0.12046  0.60316  0.91702]
  (O) 3.37 = [ 0.18614  0.20344  0.58298]
  (O) 3.38 = [ 0.37954  0.29656  0.4827 ]
  (O) 3.39 = [ 0.31386  0.89684  0.0173 ]
  (O) 3.40 = [ 0.89684  0.0173   0.31386]
  (O) 3.41 = [ 0.29656  0.4827   0.37954]
  (O) 3.42 = [ 0.20344  0.58298  0.18614]
  (O) 3.43 = [ 0.60316  0.91702  0.12046]
  (O) 3.44 = [ 0.4827   0.37954  0.29656]
  (O) 3.45 = [ 0.0173   0.31386  0.89684]
  (O) 3.46 = [ 0.91702  0.12046  0.60316]
  (O) 3.47 = [ 0.58298  0.18614  0.20344]
vacancy configurations:
v:+0.250,+0.375,+0.125
solute configurations:
s:+0.250,+0.375,+0.125
solute-vacancy configurations:
s:+0.375,+0.125,+0.250-v:+0.750,+0.125,+0.375
omega0 jumps:
omega0:v:+0.625,+0.250,+0.875^v:+0.250,-0.125,+0.625
omega1 jumps:
omega1:s:+0.875,+0.625,+0.250-v:+0.625,+0.250,-0.125^v:+0.250,-0.125,-0.375
omega1:s:+0.750,+0.625,+0.875-v:+0.625,+0.250,+0.875^v:+0.250,-0.125,+0.625
omega1:s:+0.625,+0.875,+0.750-v:+0.625,+1.250,+0.875^v:+0.250,+0.875,+0.625
omega2 jumps:
omega2:s:+0.250,+0.875,+0.625-v:+0.625,+1.250,+0.875^s:+0.625,+0.250,+0.875-v:+0.250,-0.125,+0.625

Quick analysis on our jump network:

  1. What is the connectivity, \(Z\)?

  2. What is the individual contribution to \(\mathbf{\delta x}\otimes\mathbf{\delta x}\)? And 1/3 Tr (which will be the symmetrized contribution)?

  3. What is the squared magnitude \(\delta x^2\)?

[4]:
for jlist in jumpnetwork:
    Z = 0
    dx2 = np.zeros((3,3))
    for (i,j), dx in jlist:
        if i==0:
            Z += 1
            dx2 += np.outer(dx,dx)
    print("coordination number:", Z)
    print(dx2)
    print("1/3 Tr dx dx:", dx2.trace()/3)
    print("dx^2:", np.dot(dx,dx))
coordination number: 4
[[ 0.0625   0.       0.     ]
 [ 0.       0.15625 -0.125  ]
 [ 0.      -0.125    0.15625]]
1/3 Tr dx dx: 0.125
dx^2: 0.09375

Next, we assemble our data: the energies and prefactors, for a VMg in pyrope for our representative states and transitions: these are the first states in the lists, which are also identified by the tags above. As we are computing a tracer, we make the choice to set \(\nu_0 = 1/Z\) where \(Z=4\) is the coordination number.

[5]:
nu0 = 0.25
Etrans = 0.
# we don't need to use the tags, since there's only one site and jump type, and
# we want to build a tracer.
Mgthermodict = {'preV': np.ones(len(sitelist)),
                'eneV': np.zeros(len(sitelist)),
                'preT0': nu0*np.ones(len(jumpnetwork)),
                'eneT0': Etrans*np.ones(len(jumpnetwork))}
Mgthermodict.update(Mgdiffuser.maketracerpreene(**Mgthermodict))
for k,v in Mgthermodict.items():
    print('{}: {}'.format(k, v))
eneSV: [ 0.]
eneT0: [ 0.]
preS: [ 1.]
preT1: [ 0.25  0.25  0.25]
preSV: [ 1.]
preT0: [ 0.25]
eneV: [ 0.]
eneT2: [ 0.]
eneS: [ 0.]
preV: [ 1.]
eneT1: [ 0.  0.  0.]
preT2: [ 0.25]

We compute the Onsager matrices, and look at \(-L_\text{ss}/L_\text{sv}\) to get our correlation coefficient.

Note: we can define \(f\) (for our tracer) as the ratio of \(L_\text{ss}\) to \(Z (\delta x)^2 w_2 c_\text{v}c_\text{s}/6 = \frac{1}{16}\nu_0 a_0^2\) in this case, the same as what we get for \(L_\text{vv}\) and \(-L_\text{sv}\).

[6]:
Lvv, Lss, Lsv, L1vv = Mgdiffuser.Lij(*Mgdiffuser.preene2betafree(1, **Mgthermodict))
print(Lvv)
print(Lss)
print(Lsv)
print(L1vv)
print("Correlation coefficient:", -Lss[0,0]/Lsv[0,0])
[[ 0.015625 -0.       -0.      ]
 [-0.        0.015625 -0.      ]
 [-0.       -0.        0.015625]]
[[ 0.00585895  0.          0.        ]
 [ 0.          0.00585895  0.        ]
 [ 0.          0.          0.00585895]]
[[-0.015625  0.        0.      ]
 [ 0.       -0.015625  0.      ]
 [ 0.        0.       -0.015625]]
[[ -1.17108470e-34   0.00000000e+00   0.00000000e+00]
 [  0.00000000e+00  -1.17108470e-34   0.00000000e+00]
 [  0.00000000e+00   0.00000000e+00  -1.17108470e-34]]
Correlation coefficient: 0.374972670783

Compare with tabulated GF data from Carlson and Wilson paper. They use the notation \((l,m,n)\) for a \(\mathbf{\delta x}\) vector that is \(a_0(l\hat x+m\hat y+n\hat z)/8\). We will need to find a corresponding site that lands at that displacement from our origin site.

Unfortunately, it looks like in two cases ((800), (444)) there are two distinct sites that are mapped in that displacement vector, which have different GF values; the CW reported values appear to be the averaged values. In two other cases, ((640), (420)) the reported values are half of what the computed values are here.

As Carlson and Wilson used a stochastic approach to compute their GF values, all of their other data has errors \(\sim 10^{-4}\).

[7]:
# tabulated data from paper
CarlsonWilsonGFdata = \
{(0,0,0): 2.30796022, (2,1,1): 1.30807261, (3,3,2): 0.80669536,
 (4,2,0): 0.40469085, (4,4,4): 0.50242046, (5,3,2): 0.56195744,
 (6,1,1): 0.56071092, (6,4,0): 0.22460654, (6,5,3): 0.42028488,
 (6,5,5): 0.40137897, (7,2,1): 0.44437878, (8,0,0): 0.41938675}
[8]:
print('CW index\tdx match\tGF (FT eval)\tGF(CW stoch.)\terror')
GF = Mgdiffuser.GFcalc  # get our GF calculator; should already have rates set
basis = pyrope.basis[chem]
x0 = np.dot(alatt, basis[0])
for vec,gCW in CarlsonWilsonGFdata.items():
    dx0 = np.array(vec,dtype=float)/8
    nmatch, Gave, Gmatch = 0, 0, {}
    for g in pyrope.G:
        dx = np.dot(g.cartrot, dx0)
        j = pyrope.cart2pos(x0+dx)[1]
        if j is not None and j[0]==chem and j[1]<6:
            G = GF(0, j[1], dx)
            Gmatch[tuple((8*dx).astype(int))] = G
            nmatch += 1
            Gave += G
    Gave /= nmatch
    for t,G in Gmatch.items():
        print('{}\t{}\t{:.12f}\t{:.8f}\t{:.4e}'.format(vec, t, -G, gCW, abs(G+gCW)))
    print('{}\taverage value\t{:.12f}\t{:.8f}\t{:.4e}'.format(vec, -Gave, gCW, abs(Gave+gCW)))
CW index        dx match        GF (FT eval)    GF(CW stoch.)   error
(8, 0, 0)       (0, 0, -8)      0.427361034009  0.41938675      7.9743e-03
(8, 0, 0)       (8, 0, 0)       0.403566247455  0.41938675      1.5821e-02
(8, 0, 0)       (0, 8, 0)       0.427361034009  0.41938675      7.9743e-03
(8, 0, 0)       (-8, 0, 0)      0.403566247455  0.41938675      1.5821e-02
(8, 0, 0)       (0, -8, 0)      0.427361034009  0.41938675      7.9743e-03
(8, 0, 0)       (0, 0, 8)       0.427361034009  0.41938675      7.9743e-03
(8, 0, 0)       average value   0.419429438491  0.41938675      4.2688e-05
(6, 1, 1)       (-1, 6, 1)      0.560766700022  0.56071092      5.5780e-05
(6, 1, 1)       (-1, -6, -1)    0.560766700022  0.56071092      5.5780e-05
(6, 1, 1)       (1, 1, 6)       0.560766700022  0.56071092      5.5780e-05
(6, 1, 1)       (1, -1, -6)     0.560766700022  0.56071092      5.5780e-05
(6, 1, 1)       average value   0.560766700022  0.56071092      5.5780e-05
(3, 3, 2)       (-3, 3, -2)     0.806767995595  0.80669536      7.2636e-05
(3, 3, 2)       (3, -2, 3)      0.806767995595  0.80669536      7.2636e-05
(3, 3, 2)       (3, 2, -3)      0.806767995595  0.80669536      7.2636e-05
(3, 3, 2)       (-3, -3, 2)     0.806767995595  0.80669536      7.2636e-05
(3, 3, 2)       average value   0.806767995595  0.80669536      7.2636e-05
(2, 1, 1)       (-1, -2, 1)     1.308081132926  1.30807261      8.5229e-06
(2, 1, 1)       (-1, 2, -1)     1.308081132926  1.30807261      8.5229e-06
(2, 1, 1)       (1, 1, -2)      1.308081132926  1.30807261      8.5229e-06
(2, 1, 1)       (1, -1, 2)      1.308081132926  1.30807261      8.5229e-06
(2, 1, 1)       average value   1.308081132926  1.30807261      8.5229e-06
(6, 5, 3)       (3, -6, 5)      0.420386782427  0.42028488      1.0190e-04
(6, 5, 3)       (-3, -5, 6)     0.420386782427  0.42028488      1.0190e-04
(6, 5, 3)       (5, -3, -6)     0.420386782427  0.42028488      1.0190e-04
(6, 5, 3)       (-3, 5, -6)     0.420386782427  0.42028488      1.0190e-04
(6, 5, 3)       (-5, -6, -3)    0.420386782427  0.42028488      1.0190e-04
(6, 5, 3)       (-5, 6, 3)      0.420386782427  0.42028488      1.0190e-04
(6, 5, 3)       (5, 3, 6)       0.420386782427  0.42028488      1.0190e-04
(6, 5, 3)       (3, 6, -5)      0.420386782427  0.42028488      1.0190e-04
(6, 5, 3)       average value   0.420386782427  0.42028488      1.0190e-04
(6, 4, 0)       (-6, 0, 4)      0.449091350780  0.22460654      2.2448e-01
(6, 4, 0)       (6, 4, 0)       0.449091350780  0.22460654      2.2448e-01
(6, 4, 0)       (6, -4, 0)      0.449091350780  0.22460654      2.2448e-01
(6, 4, 0)       (-6, 0, -4)     0.449091350780  0.22460654      2.2448e-01
(6, 4, 0)       average value   0.449091350780  0.22460654      2.2448e-01
(0, 0, 0)       (0, 0, 0)       2.308081141615  2.30796022      1.2092e-04
(0, 0, 0)       average value   2.308081141615  2.30796022      1.2092e-04
(4, 2, 0)       (2, 0, -4)      0.809394258097  0.40469085      4.0470e-01
(4, 2, 0)       (-2, -4, 0)     0.809394258097  0.40469085      4.0470e-01
(4, 2, 0)       (-2, 4, 0)      0.809394258097  0.40469085      4.0470e-01
(4, 2, 0)       (2, 0, 4)       0.809394258097  0.40469085      4.0470e-01
(4, 2, 0)       average value   0.809394258097  0.40469085      4.0470e-01
(5, 3, 2)       (-5, 2, -3)     0.561961239416  0.56195744      3.7994e-06
(5, 3, 2)       (3, 2, 5)       0.561961239416  0.56195744      3.7994e-06
(5, 3, 2)       (3, -2, -5)     0.561961239416  0.56195744      3.7994e-06
(5, 3, 2)       (-5, -2, 3)     0.561961239416  0.56195744      3.7994e-06
(5, 3, 2)       (5, 3, -2)      0.561961239416  0.56195744      3.7994e-06
(5, 3, 2)       (5, -3, 2)      0.561961239416  0.56195744      3.7994e-06
(5, 3, 2)       (-3, -5, -2)    0.561961239416  0.56195744      3.7994e-06
(5, 3, 2)       (-3, 5, 2)      0.561961239416  0.56195744      3.7994e-06
(5, 3, 2)       average value   0.561961239416  0.56195744      3.7994e-06
(7, 2, 1)       (-1, -2, -7)    0.444350262895  0.44437878      2.8517e-05
(7, 2, 1)       (1, 7, 2)       0.444350262895  0.44437878      2.8517e-05
(7, 2, 1)       (7, 2, -1)      0.444350262895  0.44437878      2.8517e-05
(7, 2, 1)       (-7, 1, -2)     0.444350262895  0.44437878      2.8517e-05
(7, 2, 1)       (-7, -1, 2)     0.444350262895  0.44437878      2.8517e-05
(7, 2, 1)       (7, -2, 1)      0.444350262895  0.44437878      2.8517e-05
(7, 2, 1)       (-1, 2, 7)      0.444350262895  0.44437878      2.8517e-05
(7, 2, 1)       (1, -7, -2)     0.444350262895  0.44437878      2.8517e-05
(7, 2, 1)       average value   0.444350262895  0.44437878      2.8517e-05
(4, 4, 4)       (-4, 4, 4)      0.457297218361  0.50242046      4.5123e-02
(4, 4, 4)       (-4, 4, -4)     0.547635344309  0.50242046      4.5215e-02
(4, 4, 4)       (-4, -4, -4)    0.457297218361  0.50242046      4.5123e-02
(4, 4, 4)       (4, 4, -4)      0.547635344309  0.50242046      4.5215e-02
(4, 4, 4)       (-4, -4, 4)     0.547635344309  0.50242046      4.5215e-02
(4, 4, 4)       (4, -4, -4)     0.457297218361  0.50242046      4.5123e-02
(4, 4, 4)       (4, -4, 4)      0.547635344309  0.50242046      4.5215e-02
(4, 4, 4)       (4, 4, 4)       0.457297218361  0.50242046      4.5123e-02
(4, 4, 4)       average value   0.502466281335  0.50242046      4.5821e-05
(6, 5, 5)       (-5, 6, -5)     0.401425331863  0.40137897      4.6362e-05
(6, 5, 5)       (-5, -6, 5)     0.401425331863  0.40137897      4.6362e-05
(6, 5, 5)       (5, -5, 6)      0.401425331863  0.40137897      4.6362e-05
(6, 5, 5)       (5, 5, -6)      0.401425331863  0.40137897      4.6362e-05
(6, 5, 5)       average value   0.401425331863  0.40137897      4.6362e-05