- \sum_j \Phi_j^* \frac{\partial H}{\partial \vec{R}_i} \Phi_j
-+\sum_j \frac{\partial \Phi_j^*}{\partial \vec{R}_i} H \Phi_j
-+\sum_j \Phi_j^* H \frac{\partial \Phi_j}{\partial \vec{R}_i}
+ \sum_j \langle \Phi_j | \frac{\partial H}{\partial\vec{R}_i} | \Phi_j \rangle
++\sum_j \langle \frac{\partial \Phi_j}{\partial\vec{R}_i} | H \Phi_j \rangle
++\sum_j \langle \Phi_j H | \frac{\partial \Phi_j}{\partial \vec{R}_i} \rangle