Hi ,<div><br><div>I reproduced your PRL paper on NbN monolayers https://journals.aps.org/prl/abstract/10.1103/PhysRevLett.123.037601.</div><div><br></div><div>I am following the steps as described by Quantum espresso example file on PbTiO3. But with this method, I am getting really weird numbers for polarization (really large numbers are coming). </div><div>Also QE write the results C/m^2. That means it is dividing the dipole moment by volume.</div><div>Is QE using the volume of the unit cell which includes vacuum also ?</div><div><br></div><div>Best regards</div><div>Djamel </div><div><br></div><div> </div><div><br></div><div>Could you please help me in figuring out where am I going wrong? I attached my scf input file and the non-scf input and output files to calculate the polarization by implementing berry phase rule. The files b1, and b2 corresponds to the direction of polarization along x and z-direction.</div><div><br><div id="ymail_android_signature"><a id="ymail_android_signature_link" href="https://mail.onelink.me/107872968?pid=nativeplacement&c=Global_Acquisition_YMktg_315_Internal_EmailSignature&af_sub1=Acquisition&af_sub2=Global_YMktg&af_sub3=&af_sub4=100000604&af_sub5=EmailSignature__Static_">Envoyé depuis Yahoo Mail pour Android</a></div></div></div>