I am getting a good result.

Since the psf contains water and DPC, then I should not use ?nres because it include waters.

However the timeseries is still influenced by the intermittent breakup by reaching the boundary of the unit cell.
But it is easier to exclude those values.