<div dir="ltr">Dear pwscf developers and users,<br><br>I am trying to calculate complex band structure of silicon using espresso-4.0.2. The output returns with the<br>error <br><br>&nbsp;%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%<br>
&nbsp;&nbsp;&nbsp;&nbsp; from init_orbitals : error #&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1<br>&nbsp;&nbsp;&nbsp;&nbsp; Too big atomic spheres<br>&nbsp;%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%<br><br>&nbsp;&nbsp;&nbsp;&nbsp; stopping ...<br><br>Below is the input for scf calculation as in example 05:<br>
<br>&nbsp;&amp;control<br>&nbsp;&nbsp;&nbsp; calculation=&#39;scf&#39;<br>&nbsp;&nbsp;&nbsp; restart_mode=&#39;from_scratch&#39;,<br>&nbsp;&nbsp;&nbsp; prefix=&#39;si&#39;<br>&nbsp;&nbsp;&nbsp; pseudo_dir = &#39;/home/rajan/espresso-4.0.2_P/pseudo/&#39;,<br>&nbsp;&nbsp;&nbsp; outdir=&#39;/home/rajan/tmp/&#39;<br>
&nbsp;/<br>&nbsp;&amp;system<br>&nbsp;&nbsp;&nbsp; ibrav= 2, celldm(1)= 10.2, nat= 2, ntyp= 1,<br>&nbsp;&nbsp;&nbsp; ecutwfc =18.0<br>&nbsp;/<br>&nbsp;&amp;electrons<br>&nbsp;&nbsp;&nbsp; conv_thr =&nbsp; 1.0d-8<br>&nbsp;&nbsp;&nbsp; mixing_beta = 0.7<br>&nbsp;/<br>ATOMIC_SPECIES<br>&nbsp;Si&nbsp; 28.086&nbsp; Si.vbc.UPF<br>
ATOMIC_POSITIONS<br>&nbsp;Si 0.00 0.00 0.00<br>&nbsp;Si 0.25 0.25 0.25<br>K_POINTS<br>&nbsp; 10<br>&nbsp;&nbsp; 0.1250000&nbsp; 0.1250000&nbsp; 0.1250000&nbsp;&nbsp; 1.00<br>&nbsp;&nbsp; 0.1250000&nbsp; 0.1250000&nbsp; 0.3750000&nbsp;&nbsp; 3.00<br>&nbsp;&nbsp; 0.1250000&nbsp; 0.1250000&nbsp; 0.6250000&nbsp;&nbsp; 3.00<br>&nbsp;&nbsp; 0.1250000&nbsp; 0.1250000&nbsp; 0.8750000&nbsp;&nbsp; 3.00<br>
&nbsp;&nbsp; 0.1250000&nbsp; 0.3750000&nbsp; 0.3750000&nbsp;&nbsp; 3.00<br>&nbsp;&nbsp; 0.1250000&nbsp; 0.3750000&nbsp; 0.6250000&nbsp;&nbsp; 6.00<br>&nbsp;&nbsp; 0.1250000&nbsp; 0.3750000&nbsp; 0.8750000&nbsp;&nbsp; 6.00<br>&nbsp;&nbsp; 0.1250000&nbsp; 0.6250000&nbsp; 0.6250000&nbsp;&nbsp; 3.00<br>&nbsp;&nbsp; 0.3750000&nbsp; 0.3750000&nbsp; 0.3750000&nbsp;&nbsp; 1.00<br>
&nbsp;&nbsp; 0.3750000&nbsp; 0.3750000&nbsp; 0.6250000&nbsp;&nbsp; 3.00<br><br>and the input to calculate the complex band structure is<br><br>&nbsp;&nbsp; &amp;inputcond<br>&nbsp;&nbsp;&nbsp;&nbsp; outdir = &#39;/home/rajan/tmp/&#39;<br>&nbsp;&nbsp;&nbsp;&nbsp; prefixl = &#39;si&#39;<br>&nbsp;&nbsp;&nbsp;&nbsp; band_file = &#39;<a href="http://bands.si">bands.si</a>&#39;<br>
&nbsp;&nbsp;&nbsp;&nbsp; ikind = 0<br>&nbsp;&nbsp;&nbsp;&nbsp; energy0 = 10.d0<br>&nbsp;&nbsp;&nbsp;&nbsp; denergy = -0.4d0<br>&nbsp;&nbsp;&nbsp;&nbsp; ewind = 3.d0<br>&nbsp;&nbsp;&nbsp;&nbsp; epsproj = 1.d-7<br>&nbsp;&nbsp;&nbsp;&nbsp; delgep = 1.d-10<br>&nbsp;&nbsp;&nbsp;&nbsp; cutplot = 3.d0<br>&nbsp; /<br>&nbsp;&nbsp;&nbsp;&nbsp; 1<br>&nbsp;&nbsp;&nbsp;&nbsp; 0.0 0.0 1.0<br>&nbsp;&nbsp;&nbsp;&nbsp; 60&nbsp;&nbsp;&nbsp;&nbsp; <br><br>The output returns with the following error:<br>
<br>&nbsp;&nbsp;&nbsp;&nbsp; Program POST-PROC v.4.0.2&nbsp; starts ...<br>&nbsp;&nbsp;&nbsp;&nbsp; Today is 14Sep2008 at&nbsp; 5:20:33<br><br>&nbsp;&nbsp;&nbsp;&nbsp; Parallel version (MPI)<br><br>&nbsp;&nbsp;&nbsp;&nbsp; Number of processors in use:&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1<br><br>&nbsp;&nbsp;&nbsp;&nbsp; Planes per process (thick) : nr3 = 20 npp =&nbsp; 20 ncplane =&nbsp; 400<br>
<br>&nbsp;&nbsp;&nbsp;&nbsp; Proc/&nbsp; planes cols&nbsp;&nbsp;&nbsp;&nbsp; G&nbsp;&nbsp;&nbsp; planes cols&nbsp;&nbsp;&nbsp; G&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; columns&nbsp; G<br>&nbsp;&nbsp;&nbsp;&nbsp; Pool&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; (dense grid)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; (smooth grid)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; (wavefct grid)<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp;&nbsp;&nbsp; 20&nbsp;&nbsp;&nbsp; 253&nbsp;&nbsp;&nbsp;&nbsp; 2733&nbsp;&nbsp; 20&nbsp;&nbsp;&nbsp; 253&nbsp;&nbsp;&nbsp;&nbsp; 2733&nbsp;&nbsp;&nbsp;&nbsp; 91&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 609<br><br>
<br>===== INPUT FILE containing the left lead =====<br><br>&nbsp;&nbsp;&nbsp;&nbsp; GEOMETRY:<br><br>&nbsp;&nbsp;&nbsp;&nbsp; lattice parameter (a_0)&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 10.2000&nbsp; a.u.<br>&nbsp;&nbsp;&nbsp;&nbsp; the volume&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp; 265.3020 (a.u.)^3<br>&nbsp;&nbsp;&nbsp;&nbsp; the cross section&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 26.0100 (a.u.)^2<br>
&nbsp;&nbsp;&nbsp;&nbsp; l of the unit cell&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.0000 (a_0)<br>&nbsp;&nbsp;&nbsp;&nbsp; number of atoms/cell&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 2<br>&nbsp;&nbsp;&nbsp;&nbsp; number of atomic types&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1<br><br>&nbsp;&nbsp;&nbsp;&nbsp; crystal axes: (cart. coord. in units of a_0)<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; a(1) = ( -0.5000&nbsp; 0.0000&nbsp; 0.5000 )<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; a(2) = (&nbsp; 0.0000&nbsp; 0.5000&nbsp; 0.5000 )<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; a(3) = ( -0.5000&nbsp; 0.5000&nbsp; 0.0000 )<br><br>&nbsp;&nbsp; Cartesian axes<br><br>&nbsp;&nbsp;&nbsp;&nbsp; site n.&nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; positions (a_0 units)<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Si&nbsp; tau(&nbsp; 1)=(&nbsp; 0.0000&nbsp; 0.0000&nbsp; 0.0000&nbsp; )<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Si&nbsp; tau(&nbsp; 2)=(&nbsp; 0.2500&nbsp; 0.2500&nbsp; 0.2500&nbsp; )<br><br>&nbsp;&nbsp;&nbsp;&nbsp; nr1s&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 20<br>&nbsp;&nbsp;&nbsp;&nbsp; nr2s&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 20<br>&nbsp;&nbsp;&nbsp;&nbsp; nr3s&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 20<br>&nbsp;&nbsp;&nbsp;&nbsp; nrx1s&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 20<br>
&nbsp;&nbsp;&nbsp;&nbsp; nrx2s&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 20<br>&nbsp;&nbsp;&nbsp;&nbsp; nrx3s&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 20<br>&nbsp;&nbsp;&nbsp;&nbsp; nr1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 20<br>&nbsp;&nbsp;&nbsp;&nbsp; nr2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 20<br>&nbsp;&nbsp;&nbsp;&nbsp; nr3&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 20<br>
&nbsp;&nbsp;&nbsp;&nbsp; nrx1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 20<br>&nbsp;&nbsp;&nbsp;&nbsp; nrx2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 20<br>&nbsp;&nbsp;&nbsp;&nbsp; nrx3&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 20<br><br>&nbsp;_______________________________<br>&nbsp; Radii of nonlocal spheres:<br><br>
&nbsp;&nbsp;&nbsp;&nbsp; type&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ibeta&nbsp;&nbsp;&nbsp;&nbsp; ang. mom.&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; radius (a_0 units)<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Si&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.2832<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Si&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.2978<br><br>&nbsp;%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%<br>
&nbsp;&nbsp;&nbsp;&nbsp; from init_orbitals : error #&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1<br>&nbsp;&nbsp;&nbsp;&nbsp; Too big atomic spheres<br>&nbsp;%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%<br><br>&nbsp;&nbsp;&nbsp;&nbsp; stopping ...<br><br><br>I appreciate any help/suggestion.<br>
<br>Thanks,<br><br>Rajan <br><br><br><br>&nbsp;<br></div>