<div>Good day.</div><div><br /></div><div>I wish to bring under attention </div><div>of the developers </div><div>(of this approach</div><div></div><div></div><div>https://manual.cp2k.org/trunk/methods/qm_mm/image_charges.html)</div><div>two examples, </div><div>where </div><div>their code, apparently, fails to compute</div><div>optical electronic transitions.</div><div><br /></div><div>Specifically, following the online examples,</div><div>I condition OH neutral molecule at four gold layers (facet 111),</div><div>where gold atoms are fixed under hexagonal symmetry</div><div>to provide a stable, well-defined structural reference.</div><div><br /></div><div>The neutral OH molecule has one unpaired electron.</div><div>Therefore multiplicity of the system is 2.</div><div><br /></div><div>My coworker and me, we run time-dependent DFT to compute three transitions</div><div>under regular setting</div><div>and under SAOP setting.</div><div><br /></div><div>Here, I attach inputs and outputs for both cases </div><div>and the coordinate file.</div><div><br /></div><div>If you wish to try to run the jobs yourself,</div><div>you would need Au.pot file, which is present on the CP2K website.</div><div>and you would need to comment out the restart wavefunction </div><div>and reset charges to ATOMIC. </div><div><br /></div><div>My co-worker and me tested the jobs </div><div>under linux cp2k of different installations,</div><div>including the most recent one.</div><div><br /></div><div>The results are the same:</div><div>1) SAOP job dies with the same</div><div> *** WARNING in qs_tddfpt2_utils.F:634 :: Orbital energy correction ***<br /> *** potential is an experimental feature. Use it with extreme care ***</div><div>2) regular TDDFT job dies without any notice.</div><div><br /></div><div>If possible, </div><div>would you be kind to look at the case and comment,</div><div>if the failure(s) can be amended.</div><div><br /></div><div>If the CP2K QMMM-IMAGE theory cannot solve such cases,</div><div>would you instruct kindly</div><div>if the recent CP2K version got an update for TD-DFT</div><div>to compute optical transitions </div><div>to include only a narrow sub-set of orbitals, which may be selected</div><div>to be specific to HOMO-LUMO properties of the adsorbed molecule</div><div>rather than to the metal continuum sp-d transitions.   </div><div><br /></div><div>If such TD-DFT is possible, would you be kind to hint</div><div>how to instruct this in input file.</div><div><br /></div><div>Thank you for your time and attention.</div><div><br /></div><div>Sincerely, </div><div>I wish you a wonderful and peaceful time </div><div>during the spring festival of life revival.</div><div><br /></div><div>Victor</div><div><br /></div>

<p></p>

-- <br />
You received this message because you are subscribed to the Google Groups "cp2k" group.<br />
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp2k+unsubscribe@googlegroups.com">cp2k+unsubscribe@googlegroups.com</a>.<br />
To view this discussion visit <a href="https://groups.google.com/d/msgid/cp2k/6967f5d1-7230-460b-a506-0991c7d67540n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/6967f5d1-7230-460b-a506-0991c7d67540n%40googlegroups.com</a>.<br />