Dear all,
I am trying to compare lattice energies of different zeolite structures with identical chemical compositions. For this I am generating their structures with symmetries and cation sites as reported from structural studies, but this implies split cation sites which are mutually exclusive (i.e. two closely spaced, symmetrically equivalent sites).
GULP runs just fine but the user manual specifically states that care must be taken to interpret the results in these situations, especially with partially occupancied, mutually exclusive sites. It further says that it may be necessary to define ‘exclude potentials’ to get the desired behaviour.
I have two questions:
- How do you define exclude potentials in this situation?
- Do you think it is fair to compare lattice energies calculated in this way (at least qualitatively), or does anyone have bad experiences when trying this with solid solutions?
Thanks a lot.
I add a typical input file; the split cation site which is mutually exclusive in this case would be ‘K3’
opti conp relax
name K_MER_167
dump every 1 K_MER_167.res
output cif K_MER_167.cif
cell
10.087400 14.031200 14.26750 90.0000 90.0000 90.0000
space
71
fractional
K1 core 0.000000 0.000000 0.215280 1.000000000 1.0000
K2 core 0.000000 0.352400 0.500000 1.000000000 1.0000
K3 core 0.322400 0.408000 0.000000 1.000000000 0.5000
O1 core 0.000000 0.371650 0.280970 0.869020000 1.0000
O2 core 0.000000 0.181540 0.122670 0.869020000 1.0000
O3 core 0.326000 0.000000 0.294130 0.869020000 1.0000
O4 core 0.206100 0.212150 0.000000 0.869020000 1.0000
O5 core 0.190730 0.319350 0.155470 0.869020000 1.0000
O6 core 0.245760 0.127450 0.166770 0.869020000 1.0000
O1 shel 0.000000 0.371650 0.280970 -2.869020000 1.0000
O2 shel 0.000000 0.181540 0.122670 -2.869020000 1.0000
O3 shel 0.326000 0.000000 0.294130 -2.869020000 1.0000
O4 shel 0.206100 0.212150 0.000000 -2.869020000 1.0000
O5 shel 0.190730 0.319350 0.155470 -2.869020000 1.0000
O6 shel 0.245760 0.127450 0.166770 -2.869020000 1.0000
Si1 core 0.346260 0.110120 0.255860 4.000000000 0.6250
Si2 core 0.158340 0.211340 0.111400 4.000000000 0.6250
Al1 core 0.346260 0.110120 0.255860 3.000000000 0.3750
Al2 core 0.158340 0.211340 0.111400 3.000000000 0.3750