<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 trying to calculate diffusion barriers of single atoms adsorbed on a surface using the NEB method.</div><div>
 I optimized at DFT level the initial and final positions of the 
diffusion hop and I am providing those as first and last replicas.<br></div><div>I have some questions about the set up of the calculation and the computational strategy to follow to obtain good results.</div><div> <br></div><div>1)
 My system contains more than 200 atoms. Can I fix the atoms underneath 
the first layer of the surface to save computational time?</div><div>2) Should the value of K_SPRING be varied in order to see if there are changes in the results?<br></div><div>3) Should the value of NUMBER_OF_REPLICAS be varied in order to see if there are changes in the results?</div><div>4) Are there any other important parameters to take care of?<br></div><div></div><div>5) When
 plotting the energy profile of the NEB calculations I always find that 
the last frame, which I gave as an input, changes in energy with 
respect to the result found in the geometry optimization previously 
performed.Is it due to the ALIGN_FRAMES and ROTATE_FRAMES keywords, which I did not switched off? <br></div><div><br></div><div>Is this a good strategy to obtain the diffusion energy barriers or would you suggest something else?<br></div><div><br></div><div>I've attached the input file used and the geometries of the initial and final replicas.</div><div>Thank you for your help!</div><div><br></div><div>Stefano Ferrero<br></div></div>