Open ChiCheng45 opened 2 weeks ago
I just answered in the other issue. From a physical point of view, the minimum Q that makes sense in the simulation is the (1,0,0) vector. Here we calculate S(Q) as the FT of the g(r), so the Q-range chosen is arbitrary, but any Q value below that Qmin is not meaningful. In particular, for Q=0 the exp term is 1 and you will get a number that is proportional to the sum of all the scattering lengths of the atoms in your system.
Description of the error Why does the static structure factor have a large negative value for q->0. Here is one I generated using a cp2k water trajectory. I'd imagine there are numerical issues as q=0 but the drop to a large negative value begins before it gets close to 0.
Here is an experimental static structure factor of liquid Ar (https://journals.aps.org/pra/abstract/10.1103/PhysRevA.7.2130).
it goes to zero as q->0.