[CP2K-user] CI-NEB calculation: how to control the spin of every replica?

Stefano Ferrero stefano.... at edu.unito.it
Wed Jan 29 17:02:53 UTC 2020


Dear CP2K experts,

I have just started using CP2K and I would be really thankful if anybody 
can help me.
I am studying the reaction between an oxygen atom and an hydrogen atom on 
an ice surface. I have optimized the stationary points (O and H separated 
as the initial state and the final product OH) at a PBE-D3 level. 
Now I want to find the TS of the reaction and I am doing a CI-NEB 
calculation.
For the reaction to occur, if I am not mistaken, the multiplicity of the 
system should be 2, e.g. with the oxygen atom with its two unpaired 
electrons spin up and the hydrogen atom with its unpaired electron spin 
down or vice versa and I accomplished it using a broken symmetry scheme. 

Will the broken symmetry approach be applied also in the replicas generated 
by the program or not?
if not, how can I control the spin situation in the replicas that are 
generated by the program? 
Are there any other strategies to find the TS avoiding the NEB calculation 
but remaining sure about the spin situation of the system?

I have attached my input, initial and final state.

Thank you for your help
Stefano Ferrero 
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20200129/39bbc917/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: initial_state.xyz
Type: chemical/x-xyz
Size: 10418 bytes
Desc: not available
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20200129/39bbc917/attachment.xyz>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: final_state.xyz
Type: chemical/x-xyz
Size: 10416 bytes
Desc: not available
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20200129/39bbc917/attachment-0001.xyz>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: NEB_O_hydrogenation_1_replicas8.inp
Type: chemical/x-gamess-input
Size: 3132 bytes
Desc: not available
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20200129/39bbc917/attachment.inp>


More information about the CP2K-user mailing list