Apparently cp2k assigns the same covalent radius for the MM atoms.<br>I took the input file from tests/QMMM/QS/regtest-4/wat_nacl.inp and changed the print level to medium.<br>This is what I get:<br><br><font face="courier new,monospace"> MM POINT CHARGES GENERATING THE QM/MM ELECTROSTATIC POTENTIAL<br>
<br> -------------------------------------------------------------------------<br> MM ATOM: 3 RADIUS: 1.511781 CHARGE: 0.417000CORR. RADIUS 1.511781<br> MM ATOM: 4 RADIUS: 1.511781 CHARGE: -0.834000CORR. RADIUS 1.511781<br>
MM ATOM: 5 RADIUS: 1.511781 CHARGE: 0.417000CORR. RADIUS 1.511781<br> MM ATOM: 6 RADIUS: 1.511781 CHARGE: 0.417000CORR. RADIUS 1.511781<br> MM ATOM: 7 RADIUS: 1.511781 CHARGE: -0.834000CORR. RADIUS 1.511781<br>
MM ATOM: 8 RADIUS: 1.511781 CHARGE: 0.417000CORR. RADIUS 1.511781<br> MM ATOM: 9 RADIUS: 1.511781 CHARGE: 0.417000CORR. RADIUS 1.511781<br> MM ATOM: 10 RADIUS: 1.511781 CHARGE: 1.000000CORR. RADIUS 1.511781<br>
MM ATOM: 11 RADIUS: 1.511781 CHARGE: -1.000000CORR. RADIUS 1.511781<br><br></font><font face="arial,helvetica,sans-serif">Atom 3 is an oxygen, while 4 if a hydrogen. Is this related to periodic QM/MM?<br>
Thanks<br>Marius<br></font>