<div class="__aliyun_email_body_block"><div  style="font-family:Tahoma,Arial,STHeitiSC-Light,SimSun;"><div  style="font-family:Tahoma,Arial,STHeitiSC-Light,SimSun;clear:both;"><span >Dear All,<br ></span><div  style="clear:both;"><br ></div><div  style="clear:both;">I am doing the dipole correction to my slab system, where the dipole is added along the z axis.<br ></div><div  style="clear:both;"><br ></div><div  style="clear:both;">This is my input file in the SCF calculation step.<br ></div><div  style="clear:both;">&CONTROL<br ></div><div  style="clear:both;">   title = 'Graphene_SiO2',<br ></div><div  style="clear:both;">   calculation = 'scf',<br ></div><div  style="clear:both;">   restart_mode = 'from_scratch',<br ></div><div  style="clear:both;">   tefield = .true.,<br ></div><div  style="clear:both;">   dipfield = .true.,<br ></div><div  style="clear:both;">   pseudo_dir = './'<br ></div><div  style="clear:both;">/<br ></div><div  style="clear:both;">&SYSTEM<br ></div><div  style="clear:both;">   ibrav = 0,<br ></div><div  style="clear:both;">   nat = 2,<br ></div><div  style="clear:both;">   ntyp = 1,<br ></div><div  style="clear:both;">   ecutwfc = 50,<br ></div><div  style="clear:both;">   ecutrho = 200,<br ></div><div  style="clear:both;">   nosym = .true.,<br ></div><div  style="clear:both;">   noinv = .true.,<br ></div><div  style="clear:both;">   noncolin =.true.,<br ></div><div  style="clear:both;">   lspinorb=.true.,<br ></div><div  style="clear:both;">   edir = 3,<br ></div><div  style="clear:both;">   eopreg = 0.2,<br ></div><div  style="clear:both;">   emaxpos = 0.8,<br ></div><div  style="clear:both;">   eamp = 0.001,<br ></div><div  style="clear:both;">   nbnd = 50,<br ></div><div  style="clear:both;">   occupations = 'smearing',<br ></div><div  style="clear:both;">   degauss = 1d-9,<br ></div><div  style="clear:both;">   vdw_corr = DFT-D3,<br ></div><div  style="clear:both;">   smearing = 'gaussian'<br ></div><div  style="clear:both;">/<br ></div><div  style="clear:both;">&ELECTRONS<br ></div><div  style="clear:both;">   electron_maxstep = 100,<br ></div><div  style="clear:both;">   conv_thr = 1.0d-9,<br ></div><div  style="clear:both;">   mixing_mode = 'plain',<br ></div><div  style="clear:both;">   mixing_beta = 0.8d0,<br ></div><div  style="clear:both;">   diagonalization = 'david'<br ></div><div  style="clear:both;">/<br ></div><div  style="clear:both;">&IONS<br ></div><div  style="clear:both;">    ion_positions = 'from_input'<br ></div><div  style="clear:both;">/<br ></div><div  style="clear:both;">&CELL<br ></div><div  style="clear:both;">    cell_dynamics = 'none'<br ></div><div  style="clear:both;">/<br ></div><div  style="clear:both;">ATOMIC_SPECIES<br ></div><div  style="clear:both;">  Si  14.0 Si_ONCV_PBE_fr.upf<br ></div><div  style="clear:both;">   O   8.0 O_ONCV_PBE_fr.upf<br ></div><div  style="clear:both;">   C  12.0 C_ONCV_PBE_fr.upf<br ></div><div  style="clear:both;">CELL_PARAMETERS (angstrom)<br ></div><div  style="clear:both;">  4.9130001068 0.0000000000  0.0000000000<br ></div><div  style="clear:both;"> -2.4565000534 4.2547829013  0.0000000000<br ></div><div  style="clear:both;">  0.0000000000 0.0000000000 27.7607002258<br ></div><div  style="clear:both;">ATOMIC_POSITIONS (crystal)<br ></div><div  style="clear:both;">Si  0.52990000      0.52990000     0.48013653<br ></div><div  style="clear:both;">Si  0.00000000      0.47010000     0.41523428<br ></div><div  style="clear:both;">Si  0.47010000      0.00000000     0.35033203<br ></div><div  style="clear:both;">O   0.85400000      0.58640000     0.45692101<br ></div><div  style="clear:both;">O   0.26760000      0.41360000     0.43844982<br ></div><div  style="clear:both;">O   0.14600000      0.73240000     0.37354756<br ></div><div  style="clear:both;">O   0.73240000      0.14600000     0.39201875<br ></div><div  style="clear:both;">O   0.58640000      0.85400000     0.30864531<br ></div><div  style="clear:both;">O   0.41360000      0.26760000     0.32711649<br ></div><div  style="clear:both;">C   0.52990000      0.02990000     0.58820290<br ></div><div  style="clear:both;">C   0.69656667      0.36323333     0.58820290<br ></div><div  style="clear:both;">C   0.52990000      0.52990000     0.58820290<br ></div><div  style="clear:both;">C   0.69656667      0.86323333     0.58820290<br ></div><div  style="clear:both;">C   0.02990000      0.02990000     0.58820290<br ></div><div  style="clear:both;">C   0.19656667      0.36323333     0.58820290<br ></div><div  style="clear:both;">C   0.02990000      0.52990000     0.58820290<br ></div><div  style="clear:both;">C   0.19656667      0.86323333     0.58820290<br ></div><div  style="clear:both;">K_POINTS automatic<br ></div><div  style="clear:both;">15 15 1 1 1 1<br ></div><div  style="clear:both;"><br ></div><div  style="clear:both;">This is my input file in the pp.x calculation step.<br ></div><div  style="clear:both;">&INPUTPP<br ></div><div  style="clear:both;">prefix = 'pwscf',<br ></div><div  style="clear:both;">outdir = './',<br ></div><div  style="clear:both;">plot_num = 11,<br ></div><div  style="clear:both;">/<br ></div><div  style="clear:both;">&PLOT<br ></div><div  style="clear:both;">fileout = 'dipole.dat',<br ></div><div  style="clear:both;">/<br ></div><div  style="clear:both;"><br ></div><div  style="clear:both;">This is my input file in the average.x calculation step.<br ></div><div  style="clear:both;">1<br ></div><div  style="clear:both;">tmp.pp<br ></div><div  style="clear:both;">1.0<br ></div><div  style="clear:both;">3000<br ></div><div  style="clear:both;">3<br ></div><div  style="clear:both;">3.0<br ></div><div  style="clear:both;"><br ></div><div  style="clear:both;">Then I obtain the avg.dat file. When I plot the dipole, I see no correction to the system at all.<br ></div><div  style="clear:both;"><img  height="737" src="cid:__aliyun172906889518529251" width="1179" style="vertical-align:bottom;"></div><div  style="clear:both;">Would anyone please give me some suggestions on how to apply the correct the dipole to the system?<br ></div><div  style="clear:both;"><br ></div><span >Thank you in advance.</span></div></div></div>