<div dir="ltr">I need to post it again because I forgot to add input files and the version number.<br>Dear QE Users,<div><div><br>I am a PhD student working with Berry phase polarization calculations of single layer 2D material hBN in quantum espresso (applying vacuum in the z-direction).<br><br>I found -0.0000002 C/m^2 polarization in the z-direction in Quantum ESPRESSO 7.0 without electric field; however, total phases belonging to electronic and ionic contribution found as equal to each other and given as 0.40381 and -0.40381 with 5 significant digits.<br>Is it OK to take -0.0000002 C/m^2 as the polarization of the system? or Is this a numerical error?<br><br>The material in the z-direction is not polar, is it possible that the polarization operator for a non-polar material is ill defined for this calculation? Should I use berry phase polarization for the polar materials only?<br><br>Also, when I change the dimension of the unit cell in the z-direction and position of the material inside the cell, this polarization value changes even though I use the same material.<br><br>Most importantly, this polarization value changes with the number of processors used with the mpirun -np. This is very strange, should I stick to -np 1. I should not make it parallel?</div><div><br></div><div>I also saw that this polarization is given with microC/cm^2 in the papers. Does this mean the last two digits from  " P =  -0.0000002  (mod   2.9386420)  C/m^2"  are not considered from the results obtained in the Quantum ESPRESSO run?</div><div><br>Are all those numerical errors or physical phenomena because I know that there are vacuum fluctuations and Lamb modes, and so on? Can we probe those with QE?<br><br>I look forward to hearing from you soon.<br><br>Best regards,<br>Aykut Turfanda<br>PhD student at Istanbul Technical University<br><br><br>###############################<br> SUMMARY OF PHASES<br>                               ~~~~~~~~~~~~~~~~~<br><br>                          Ionic Phase: -0.40381<br>                     Electronic Phase:  0.40381<br>                          TOTAL PHASE: -0.00000 MOD_TOT: 1<br><br><br>                             VALUES OF POLARIZATION<br>                             ~~~~~~~~~~~~~~~~~~~~~~<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><br>           P =  -0.0000030  (mod  40.7880378)  (e/Omega).bohr<br><br>           P =  -0.0000000  (mod   0.0514003)  e/bohr^2<br><br>           P =  -0.0000002  (mod   2.9386420)  C/m^2<br><br>        The polarization direction is:  ( 0.00000 , 0.00000 , 1.00000 )<br><br>################################<br>&control<br>calculation = 'nscf',<br>restart_mode='from_scratch'<br>prefix = 'bn'<br>outdir = './tempdir'<br>pseudo_dir = '~/pseudo_my/'<br>etot_conv_thr = 1d-05<br>forc_conv_thr =   1d-04<br>lberry=.true.<br>gdir          = 3<br>nppstr        = 10<br>/<br>&system<br>ibrav =  0,<br>nat =  2,<br>ntyp = 2,<br>ecutwfc = 80<br>occupations='fixed'<br>smearing='gaussian'<br>degauss=0.02<br>/<br>&electrons<br> conv_thr =   1d-08<br>/<br>&ions<br>/<br>&cell<br>/<br>ATOMIC_SPECIES<br>B    10.811  B.pbe.nc.sr.upf<br>N    14.0067  N.pbe.nc.sr.upf<br>ATOMIC_POSITIONS angstrom<br>B            -0.0000000000        1.4480800000        1.6085249900<br>N             1.2540740000        0.7240400000        1.6085249900<br>K_POINTS automatic<br>    10 10 10 0 0 0<br>CELL_PARAMETERS angstrom<br>2.5081477165 0.0000000000 0.0000000000<br>-1.2540738583 2.1721196390 0.0000000000<br>0.0000000000 0.0000000000 21.5841000986<br>##############################<br>&control<br>calculation = 'scf',<br>restart_mode='from_scratch'<br>prefix = 'bn'<br>outdir = './tempdir'<br>pseudo_dir = '~/pseudo_my/'<br>etot_conv_thr = 1d-05<br>forc_conv_thr =   1d-04<br>/<br>&system<br>ibrav =  0,<br>nat =  2,<br>ntyp = 2,<br>ecutwfc = 80<br>occupations='fixed'<br>smearing='gaussian'<br>degauss=0.02<br>/<br>&electrons<br> conv_thr =   1d-08<br>/<br>&ions<br>/<br>&cell<br>/<br>ATOMIC_SPECIES<br>B    10.811  B.pbe.nc.sr.upf<br>N    14.0067  N.pbe.nc.sr.upf<br>ATOMIC_POSITIONS angstrom<br>B            -0.0000000000        1.4480800000        1.6085249900<br>N             1.2540740000        0.7240400000        1.6085249900<br>K_POINTS automatic<br>    10 10 1 0 0 0<br>CELL_PARAMETERS angstrom<br>2.5081477165 0.0000000000 0.0000000000<br>-1.2540738583 2.1721196390 0.0000000000<br>0.0000000000 0.0000000000 21.5841000986<br>#################################<br></div></div></div>