<div dir="ltr">Dear cp2k developers,<div><br></div><div>I have a question about the sign change of atomic virial between qs_dftb_coulomb.F and qs_dftb_matrices.F for 7.1 version.  In general, the atomic virial is -0.5*r_{ij}*f_{ij} and it seems to be the convention (i.e. negative sign) used in qs_dftb_matrices.F.  However, it is the opposite in qs_dftb_coulomb.F.</div><div><br></div><div>For example, lines 234 to 239 read:</div><div><div>               fi = 1.0_dp</div><div>               IF (iatom == jatom) fi = 0.5_dp</div><div>               CALL virial_pair_force(virial%pv_virial, fi, fij, rij)</div><div>               IF (atprop%stress) THEN</div><div>                  CALL virial_pair_force(atprop%atstress(:, :, iatom), fi*0.5_dp, fij, rij)</div><div>                  CALL virial_pair_force(atprop%atstress(:, :, jatom), fi*0.5_dp, fij, rij)</div></div><div>Other places in qs_dftb_coulomb.F also used the positive sign.</div><div><br></div><div>I am wondering if the coulombic force was using a different sign or it was something else I am missing.</div><div><br></div><div>Thanks, Wei</div></div>