<div dir="ltr">Dear Jana, Jing,<div><br></div><div style>probably by now you have already solved the problem, but I was just faced with the same</div><div style>issue and I've just figured out what happened. Apparently, there was a minor change from 2.3, 2.4 (development version) and 2.4 (trunk version), my guess is that you're using (2.4 trunk). On this version the location of the subsections related to printing dipole moments and wannier centers have changed on the input file. Now, you don't need the LOCALIZATION section anymore, you can just start the subsection TOTAL_DIPOLE inside LOCALIZE%PRINT. </div>
<div style><br></div><div style>In your case, just delete the following lines:</div><div style><span style="font-family:arial,sans-serif;font-size:13px"> &LOCALIZATION</span><br></div><div style><span style="font-family:arial,sans-serif;font-size:13px">and</span></div>
<div style><span style="font-family:arial,sans-serif;font-size:13px">&END (refering to LOCALIZATION)</span><span style="font-family:arial,sans-serif;font-size:13px"><br></span></div><div style><span style="font-family:arial,sans-serif;font-size:13px"><br>
</span></div><div style><span style="font-family:arial,sans-serif;font-size:13px">and move &END LOCALIZE to after </span><span style="font-family:arial,sans-serif;font-size:13px">&END PRINT</span></div><div style>
<span style="font-family:arial,sans-serif;font-size:13px"><br></span></div><div style><span style="font-family:arial,sans-serif;font-size:13px">cheers,</span></div><div style><span style="font-family:arial,sans-serif;font-size:13px">Wagner.</span></div>
<div style><span style="font-family:arial,sans-serif;font-size:13px"><br></span></div><div style><font face="arial, sans-serif">ps. if you keep the old sintax the code won't complaint, it'll just run normally without computing it, maybe it's something for the developers to fix.</font></div>
<div style><span style="font-family:arial,sans-serif;font-size:13px"><br></span></div><div style><br></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Thu, Jul 4, 2013 at 10:15 PM, JQ <span dir="ltr"><<a href="mailto:gohjin...@gmail.com" target="_blank">gohjin...@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="im HOEnZb">Hi Jana,<div><br></div><div>I have just done a test run on my pc based on your input. Your input seems working fine that I manage to get the output of delta dipole moment and etc, saved into the files as defined. You may refer to the attached input file and output file (only few MD steps), for your reference. I am using CP2K 2.3 and you can find its details in the output file</div>
<div><br></div></div><span class="HOEnZb"><font color="#888888"><div>JQ</div></font></span><div class="HOEnZb"><div class="h5"><br>On Friday, June 28, 2013 3:28:39 AM UTC+3, Jana wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex">
Dear all,<div>I've been told I can get a total dipole file containing delta dipoles and these can be postprocessed via FT from a CPMD utility (and unit conversion) to obtain vibrational frequencies.</div><div><br></div>
<div>(This is an attempt to get a comparison as I've got reasonable IR spectrum simply from NORMAL_MODES where all the modes fit but for some reason the only one most characteristic frequency is strangely shifted like 500 cm-1 to low frequencies...this happens only in one of my several model systems. I've been trying to improve convergence criteria or basis sets or adding vdw correction but either no help or no convergence even after like more than one week on the wall...)</div>
<div><br></div><div>To get this I should run MD and print out TOTAL_DIPOLE, with reference to centers of atomic charges. I've tried to get dipoles and wannier centres according to what I've found while searching this google group, but I have only managed to get it run with ENERGY calculation, not MD, and as it is just a trial and error I got a bit lost and cannot really see what am I doing wrong. Could you please somebody help?</div>
<div><br></div><div><div>...</div><div>&MOTION</div><div> &MD</div><div> ENSEMBLE NVT</div><div> STEPS 1000</div><div> TIMESTEP 0.5</div><div> TEMPERATURE 300.0</div><div> &THERMOSTAT</div><div>
TYPE NOSE</div><div> REGION GLOBAL</div><div> &NOSE</div><div> LENGTH 3<br></div><div> YOSHIDA 3</div><div> MTS 2</div><div> &END</div><div> &END<br></div><div>
&END MD</div><div>&END MOTION<br></div><div><br></div><div>&FORCE_EVAL<br></div><div> METHOD QS</div><div> &DFT<br></div><div>&LOCALIZE<br></div><div><div> METHOD CRAZY</div><div> OPERATOR BERRY</div>
<div> &END LOCALIZE</div><div>&PRINT<br></div><div> &LOCALIZATION</div><div> &WANNIER_CENTERS MEDIUM</div><div> FILENAME =water2.wannier</div><div> &EACH</div><div>
MD 1</div><div> &END EACH</div><div> ADD_LAST<br></div><div> COMMON_ITERATION_LEVELS 3</div><div> &END</div><div><br></div><div> &MOLECULAR_DIPOLES HIGH</div>
<div> &END MOLECULAR_DIPOLES</div><div><br></div><div> &TOTAL_DIPOLE MEDIUM</div><div> FILENAME =water2.dip</div><div> &EACH</div><div> MD 1</div><div>
&END EACH</div><div> REFERENCE COM</div><div> PERIODIC TRUE</div><div> COMMON_ITERATION_LEVELS 3</div><div> &END TOTAL_DIPOLE</div><div> &END</div><div>&END PRINT</div>
</div><div><br></div><div> BASIS_SET_FILE_NAME GTH_BASIS_SETS<br></div><div> POTENTIAL_FILE_NAME GTH_POTENTIALS</div><div> &MGRID</div><div> CUTOFF [Ry] 450</div><div> &END MGRID</div><div> &QS</div>
<div> METHOD GPW</div><div> &END QS</div><div> &SCF</div><div> SCF_GUESS RESTART</div><div> MAX_SCF 100</div><div> EPS_SCF 1.0e-7</div><div> &OUTER_SCF</div><div> MAX_SCF 5</div>
<div> EPS_SCF 1.0E-7</div><div> &END</div><div> &END SCF</div><div> &XC</div><div> &XC_FUNCTIONAL PBE</div><div> &END XC_FUNCTIONAL</div><div> &XC_GRID<br></div>
<div> XC_SMOOTH_RHO NN50</div><div> XC_DERIV NN50_SMOOTH</div><div> &END XC_GRID</div><div> &END XC</div><div>&END DFT<br></div><div>...</div></div><div><br></div><div>I've also tried it with no specification of filename or each or common iteration levels, also for nonperiodic case with MT poisson solver but I never get any file or part of the output mentioning dipole or wannier centres. Quite sure the mistake is something really stupid, but I've been just spending few nights on it and still cannot figure it out. </div>
<div>Any help would be much appreciated!</div><div>J.</div><div><br></div></blockquote>
<p></p>
-- <br>
You received this message because you are subscribed to the Google Groups "cp2k" group.<br>
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp2k%2Bun...@googlegroups.com" target="_blank">cp2k+uns...@googlegroups.com</a>.<br>
To post to this group, send email to <a href="mailto:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a>.<br>
Visit this group at <a href="http://groups.google.com/group/cp2k" target="_blank">http://groups.google.com/group/cp2k</a>.<br>
For more options, visit <a href="https://groups.google.com/groups/opt_out" target="_blank">https://groups.google.com/groups/opt_out</a>.<br>
<br>
<br>
</div></div></blockquote></div><br></div>