ISISNeutronMuon / MDANSE

MDANSE: Molecular Dynamics Analysis for Neutron Scattering Experiments
https://www.isis.stfc.ac.uk/Pages/MDANSEproject.aspx
GNU General Public License v3.0
18 stars 4 forks source link

Correlation update for Angular correlation, DOS, MSD, PACF and VACF #444

Closed ChiCheng45 closed 1 month ago

ChiCheng45 commented 1 month ago

Description of work Applied the correlation update to some dynamics jobs following #435. The approach that I'm that is proposed in the PR is similar to the correlation function defined in mcquarrie statistical mechanics.

image

An alternative truncation-type procedure would be different from the above (before taking the limit) since the window size $-T$ to $T$ decreases for increasing $\tau$.

A new option in the frames setting has been added "correlation frames". This gives you the resulting number of frames in your correlations function. E.g. With VACF if "correlation frames" is set to 5 you will have 5 time steps in your VACF.

image

For the default setting I set the correlation frames used to around half of the total frames available.

Angular Correlation Spurious oscillation at long time are removed.

Before image

After image

DOS Overall noise of the DOS is reduced.

Before image

After image

MSD Some small oscillation at long time are removed. The new MSD appears linear at long times while the old scheme doesn't look straight for long times.

Before image

After image

PACF Spurious oscillation at long time are removed.

Before image

After image

VACF Spuriously increasing correlation a long times are removed.

Before image

After image

To test Run one of the dynamics jobs and check everything is working correctly. Try changing the number of correlation frames and check that it works as expected.