<div dir="ltr"><div>Dear CP2K experts,</div><div><br></div><div>I have just started using CP2K and I would be really thankful if anybody can help me.</div><div>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. <br></div><div>Now I want to find the TS of the reaction and I am doing a CI-NEB calculation.</div><div> </div><div>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. <br></div><div><br></div><div>Will the broken symmetry approach be applied also in the replicas generated by the program or not?</div><div>if not, how can I control the spin situation in the replicas that are generated by the program? <br></div><div>Are there any other strategies to find the TS avoiding the NEB calculation but remaining sure about the spin situation of the system?<br></div><div><br></div><div>I have attached my input, initial and final state.<br></div><div><br></div><div>Thank you for your help</div><div>Stefano Ferrero <br></div></div>