R中Debye散射方程对Au纳米粒子的衍射规律。

[poster]

我想通过使用德拜散射方程绘制金纳米颗粒的衍射图。我有包含原子位置的数据。这是等式，I（q）= f 2 ∑ i N ∑ j N sin（2 pi r ij q）/（2 pi r ij q）。 i = j = 1，N =原子数。

我尝试自己做，但获得的图与参考图不同。我认为我在代码的数学部分出现了错误。如果有人帮助我，我会很高兴。 read_delim("/Users/batyr/Downloads/0.151.dat",col_names = F,delim=" ") %>% rename(x=X1,y=X2,z=X3) %>% as.matrix() ->atomic_arrangement # vector position of atoms, x,y,z. x