Код решетки FCC 100 и 111

Question

Я пытаюсь сделать кристаллическую решетку FCC 100 и FCC 111, используя код C. Я настроил решетку FCC 100, и я сделал это правильно. Однако, когда я попытался сделать решетку FCC 111, она выглядела иначе, чем я ожидал. Я использую эталонную решетку FCC 111 (созданную другой программой с графическим интерфейсом "Амстердамская функция плотности"), чтобы сравнить мою решетку FCC 111 (используя код C). Мой вопрос: параметры, которые я использовал, чтобы сделать FCC (111) в моем коде, верны? Если они не верны, пожалуйста, помогите мне построить решетку FCC 111. заранее спасибо

/**** To create FCC 100 lattice*********/
for(i=0; i<Xatom; i++){            // Number of Atoms in the X direction
for(j=0; j<Yatom; j++){        // Number of Atoms in the Y direction
    for(k=0; k<nL; k++) {       // Number of layers in the Z direction

x[i]= I * a_lattice;   y[j]= j * a_lattice;   z[k]= k * a_lattice;
fprintf(fpo," %s\t %lf\t %lf\t %lf\n", AtomName, x[i], y[j], z[k]);

x[i]= I * a_lattice;   y[j]= 0.5 * a_lattice +  j * a_lattice;    z[k]= 0.5 * a_lattice + k * a_lattice;
fprintf(fpo," %s\t %lf\t %lf\t %lf\n", AtomName, x[i], y[j], z[k]);

x[i]= 0.5 * a_lattice + I * a_lattice;  y[j]= j * a_lattice;   z[k]= 0.5 * a_lattice + k *a_lattice;
fprintf(fpo," %s\t %lf\t %lf\t %lf\n", AtomName, x[i], y[j], z[k]);

x[i]= 0.5 * a_lattice + i*a_lattice;  y[j] = 0.5 * a_lattice + j * a_lattice;   z[k]= k * a_lattice;
fprintf(fpo," %s\t %lf\t %lf\t %lf\n", AtomName, x[i], y[j], z[k]);
        }   
    }
}

/**** To create FCC 111 lattice*********/
ax = a_lattice * sqrt(2)/2;
ay = a_lattice * sqrt(6)/2;
az = a_lattice * sqrt(3);

x2 = sqrt(2)/4 * a_lattice;
y2 = sqrt(6)/4 * a_lattice;
y3 = sqrt(6)/6 * a_lattice;
y4 = sqrt(6)*5/12 * a_lattice;
y5 = sqrt(6)*2/6 * a_lattice;
y6 = sqrt(6)/12 * a_lattice;

for(i=0; i<Xatom; i++){
for(j=0; j<Yatom; j++){
  layer = 0;
for(k=0; k<nL; k++){

x[i] = I * ax;  y[j] = j * ay;  z[k]= layer/ 3.0 * az;
fprintf(fpo," %s\t %lf\t %lf\t %lf\n", AtomName, x[i], y[j], z[k]);


x[i] = x2 + I * ax;   y[j] = y2 + j * ay;   z[k] = layer/3.0 * az;
fprintf(fpo," %s\t %lf\t %lf\t %lf\n", AtomName, x[i], y[j], z[k]);

layer = layer + 1;

x[i] = I * ax;  y[j] = y3 + j*ay;  z[k] = layer/3.0 * az;
fprintf(fpo," %s\t %lf\t %lf\t %lf\n", AtomName, x[i], y[j], z[k]);

x[i] = x2 + I * ax;   y[j] = y4 + j * ay;   z[k]=layer/3.0 * az;
fprintf(fpo," %s\t %lf\t %lf\t %lf\n", AtomName, x[i], y[j], z[k]);

layer = layer + 1;

x[i] = I * ax;  y[j] = y5 + j * ay;     z[k] = layer/3.0 * az;
fprintf(fpo," %s\t %lf\t %lf\t %lf\n", AtomName, x[i], y[j], z[k]);

x[i] = x2 + I * ax;   y[j] = y6 + j * ay;   z[k] = layer/3.0 * az;
fprintf(fpo," %s\t %lf\t %lf\t %lf\n", AtomName, x[i], y[j], z[k]);
layer = layer + 1; 
        }
    }
}

Здесь я также показал справочную структуру и мою структуру вывода из моего кода

Answer 1

Было бы полезно, если бы кто-то мог взглянуть на код C, который используется для создания этой решетки FCC (111) с прямоугольным поперечным сечением, и предложить, что можно изменить в этом коде, чтобы сделать решетку FCC (111) с поперечное сечение алмаза, аналогичное GUI-REF, представленному в этом вопросе.