<div dir="ltr">Thank you very much.<br><br>在 2018年3月14日星期三 UTC+8下午7:43:41，ganta....@gmail.com写道：<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">Hi Zuohen,<div><br></div><div>I think it is not implemented in CP2K yet (not sure), so I am using different water model that does not require coulombic subtraction. </div><div><br></div><div>Have a nice day.</div><div><br></div><div>Thanks and Regards,</div><div>Prasanth. <br><br>On Monday, March 12, 2018 at 1:07:22 PM UTC+1, <a>zuoh...@gmail.com</a> wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">I have the same question. Did you implement the coulombic subtraction for morse potential successfully?<br><br>在 2017年11月20日星期一 UTC+8上午6:38:15，<a>ganta.pra...@<wbr>gmail.com</a>写道：<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Dear CP2K team/users,<div><br></div><div>I am trying to use FIST to do GEO_OPT of a box of water molecules. I have following doubts.</div><div><br></div><div>1) How to implement coulombic subtraction in &BONDED --> Morse potential.</div><div><br></div><div>In CP2K the morse potential is in the form</div><div><br></div><div><em style="color:rgb(0,0,0);font-family:"Times New Roman";font-size:medium">K1*[(1-exp(-K2*(R-R0)))^2-1])</em><span style="color:rgb(0,0,0);font-family:"Times New Roman";font-size:small"> </span><br></div><div><span style="color:rgb(0,0,0);font-family:"Times New Roman";font-size:small"><br></span></div><div><span style="color:rgb(0,0,0);font-family:"Times New Roman";font-size:small">But I want to include coulombic subtraction and implement the below form. </span></div><div><br></div><div><em style="color:rgb(0,0,0);font-family:"Times New Roman";font-size:medium">K1*[(1-exp(-K2*(R-R0)))^2-1]) - D*Q1*Q2/R</em><br></div><div><em style="color:rgb(0,0,0);font-family:"Times New Roman";font-size:medium"><br></em></div><div><span style="color:rgb(0,0,0)"><font face="arial, sans-serif" size="2">where Q1 and Q2 are charges of the atoms.</font></span></div><div><span style="color:rgb(0,0,0)"><font face="arial, sans-serif" size="2"><br></font></span></div><div><span style="color:rgb(0,0,0)"><font face="arial, sans-serif" size="2">Please let me know if it is possible. </font></span></div><div><span style="color:rgb(0,0,0)"><font face="arial, sans-serif" size="2"><br></font></span></div><div><span style="color:rgb(0,0,0)"><font face="arial, sans-serif" size="2">2) What will be the units of the FUNCTION for below case. Will it be eV or CP2K internal  units? Do we need to adjust the FUNCTION expression manually to match the CP2K internal units (atomic units) ? Please let me know. </font></span></div><div><span style="color:rgb(0,0,0)"><font face="arial, sans-serif" size="2"><br></font></span></div><div><span style="color:rgb(0,0,0)"><font face="arial, sans-serif" size="2"> </font></span><font color="#000000" face="arial, sans-serif" size="2">        &GENPOT                                                                                                                        </font></div><div><font color="#000000" face="arial, sans-serif" size="2">          ATOMS OWAT OWAT                                                                                                              </font></div><div><font color="#000000" face="arial, sans-serif" size="2">          FUNCTION ((A/(r^12))-(B/(r^6)))                                                                                              </font></div><div><font color="#000000" face="arial, sans-serif" size="2">          PARAMETERS A B                                                                                                               </font></div><div><font color="#000000" face="arial, sans-serif" size="2">          UNITS eV*angstrom^12 eV*angstrom^6                                                                                           </font></div><div><font color="#000000" face="arial, sans-serif" size="2">          VALUES  39344.98  42.15                                                                                                      </font></div><div><font color="#000000" face="arial, sans-serif" size="2">          VARIABLES r                                                                                                                  </font></div><div><font color="#000000" face="arial, sans-serif" size="2">          RCUT [angstrom] $WATERCUTOFF                                                                                                 </font></div><div><font color="#000000" face="arial, sans-serif" size="2">        &END GENPOT     </font></div><div><font color="#000000" face="arial, sans-serif" size="2"><br></font></div><div><font color="#000000" face="arial, sans-serif" size="2">Have a nice day.</font></div><div><font color="#000000" face="arial, sans-serif" size="2"><br></font></div><div><font color="#000000" face="arial, sans-serif" size="2">Thanks and Regards,</font></div><div><font color="#000000" face="arial, sans-serif" size="2">Prasanth.                                                 </font></div></div></blockquote></div></blockquote></div></div></blockquote></div>