<div>Dear Jörg,</div><div>maybe, I did not understand  the question; </div><div>if you "decouple the cell" with MT scheme you are going to model the system in gas phase like in GAMESS.</div><div>Anyway, if you do not use periodic system, you have to check that electron density goes to zero prior to the cell border.</div><div>When I study a reaction with cp2k, and I need to  characterize the TS, I use dimer-method or NEB. </div><div>If I need to model the system in gas phase</div><div>(no periodic system) and the system is "spheric-like", </div><div>I prefer to use  WAVELET.</div><div>I hope it is helpful to you.<br></div><div>Cheers,</div><div><br></div>Luca<div><br></div><div><br><div><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">Il giorno domenica 13 dicembre 2020 alle 12:26:24 UTC+1 sassy ha scritto:<br></div><blockquote class="gmail_quote" style="margin: 0 0 0 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">Dear all,
<br>
<br>I am normally calculating molecules in the gas phase, so my knowledge of cp2k 
<br>is still a but rusty as I don't use it as much as I would love to. 
<br>
<br>I am currently looking into an observed solid state reaction of the type:
<br>
<br>A -> B + C
<br>
<br>We got an x-ray of A and we observe the product B + C per x-ray as well as a ~ 
<br>10% "contamination".
<br>If I am using the usual DFT tools like GAMESS-US or Gaussian and put the 
<br>coordinates of A + B in the input file, they are flying away from each other. 
<br>That makes sense in the gas phase but not in the solid state. 
<br>Thus, I thought of calculating the unit cell in cp2k which should mimic better 
<br>the conditions where we are observing the products. 
<br>
<br>I was thinking of using PBE-D3 for that, together with the Molopt basis set 
<br>which was working well for me in other calculations I done before. 
<br>
<br>I am a bit unsure regarding the affect the size of the cell and which 
<br>decoupler to use. I usually go for MT but as there has been quite some 
<br>development recently it seems, I am not sure if that is still the best current 
<br>advice. 
<br>
<br>For example, what will happen if for my specific problem I make the cell a bit 
<br>bigger than what you usually would do for the chosen decoupler? 
<br>
<br>The overall aim of the project is more ambitious, I want to model the reaction 
<br>path of this reaction if that is possible. However, that is for the future. 
<br>
<br>Any suggestions are much appreciated.
<br>
<br>All the best from Lyttle Britain.
<br>
<br>Jörg
<br>
<br>
<br>
<br></blockquote></div></div>