water_dist = md.compute_distances(self.traj[current_frame],water_pairs)[0]# unit in nmN_QR.append(self.single_frame_N_QR(Q,R,water_dist))# equation (4) in ref [2]err_Sn_QR = np.std(N_QR)/len(N_QR)ifQ ==0:# equation (7) in ref [2]Sn_QR = np.mean(N_QR)-4./3....