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>