<div dir="ltr"><div>CP2K is a perfectly fine force / energy calculator. Lack of k-points and symmetry might make things cumbersome. There is / used to be a wrapper for phonopy which can be used to drive the phonon calculations. Lev Kantorovich's Tetr / Lev00 also have some functionality.</div><div><br></div><div>Have a look at this:<br></div><div><br></div><div><a href="https://www.cp2k.org/exercises:2017_uzh_cmest:phonon_calculation">https://www.cp2k.org/exercises:2017_uzh_cmest:phonon_calculation</a><br><br>it might be useful.</div><div><br></div><div>Matt<br><br>On Tuesday, April 30, 2019 at 2:51:00 AM UTC+1, bobs your uncle wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">Hi folks,<div><br></div><div>Just a general question regarding vibrational analysis in an organic crystal using CP2K (just solid state DFT ignoring AIMD for now).</div><div>I thought about using more specialized software like quantum espresso - although I would prefer to use Crystal which I don't have access to currently - but the possibility of using CP2K for both solid state DFT, and later AIMD was more appealing to me for obvious reasons.</div><div>So first off, would you recommend CP2K for vibrational analysis of periodic systems? It would seem to me that it is uncommon, as I am having issues finding examples in the literature, and there aren't any guides either.</div><div><br></div><div>To get to grips with CP2K, I first completed some of the tutorials as recommended before moving onto a "toy" system, where I have been attempting to reinvent the wheel with vibrational analysis. </div><div>Of course, I would really like to know if anyone could advise me on the "right" way to do this. </div><div>Currently, I start with a cell optimisation on an experimental crystal structure - following the sodium chloride guide on the CP2K website - then reoptimize the geometry with tighter convergence threshholds etc. </div><div>I take that geometry and perform a vibrational analysis, and have succeeded in producing a complete set of real harmonics this way - which is a start - but I haven't been able to simulate some of the low frequency "terahertz" modes that have been reproduced in the literature (although for my system, the papers I can find don't talk at all about the software or process followed).</div><div>No one in my immediate academic circle is trained in these calculations, so I would really appreciate some of your wisdom if you could share it!</div><div><br></div><div>Sorry, that was a bit of a read, but I'll put input files up if you are interested.</div><div><br></div><div>Cheers</div></div></blockquote></div></div>