<div dir="ltr"><div>Dear Users,</div><div><div><div><p>I am attempting 
to calculate the piezoelectric coefficients of 2D MoS2 using the Berry 
Phase method on a strained structure (0.5% in the armchair direction). </p><p>I faced two issues during the process.</p></div></div></div><div>1) In the nscf calculations, I added the 'nbnd' tag in the 
input file (<span style="color:rgb(0,0,0)"><a href="http://nscf1.in" target="_blank">nscf1.in</a></span>). However, when examining the output file 
(nscf1.out), the polarization value was displayed as NaN (see below). I have attached 
the output file for your reference. <br></div><div><br></div><div>                            VALUES OF POLARIZATION<br>                             ~~~~~~~~~~~~~~~~~~~~~~<br><br>        The calculation of phases done along the direction of vector 2<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 =         NaN  (mod  10.4913000)  (e/Omega).bohr<br><br>           P =         NaN  (mod   0.0037363)  e/bohr^2<br><br>           P =         NaN  (mod   0.2136113)  C/m^2<br><br>        The polarization direction is:  (-0.00000 , 1.00000 , 0.00000 )<br></div><div><br></div><div>However, after removing the 'nbnd' tag in the nscf input file (<a href="http://nscf.in" target="_blank">nscf.in</a>), I successfully 
obtained the polarization values (see below), and also attached the output file (nscf.out)). <br></div><div><br></div><div>                             VALUES OF POLARIZATION<br>                             ~~~~~~~~~~~~~~~~~~~~~~<br><br>        The calculation of phases done along the direction of vector 2<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 =  14.0341794  (mod  10.4913000)  (e/Omega).bohr<br><br>           P =   0.0049981  (mod   0.0037363)  e/bohr^2<br><br>           P =   0.2857471  (mod   0.2136113)  C/m^2<br><br>        The polarization direction is:  (-0.00000 , 1.00000 , 0.00000 )<br></div><div><br></div><div><br></div><div>2) At a strain of 0.5% in the armchair direction, the calculated 
polarization vector is P = 0.2857471 C/m². To convert it to C/m, I 
multiplied it by the armchair lattice parameter (10.5 angstroms) and 
obtained the piezoelectric coefficient e22 = 60006.89 pC/m. However, the
 published value for e22 is 362 pC/m (10.1021/acsnano.5b03394).</div><div><br></div><div>I have scavenged through the mailing list but found no tutorial on this topic. <br></div><br><div><span class="gmail_signature_prefix">-- </span></div><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><font><font>Best Regards</font></font></div><font color="#888888"><div><span style="color:rgb(68,68,68)"><span style="background-color:rgb(255,255,255)"><font><font>Zimmi Singh</font></font></span></span></div><div><span style="color:rgb(68,68,68)"><span style="background-color:rgb(255,255,255)"><font><font><i>Ph.D. Candidate</i></font></font></span></span></div></font></div></div><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><span style="color:rgb(68,68,68)"><i>Department of Metallurgical and Materials Engineering<br>
Indian Institute of Technology, Kharagpur<br>
Kharagpur, India</i></span></div><div><span style="color:rgb(68,68,68)"><i>Mob No. 9935804727</i></span><br></div><div><font><font><br></font></font></div></div></div></div>