[CP2K-user] [CP2K:10844] Bond Dissociation Energy Calculation

Krack Matthias (PSI) matthia... at psi.ch
Sun Oct 14 20:05:34 UTC 2018


Dear Aniruddha

UKS is certainly a good idea. There is also the possibility of an open shell singlet. I would check the literature concerning the choice of the functional.

Best

Matthias

Von: cp... at googlegroups.com <cp... at googlegroups.com> Im Auftrag von Aniruddha Dive
Gesendet: Sonntag, 14. Oktober 2018 01:13
An: cp2k <cp... at googlegroups.com>
Betreff: Re: [CP2K:10844] Bond Dissociation Energy Calculation

Thanks Matthais,

Does it mean that I need to use Unrestricted Kohn-Sham? Also currently I am doing a singlet calculations, should I do an triplet calculations for better results?

Is REVPBE okay or should I go ahead with B3LYP XC and BLYP pseudopotentials?

Regards,
Aniruddha M Dive

On Friday, October 12, 2018 at 12:49:40 PM UTC-7, Matthias Krack wrote:
Hi Aniruddha

It seems that you are using restricted Kohn-Sham (RKS) which constrains the state of the fragments. Consider the simple textbook case of H2. RKS allows only for the fragments H+ and H- if you pull the H atoms apart.

HTH

Matthias

Von: cp... at googlegroups.com<javascript:> <cp... at googlegroups.com<javascript:>> Im Auftrag von Aniruddha Dive
Gesendet: Freitag, 12. Oktober 2018 20:19
An: cp2k <cp... at googlegroups.com<javascript:>>
Betreff: [CP2K:10842] Bond Dissociation Energy Calculation

Hi All,

I am trying to calculate the bond dissociation energy of Si-S bond using CP2K. Please find attached my input file for the same.

I am initially performing a geometry optimization of a SiH3-SH molecule to calculate the equilibrium Si-S bond length. Once I have the equilibrium bond length, I am varying the Si-S distance from 1.5 A to 7.5 A. I rerun geometry optimization for these different structures keeping the Si and S atoms fixed while allowing the other H atoms to relax.

Kindly let me know if this input files and approach looks good. I am getting the bond dissociation energy to be 100 kcal/mol whereas the literature value is ~ 145 kcal/mol

Regards,
Aniruddha M Dive
Phd Candidate
Washington State University
--
You received this message because you are subscribed to the Google Groups "cp2k" group.
To unsubscribe from this group and stop receiving emails from it, send an email to cp2k+... at googlegroups.com<javascript:>.
To post to this group, send email to cp... at googlegroups.com<javascript:>.
Visit this group at https://groups.google.com/group/cp2k.
For more options, visit https://groups.google.com/d/optout.
--
You received this message because you are subscribed to the Google Groups "cp2k" group.
To unsubscribe from this group and stop receiving emails from it, send an email to cp2k+uns... at googlegroups.com<mailto:cp2k+uns... at googlegroups.com>.
To post to this group, send email to cp... at googlegroups.com<mailto:cp... at googlegroups.com>.
Visit this group at https://groups.google.com/group/cp2k.
For more options, visit https://groups.google.com/d/optout.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20181014/698dfdd1/attachment.htm>


More information about the CP2K-user mailing list