<div dir="ltr">Hi,<div><br></div><div>I think, for neutral systems, if you are running a slab, then the intensity in the non-periodic direction should be OK - make sure you centre the slab in the cell. There is no way of reliably getting the in-plane intensities atm, to my knowledge.</div><div><br></div><div>Whether things are reliable for your charged system, I really do not know. In principle the dipole is not so well defined in this case, only the full multipole expansion ... but changes in dipole might be OK.<br><br>Matt<br><br>On Friday, November 25, 2016 at 6:16:19 AM UTC, Chanwoo Noh wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">Dear all,<div><br></div><div>I'm trying to calculate the vibrational spectrum for a periodic metal slab system with the finite electronic temperature smearing option.</div><div>There are two methods for calculating dipole moment in CP2K.</div><div>One is the Berry phase formula method for periodic systems, and the other is the simple operator method for non-periodic systems.</div><div>However, Berry phase formula method is not applicable to the system with smearing option.</div><div>Thus, I used the simple operator option for calculating dipole moment in the periodic system.</div><div>Is it okay that the simple operator method is used for periodic system?</div><div>If it is not, how can I calculate the dipole moment of the periodic system with smearing option?</div><div>Please give me any comment or advice.</div><div><br></div><div><br></div><div>Chanwoo Noh<br></div><div><br></div><div><br></div><div>Following is my input file.</div><div><br></div><div><div>&GLOBAL</div><div>  PROJECT Pt_H3O </div><div>  RUN_TYPE VIBRATIONAL_ANALYSIS</div><div>  PRINT_LEVEL MEDIUM</div><div>  EXTENDED_FFT_LENGTHS</div><div>&END GLOBAL</div><div>&FORCE_EVAL</div><div>  METHOD Quickstep</div><div>  &DFT</div><div>    BASIS_SET_FILE_NAME  BASIS_MOLOPT</div><div>    POTENTIAL_FILE_NAME  GTH_POTENTIALS</div><div>    &PRINT</div><div>        &MOMENTS</div><div>            PERIODIC FALSE </div><div>        &END</div><div>    &END</div><div>    &MGRID</div><div>      CUTOFF 500</div><div>      REL_CUTOFF 50</div><div>    &END MGRID</div><div>    &QS</div><div>      EPS_DEFAULT 1.0E-14</div><div>    &END QS</div><div><br></div><div>    CHARGE 1</div><div>    MULTIPLICITY 1 </div><div><br></div><div>    &SCF</div><div>      SCF_GUESS RESTART</div><div>      MAX_SCF   500</div><div>      EPS_SCF   5.0E-7</div><div>      ADDED_MOS 500</div><div>      CHOLESKY INVERSE </div><div>      &SMEAR ON</div><div>          METHOD FERMI_DIRAC</div><div>          ELECTRONIC_TEMPERATURE [K] 300</div><div>      &END SMEAR</div><div>      &DIAGONALIZATION</div><div>          ALGORITHM STANDARD</div><div>      &END DIAGONALIZATION</div><div>      &MIXING</div><div>          METHOD BROYDEN_MIXING</div><div>          ALPHA 0.1</div><div>          BETA 1.5</div><div>          NBROYDEN 8</div><div>      &END</div><div><br></div><div>      &PRINT</div><div>          &RESTART</div><div>              &EACH</div><div>                  QS_SCF 0</div><div>                  GEO_OPT 1</div><div>              &END</div><div>              ADD_LAST NUMERIC</div><div>              FILENAME RESTART</div><div>          &END</div><div>          &RESTART_HISTORY OFF</div><div>          &END</div><div>      &END</div><div>    &END SCF</div><div><br></div><div>    &XC</div><div>      &XC_FUNCTIONAL PBE</div><div>      &END XC_FUNCTIONAL</div><div><br></div><div>      &vdW_POTENTIAL</div><div>          POTENTIAL_TYPE PAIR_POTENTIAL</div><div>          &PAIR_POTENTIAL</div><div>              TYPE DFTD3</div><div>              CALCULATE_C9_TERM .TRUE.</div><div>              PARAMETER_FILE_NAME dftd3.dat</div><div>              REFERENCE_FUNCTIONAL PBE</div><div>              R_CUTOFF [angstrom] 10.0</div><div>          &END PAIR_POTENTIAL</div><div>      &END vdW_POTENTIAL</div><div><br></div><div>    &END XC</div><div>  &END DFT</div><div><br></div><div>  &SUBSYS</div><div>    &KIND Pt</div><div>      ELEMENT<span style="white-space:pre">        </span>Pt</div><div>      BASIS_SET DZVP-MOLOPT-SR-GTH</div><div>      POTENTIAL GTH-PBE-q18</div><div>    &END KIND</div><div>    &KIND H</div><div>      ELEMENT<span style="white-space:pre"> </span>H</div><div>      BASIS_SET DZVP-MOLOPT-GTH</div><div>      POTENTIAL GTH-PBE-q1</div><div>    &END KIND</div><div>    &KIND O</div><div>      ELEMENT<span style="white-space:pre">      </span>O</div><div>      BASIS_SET DZVP-MOLOPT-GTH</div><div>      POTENTIAL GTH-PBE-q6</div><div>    &END KIND</div><div>    </div><div> </div><div>    &CELL</div><div>      A     16.752000000    0.000000000    0.000000000</div><div>      B     8.376000000    14.50766000    0.000000000</div><div>      C     0.000000000    0.000000000    70.00000000</div><div>    &END CELL</div><div>    &TOPOLOGY</div><div>        COORD_FILE_NAME     Pt_H3O.xyz</div><div>        COORD_FILE_FORMAT   XYZ</div><div>    &END TOPOLOGY</div><div>  &END SUBSYS</div><div>  &PRINT</div><div>    &TOTAL_NUMBERS  ON</div><div>    &END TOTAL_NUMBERS</div><div>  &END PRINT</div><div>&END FORCE_EVAL</div><div><br></div><div>&MOTION</div><div>    &GEO_OPT</div><div>        OPTIMIZER   LBFGS</div><div>        MAX_FORCE   0.00005</div><div>        MAX_ITER    1600</div><div>    &END</div><div>    &CONSTRAINT</div><div>        &FIXED_ATOMS</div><div>            COMPONENTS_TO_FIX XYZ</div><div>            LIST 1..144</div><div>        &END</div><div>    &END</div><div>&END MOTION</div><div><br></div><div>&VIBRATIONAL_ANALYSIS</div><div><br></div><div>    INTENSITIES</div><div>    NPROC_REP 72 </div><div>    DX 0.001</div><div>    &PRINT</div><div>        &PROGRAM_RUN_INFO ON</div><div>        &END</div><div>    &END</div><div><br></div><div>&END</div></div><div><br></div><div><br></div><div><br></div></div></blockquote></div></div>