[CP2K-user] Binding Enery Curve for Si-S bond

[Aniruddha Dive]

Hi All,

I want to generate a binding energy curve for a Si-S bond with bond lengths 
varying from 1.5 A - 7.5 A. To calculate the binding energy curve I have 
generated a SiH3-SH molecule where I vary the Si-S bond lengths from 1.5 A 
to 7.5 A over 25 different geometries. I run a geometry optimization for 
each geometry keeping the Si and S atoms fixed while allowing other atoms 
to relax. Based on the energies for all these geometries I obtain a binding 
energy curve for Si-S which flattens out giving a bond dissociation energy 
of 105 kcal/mol. The experimental Si-S bond dissociation energy reported is 
~ 145 kcal/mol which is significantly higher than what I obtain. To 
generate the binding energy curve I take all the energies to be relative to 
the minima which corresponds to zero energy.

I am thinking that this difference could be due to the BSSE effects when 
the Si-S atoms are at a larger separation leading to generation of two 
fragments SiH3 and SH. To determine the BSSE, I ran a BSSE calculations for 
the SiH3-SH molecule where the Si-S distance is 7.5 A. I obtain a 
CP-corrected total energy from these calculatons. Now, while generating the 
binding energy curve do I need to use the CP-corrected total energy as the 
total energy for that geometry? Or should I subtract this CP-corrected 
total energy from the total energy obtained during geometry optimization 
which will give me the BSSE. Further should I add the BSSE obtained to the 
relative energies for these geometries where the Si-S distance is larger 
than 3.0 A.

I have attached input file providing details about the GEO_OPT as well as 
BSSE and the coordinate file for one case.

Kindly let me know how can I generate an accurate binding energy curve for 
these molecules?

Thanks,
Aniruddha M Dive
PhD Candidate, Washington State University
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20181017/73758ba7/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: bond_diss.inp
Type: chemical/x-gamess-input
Size: 3858 bytes
Desc: not available
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20181017/73758ba7/attachment.inp>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Si_S_conf26.pdb
Type: application/vnd.palm
Size: 1045 bytes
Desc: not available
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20181017/73758ba7/attachment.bin>