<div dir="ltr"><div>Dear developers,</div><div><br></div><div>it seems  that there is a bug when using tersoff potentials with more than one element.</div><div>The SUBROUTINE read_tersoff_section in file  force_fields_input.F fills the cp2k variables of the various potentials pairs,</div><div>always with the same (i.e. the first) data set of tersoff potentials read from the input file.</div><div>I think that the lines of the kind:</div><div>CALL section_vals_val_get(tersoff_section, "A", r_val=nonbonded%pot(start+isec)%pot%set(1)%tersoff%A)</div><div>should be changed in </div><div>CALL section_vals_val_get(tersoff_section, "A", i_rep_section=isec,</div><div>r_val=nonbonded%pot(start+isec)%pot%set(1)%tersoff%A)</div><div><br></div><div>Could you please check if I am correct ? (note that if this is correct the same should be done also for siepmann potentials).</div><div>I have attached the input file I have used for this test.</div><div>Thank you for your kind help.</div><div><br></div><div>Best regards,</div><div>Michele Gusso</div></div>