Dear CP2K community,<div><br /></div><div>I'm trying to learn how to calculate polarization using CP2K. I'm using 2x2x2 supercell of BaTiO3 to benchmark the polarization values compared to other references and Quantum Espresso. </div><div><br /></div><div>I'm using Berry phase method by setting <font color="#000000"><b>PERIODIC .TRUE. in MOMENTS</b><b> </b></font>as it's periodic system. And the output for the cubic structure looks like:</div><div><font color="#808080">  Dipole vectors are based on the periodic (Berry phase) operator.<br />  They are defined modulo integer multiples of the cell matrix [Debye].<br />  [X] [ Â  38.42563418 Â  Â  0.00000000 Â  Â  0.00000000 ] [i]<br />  [Y]=[ Â  Â 0.00000000 Â  Â 38.42563418 Â  Â  0.00000000 ]*[j]<br />  [Z] [ Â  Â 0.00000000 Â  Â  0.00000000 Â  Â </font><font color="#008000">38.65618798</font><font color="#808080"> ] [k]<br />  Dipole moment [Debye]<br />  Â  X= -0.1123396E-09 Y= -0.1123351E-09 Z= </font><font color="#008000">-0.1414722E-09</font><font color="#808080"> Â Total= Â  Â  Â 0.0000000</font><br /></div><div><br /></div><div>How I have been calculating the polarization value is: 1) As I'm only interested in the polarization in z direction, I combined the two values highlighted in green. And 2) I divide the values by the cell volume. 3) I convert the unit to uC/cm^2. The issue with this method is that the polarization of the cubic phase is not zero. I've been trying to figure this out because the polarization value for cubic phase should be zero or close to zero. And the zero polarization for cubic phase is reproduced by Quantum Espresso:</div><div>  Â  Â <font color="#808080"> Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â VALUES OF POLARIZATION<br />  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â ~~~~~~~~~~~~~~~~~~~~~~<br />  Â  Â  Â  The calculation of phases done along the direction of vector 3<br />  Â  Â  Â  of the reciprocal lattice gives the following contribution to<br />  Â  Â  Â  the polarization vector (in different units, and being Omega<br />  Â  Â  Â  the volume of the unit cell):<br />  Â  Â  Â  Â  Â P = Â  0.0000000 Â (mod Â 30.0953409) Â (e/Omega).bohr<br />  Â  Â  Â  Â  Â P = Â  0.0000000 Â (mod Â  0.0088327) Â e/bohr^2<br />  Â  Â  Â  Â  Â P = Â  0.0000000 Â (mod Â  0.5053585) Â C/m^2<br />  Â  Â  Â  The polarization direction is: Â ( 0.00000 , 0.00000 , 1.00000 )</font><br /></div><div><br /></div><div><br /></div><div>What I'm confused with is that the difference between the values in the dipole vectors and the dipole moment. Is either of them electronic polarization (from Berry phase) and ionic polarization? Or is the dipole moment in the bottom line calculated from the addition of the electronic polarization (Berry phase) and the ionic polarization? Is there any reference/literature that describes how CP2K operates the Berry phase calculation?</div><div><br /></div><div><br /></div><div>Here's my input script:</div><div><font color="#808080">&GLOBAL<br />  PROJECT BTO<br />  RUN_TYPE CELL_OPT<br />  PRINT_LEVEL MEDIUM<br />&END GLOBAL<br /><br />&FORCE_EVAL<br />  ! This one is used for cell optimization<br />  STRESS_TENSOR ANALYTICAL<br />  METHOD Â QS<br />  &SUBSYS<br />  Â  &KIND Ba<br />  Â  Â  ELEMENT Â  Ba<br />  Â  Â  BASIS_SET DZVP-MOLOPT-PBE-GTH-q10<br />  Â  Â  POTENTIAL GTH-PBE-q10<br />  Â  &END KIND<br />  Â  &KIND Ti<br />  Â  Â  ELEMENT Â  Ti<br />  Â  Â  BASIS_SET DZVP-MOLOPT-PBE-GTH-q12<br />  Â  Â  POTENTIAL GTH-PBE-q12<br />  Â  &END KIND<br />  Â  &KIND O<br />  Â  Â  ELEMENT Â  O<br />  Â  Â  BASIS_SET DZVP-MOLOPT-PBE-GTH-q6<br />  Â  Â  POTENTIAL GTH-PBE-q6<br />  Â  &END KIND<br />  Â  &CELL<br />  Â  Â  A Â  Â  8.000000000 Â  Â 0.000000000 Â  Â 0.000000000<br />  Â  Â  B Â  Â  0.000000000 Â  Â 8.000000000 Â  Â 0.000000000<br />  Â  Â  C Â  Â  0.000000000 Â  Â 0.000000000 Â  Â 8.000000000<br />  Â  Â  PERIODIC XYZ<br />  Â  &END CELL<br />  Â  &TOPOLOGY<br />  Â  Â  COORD_FILE_FORMAT pdb<br />  Â  Â  COORD_FILE_NAME BTO.pdb<br />  Â  &END<br />  &END SUBSYS<br />  &DFT<br />  Â  BASIS_SET_FILE_NAME  BASIS_MOLOPT_UZH<br />  Â  POTENTIAL_FILE_NAME  GTH_POTENTIALS<br />  Â  &QS<br />  Â  Â  METHOD GPW<br />  Â  Â  EPS_DEFAULT 1.0E-12<br />  Â  &END QS<br />  Â  &POISSON<br />  Â  Â  PERIODIC XYZ<br />  Â  Â  PSOLVER PERIODIC<br />  Â  &END<br />  Â  &XC<br />  Â  Â  &XC_FUNCTIONAL<br />  Â  Â  Â  Â &PBE<br />  Â  Â  Â  Â &END PBE<br />  Â  Â  &END<br />  Â  Â  &VDW_POTENTIAL<br />  Â  Â  Â  Â POTENTIAL_TYPE PAIR_POTENTIAL<br />  Â  Â  Â  Â &PAIR_POTENTIAL<br />  Â  Â  Â  Â  Â  PARAMETER_FILE_NAME ../../dftd3.dat<br />  Â  Â  Â  Â  Â  TYPE DFTD3(BJ)<br />  Â  Â  Â  Â  Â  REFERENCE_FUNCTIONAL PBE<br />  Â  Â  Â  Â &END<br /></font></div><div><font color="#808080">  Â  Â  &END VDW_POTENTIAL<br />  Â  &END XC<br />  Â  &MGRID<br />  Â  Â  NGRIDS 4<br />  Â  Â  CUTOFF 600<br />  Â  &END MGRID<br />  Â  &SCF<br />  Â  Â  SCF_GUESS ATOMIC<br />  Â  Â  EPS_SCF 1.0E-8<br />  Â  Â  MAX_SCF 300<br />  Â  Â  &MIXING<br />  Â  Â  Â  Â METHOD BROYDEN_MIXING<br />  Â  Â  Â  Â ALPHA 0.2<br />  Â  Â  Â  Â BETA 1.5<br />  Â  Â  Â  Â NBROYDEN 8<br />  Â  Â  &END MIXING<br />  Â  &END SCF<br />  Â  &PRINT<br />  Â  Â  &MO_CUBES<br />  Â  Â  Â  WRITE_CUBE .FALSE.<br />  Â  Â  Â  NHOMO 1<br />  Â  Â  Â  NLUMO 1<br />  Â  Â  &END MO_CUBES<br />  Â  Â  &MOMENTS MEDIUM<br />  Â  Â  Â  PERIODIC T<br />  Â  Â  &END MOMENTS<br />  Â  &END PRINT<br />  &END DFT<br />  &PRINT<br />  Â  &FORCES ON<br />  Â  &END FORCES<br />  &END PRINT<br />&END FORCE_EVAL<br /><br />&MOTION<br />  &CELL_OPT<br />  Â  TYPE DIRECT_CELL_OPT<br />  Â  MAX_ITER 300<br />  Â  OPTIMIZER BFGS Â  Â #Most efficient minimizer, but only for 'small' systems<br />  &END CELL_OPT<br />  &GEO_OPT<br />  Â  MAX_DR 0.003<br />  Â  MAX_FORCE 0.00045<br />  Â  RMS_DR 0.0015<br />  Â  RMS_FORCE 0.0003<br />  Â  MAX_ITER 200<br />  Â  OPTIMIZER BFGS<br />  Â  &BFGS<br />  Â  Â  TRUST_RADIUS 0.25<br />  Â  &END BFGS<br />  &END GEO_OPT<br />&END MOTION</font><br /></div><div><br /></div><div><div>Here's the input structure:</div><div><font color="#808080">CRYST1 Â  Â 8.000 Â  Â 8.000 Â  Â 8.000 Â 90.00 Â 90.00 Â 90.00 P 1 Â  Â  Â  Â  Â  1<br />ATOM Â  Â  Â 1 Â Ba Â BTO A Â  1 Â  Â  Â  0.000 Â  0.000 Â  0.000 Â 1.00 Â 0.00 Â  Â  Â  Â  Â Ba<br />ATOM Â  Â  Â 2 Â Ti Â BTO A Â  1 Â  Â  Â  2.000 Â  2.000 Â  2.000 Â 1.00 Â 0.00 Â  Â  Â  Â  Â Ti<br />ATOM Â  Â  Â 3 Â O Â  BTO A Â  1 Â  Â  Â  0.000 Â  2.000 Â  2.000 Â 1.00 Â 0.00 Â  Â  Â  Â  Â  O<br />ATOM Â  Â  Â 4 Â O Â  BTO A Â  1 Â  Â  Â  2.000 Â  0.000 Â  2.000 Â 1.00 Â 0.00 Â  Â  Â  Â  Â  O<br />ATOM Â  Â  Â 5 Â O Â  BTO A Â  1 Â  Â  Â  2.000 Â  2.000 Â  0.000 Â 1.00 Â 0.00 Â  Â  Â  Â  Â  O<br />ATOM Â  Â  Â 6 Â Ba Â BTO A Â  1 Â  Â  Â  0.000 Â  0.000 Â  4.000 Â 1.00 Â 0.00 Â  Â  Â  Â  Â Ba<br />ATOM Â  Â  Â 7 Â Ti Â BTO A Â  1 Â  Â  Â  2.000 Â  2.000 Â  6.000 Â 1.00 Â 0.00 Â  Â  Â  Â  Â Ti<br />ATOM Â  Â  Â 8 Â O Â  BTO A Â  1 Â  Â  Â  0.000 Â  2.000 Â  6.000 Â 1.00 Â 0.00 Â  Â  Â  Â  Â  O<br />ATOM Â  Â  Â 9 Â O Â  BTO A Â  1 Â  Â  Â  2.000 Â  0.000 Â  6.000 Â 1.00 Â 0.00 Â  Â  Â  Â  Â  O<br />ATOM Â  Â  10 Â O Â  BTO A Â  1 Â  Â  Â  2.000 Â  2.000 Â  4.000 Â 1.00 Â 0.00 Â  Â  Â  Â  Â  O<br />ATOM Â  Â  11 Â Ba Â BTO A Â  1 Â  Â  Â  0.000 Â  4.000 Â  0.000 Â 1.00 Â 0.00 Â  Â  Â  Â  Â Ba<br />ATOM Â  Â  12 Â Ti Â BTO A Â  1 Â  Â  Â  2.000 Â  6.000 Â  2.000 Â 1.00 Â 0.00 Â  Â  Â  Â  Â Ti<br />ATOM Â  Â  13 Â O Â  BTO A Â  1 Â  Â  Â  0.000 Â  6.000 Â  2.000 Â 1.00 Â 0.00 Â  Â  Â  Â  Â  O<br />ATOM Â  Â  14 Â O Â  BTO A Â  1 Â  Â  Â  2.000 Â  4.000 Â  2.000 Â 1.00 Â 0.00 Â  Â  Â  Â  Â  O<br />ATOM Â  Â  15 Â O Â  BTO A Â  1 Â  Â  Â  2.000 Â  6.000 Â  0.000 Â 1.00 Â 0.00 Â  Â  Â  Â  Â  O<br />ATOM Â  Â  16 Â Ba Â BTO A Â  1 Â  Â  Â  0.000 Â  4.000 Â  4.000 Â 1.00 Â 0.00 Â  Â  Â  Â  Â Ba<br />ATOM Â  Â  17 Â Ti Â BTO A Â  1 Â  Â  Â  2.000 Â  6.000 Â  6.000 Â 1.00 Â 0.00 Â  Â  Â  Â  Â Ti<br />ATOM Â  Â  18 Â O Â  BTO A Â  1 Â  Â  Â  0.000 Â  6.000 Â  6.000 Â 1.00 Â 0.00 Â  Â  Â  Â  Â  O<br />ATOM Â  Â  19 Â O Â  BTO A Â  1 Â  Â  Â  2.000 Â  4.000 Â  6.000 Â 1.00 Â 0.00 Â  Â  Â  Â  Â  O<br />ATOM Â  Â  20 Â O Â  BTO A Â  1 Â  Â  Â  2.000 Â  6.000 Â  4.000 Â 1.00 Â 0.00 Â  Â  Â  Â  Â  O<br />ATOM Â  Â  21 Â Ba Â BTO A Â  1 Â  Â  Â  4.000 Â  0.000 Â  0.000 Â 1.00 Â 0.00 Â  Â  Â  Â  Â Ba<br />ATOM Â  Â  22 Â Ti Â BTO A Â  1 Â  Â  Â  6.000 Â  2.000 Â  2.000 Â 1.00 Â 0.00 Â  Â  Â  Â  Â Ti<br />ATOM Â  Â  23 Â O Â  BTO A Â  1 Â  Â  Â  4.000 Â  2.000 Â  2.000 Â 1.00 Â 0.00 Â  Â  Â  Â  Â  O<br />ATOM Â  Â  24 Â O Â  BTO A Â  1 Â  Â  Â  6.000 Â  0.000 Â  2.000 Â 1.00 Â 0.00 Â  Â  Â  Â  Â  O<br />ATOM Â  Â  25 Â O Â  BTO A Â  1 Â  Â  Â  6.000 Â  2.000 Â  0.000 Â 1.00 Â 0.00 Â  Â  Â  Â  Â  O<br />ATOM Â  Â  26 Â Ba Â BTO A Â  1 Â  Â  Â  4.000 Â  0.000 Â  4.000 Â 1.00 Â 0.00 Â  Â  Â  Â  Â Ba<br />ATOM Â  Â  27 Â Ti Â BTO A Â  1 Â  Â  Â  6.000 Â  2.000 Â  6.000 Â 1.00 Â 0.00 Â  Â  Â  Â  Â Ti<br />ATOM Â  Â  28 Â O Â  BTO A Â  1 Â  Â  Â  4.000 Â  2.000 Â  6.000 Â 1.00 Â 0.00 Â  Â  Â  Â  Â  O<br />ATOM Â  Â  29 Â O Â  BTO A Â  1 Â  Â  Â  6.000 Â  0.000 Â  6.000 Â 1.00 Â 0.00 Â  Â  Â  Â  Â  O<br />ATOM Â  Â  30 Â O Â  BTO A Â  1 Â  Â  Â  6.000 Â  2.000 Â  4.000 Â 1.00 Â 0.00 Â  Â  Â  Â  Â  O<br />ATOM Â  Â  31 Â Ba Â BTO A Â  1 Â  Â  Â  4.000 Â  4.000 Â  0.000 Â 1.00 Â 0.00 Â  Â  Â  Â  Â Ba<br />ATOM Â  Â  32 Â Ti Â BTO A Â  1 Â  Â  Â  6.000 Â  6.000 Â  2.000 Â 1.00 Â 0.00 Â  Â  Â  Â  Â Ti<br />ATOM Â  Â  33 Â O Â  BTO A Â  1 Â  Â  Â  4.000 Â  6.000 Â  2.000 Â 1.00 Â 0.00 Â  Â  Â  Â  Â  O<br />ATOM Â  Â  34 Â O Â  BTO A Â  1 Â  Â  Â  6.000 Â  4.000 Â  2.000 Â 1.00 Â 0.00 Â  Â  Â  Â  Â  O<br />ATOM Â  Â  35 Â O Â  BTO A Â  1 Â  Â  Â  6.000 Â  6.000 Â  0.000 Â 1.00 Â 0.00 Â  Â  Â  Â  Â  O<br />ATOM Â  Â  36 Â Ba Â BTO A Â  1 Â  Â  Â  4.000 Â  4.000 Â  4.000 Â 1.00 Â 0.00 Â  Â  Â  Â  Â Ba<br />ATOM Â  Â  37 Â Ti Â BTO A Â  1 Â  Â  Â  6.000 Â  6.000 Â  6.000 Â 1.00 Â 0.00 Â  Â  Â  Â  Â Ti<br />ATOM Â  Â  38 Â O Â  BTO A Â  1 Â  Â  Â  4.000 Â  6.000 Â  6.000 Â 1.00 Â 0.00 Â  Â  Â  Â  Â  O<br />ATOM Â  Â  39 Â O Â  BTO A Â  1 Â  Â  Â  6.000 Â  4.000 Â  6.000 Â 1.00 Â 0.00 Â  Â  Â  Â  Â  O<br />ATOM Â  Â  40 Â O Â  BTO A Â  1 Â  Â  Â  6.000 Â  6.000 Â  4.000 Â 1.00 Â 0.00 Â  Â  Â  Â  Â  O<br />END</font><br /></div><div><br /></div><div><br /></div><div>Any of your inputs would be greatly appreciated. Thank you for your time and consideration in advance.</div><div><br /></div><div>Best, </div><div>Sue</div></div>

<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+unsubscribe@googlegroups.com">cp2k+unsubscribe@googlegroups.com</a>.<br />
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/d69b786c-2b3e-47be-bada-caef6b123f7dn%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/d69b786c-2b3e-47be-bada-caef6b123f7dn%40googlegroups.com</a>.<br />