<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><br><div><div>Il giorno 29/ago/2010, alle ore 16.18, reza shidpoor ha scritto:</div><br><blockquote type="cite"><b>2 ) The RAM of my system is sufficient for 13 atoms (RAM = 8 Gigabyte DDR3).<br></b></blockquote><div><br></div><div>How can you be sure? It depends on a lot of factors, besides the number of atoms: cutoffs, number of calculated bands, cell size, diagonalization method, ..., and also pseudopotential type.</div><div><br></div><div>Indeed, US PPs have usually many more projector functions than NC PPs.</div><br><blockquote type="cite"><b>3 ) I tested a case . I used informations related to 01 example ( Si - scf ).<br> I constructed a input file from a imaginary 13-atom Si cluster that is different from<br>
previous input file in type of atoms and pseudopotential. Then I run scf calculation both of them.</b><br><br><i>Input File related to Au :<br><br> &control<br> calculation='scf'<br> title='Au13Cluster'<br>
restart_mode='from_scratch'<br> outdir='./tmp'<br> prefix='Au13'<br> etot_conv_thr=1.0D-4<br> pseudo_dir = './'<br> /<br> &system<br> ibrav= 2, celldm(1) = 38,nat= 13, ntyp= 1,<br>
ecutwfc =40.0, ecutrho=400.0,<br> occupations='smearing',smearing='gaussian',degauss=0.01,<br> nspin=1,<br> /<br> &electrons<br> electron_maxstep=100,<br> diagonalization='cg',<br>
conv_thr=1.0D-6,<br> mixing_beta=0.5,<br> /<br>&ions<br> ion_dynamics='bfgs'<br>/<br>&cell<br>/<br><br>ATOMIC_SPECIES<br> Au 196.97 Au.pbe-nd-van.UPF<br>ATOMIC_POSITIONS {angstrom}<br>Au 4.844 4.844 4.844<br>
Au 0.000 4.844 7.833<br>Au 4.844 7.833 0.000<br>Au 4.844 7.833 9.688<br>Au 9.688 4.844 1.855<br>Au 1.855 0.000 4.844<br>Au 1.855 9.688 4.844<br>Au 4.844 1.855 9.688<br>
Au 7.833 9.688 4.844<br>Au 4.844 1.855 0.000<br>Au 9.688 4.844 7.833<br>Au 7.833 0.000 4.844<br>Au 0.000 4.844 1.855<br>K_POINTS {automatic}<br>2 2 2 1 1 1 </i><br><br>
<br><b>After of several minutes, the result was </b><br><br><br> <i> Largest allocated arrays est. size (Mb) dimensions<br> Kohn-Sham Wavefunctions 110.74 Mb ( 84386, 86)<br> NL pseudopotentials 217.61 Mb ( 84386, 169)<br>
Each V/rho on FFT grid 51.50 Mb (3375000)<br> Each G-vector array 9.18 Mb (1203541)<br> G-vector shells 0.03 Mb ( 3661)<br> Largest temporary arrays est. size (Mb) dimensions<br>
Each subspace H/S matrix 0.11 Mb ( 86, 86)<br> Each <psi_i|beta_j> matrix 0.22 Mb ( 169, 86)<br> Arrays for rho mixing 411.99 Mb (3375000, 8)<br><br> Initial potential from superposition of free atoms<br>
Check: negative starting charge= -0.356687<br><br> starting charge 142.84904, renormalised to 143.00000<br><br> negative rho (up, down): 0.357E+00 0.000E+00<br> Starting wfc are 117 atomic wfcs<br><br>
total cpu time spent up to now is 255.53 secs<br><br> per-process dynamical memory: 1035.0 Mb<br></i></blockquote><div><br></div><div>This estimate is usually a lower bound. I found that the memory requirement at some points of the run can as big as the double of that number. On how many processors was this run done?</div><br><blockquote type="cite"><i><br> Self-consistent Calculation<br><br> iteration # 1 ecut= 51.00 Ry beta=0.50<br> CG style diagonalization<br>
ethr = 1.00E-02, avg # of iterations = 2.1<br>forrtl: severe (174): SIGSEGV, segmentation fault occurred<br>Image PC Routine Line Source <br>pw.x 000000000052089C addusdens_ 27 addusdens.f90<br>
pw.x 00000000004DCB76 sum_band_ 153 sum_band.f90<br>pw.x 000000000044B7D7 electrons_ 287 electrons.f90<br>pw.x 000000000040BA6E MAIN__ 92 pwscf.f90<br>
pw.x 000000000040B8DC Unknown Unknown Unknown<br>libc.so.6 00002B500B54CABD Unknown Unknown Unknown<br>pw.x 000000000040B7D9 Unknown Unknown Unknown</i><br></blockquote><div><br></div><div>You may obtain this kind of messages when you hit the maximum available RAM memory or when the code tries to access a memory address not reserved to it. To exclude the latter you could try to recompile with -g -check bounds -traceback and rerun. To exclude the former, try a smaller system or run on more processors on a bigger machine.</div><br><blockquote type="cite">
<br><br> <b>Now , for input file related Si </b><br><br><i>&control<br> calculation='scf<br> title='Si13Cluster<br> restart_mode='from_scratch<br> outdir='./tmp<br> prefix='Si13<br> etot_conv_thr=1.0D-4<br>
pseudo_dir = './<br> /<br> &system<br> ibrav= 2, celldm(1) = 38,nat= 13, ntyp= 1, <br> ecutwfc =40.0, ecutrho=400.0, <br> occupations='smearing',smearing='mp',degauss= 0.01, <br> nspin=1,<br>
/<br> &electrons<br> electron_maxstep=100,<br> diagonalization='cg', <br> conv_thr=1.0D-6,<br> mixing_beta=0.5,<br> /<br>&ions<br> <br>/<br>&cell<br>/<br><br>ATOMIC_SPECIES<br> Si 28.086 Si.pz-vbc.UPF<br>
ATOMIC_POSITIONS {angstrom}<br>Si 0.000 0.000 0.000<br>Si 0.000 4.844 7.833<br>Si 4.844 7.833 0.000<br>Si 4.844 7.833 9.688<br>Si 9.688 4.844 1.855<br>Si 1.855 0.000 4.844<br>
Si 1.855 9.688 4.844<br>Si 4.844 1.855 9.688<br>Si 7.833 9.688 4.844<br>Si 4.844 1.855 0.000<br>Si 9.688 4.844 7.833<br>Si 7.833 0.000 4.844<br>Si 0.000 4.844 1.855<br>
K_POINTS {automatic}<br>2 2 2 1 1 1<br></i><br><b>And the result was :</b><br><br><i>...............<br>...............<br>iteration # 50 ecut= 40.00 Ry beta=0.50<br> CG style diagonalization<br> ethr = 2.38E-05, avg # of iterations = 3.8<br>
<br> negative rho (up, down): 0.442E-07 0.000E+00<br><br> total cpu time spent up to now is 7717.21 secs<br><br> total energy = -97.69548088 Ry<br> Harris-Foulkes estimate = -97.70290819 Ry<br>
estimated scf accuracy < 0.46657030 Ry<br><br> iteration # 51 ecut= 40.00 Ry beta=0.50<br> CG style diagonalization<br> ethr = 2.38E-05, avg # of iterations = 3.0<br><br> negative rho (up, down): 0.250E-07 0.000E+00<br>
<br> total cpu time spent up to now is 7843.25 secs<br><br> total energy = -97.69672289 Ry<br> Harris-Foulkes estimate = -97.69637710 Ry<br> estimated scf accuracy < 0.33411951 Ry<br>
...................<br>...................<br>...................</i><br><b>and continues : The simulation was done simply !!!!!!</b><br></blockquote><div><br></div><div>It didn't converge after 51 iterations, this may be suspicious...</div><br><blockquote type="cite"><br><br><b>Dear Gabriele<br><br>I have almost waste one month because of this problem and I diddnot solve problem.<br></b></blockquote><div><br></div><div>Perhaps you are not approaching it in the correct way.</div><div><br></div><blockquote type="cite"><b>
<br>With regard to 1,2,3 points , Are you agree with me that the source of problem is only and only pseudopotential ?<br><br>Best wishes.<br></b></blockquote><div><br></div><div>Regards,</div><div><br></div><div>Gabriele</div><div><br></div><div><br></div><blockquote type="cite"><b><br>Reza.Shidpour<br>Institute for Nanoscienec and Nanotechnology,<br>Sharif University of Technology (SUT)<br>
</b><br>
_______________________________________________<br>Pw_forum mailing list<br><a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>http://www.democritos.it/mailman/listinfo/pw_forum<br></blockquote></div><br><div>
<span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><div><span class="Apple-style-span" style="color: rgb(126, 126, 126); font-size: 16px; font-style: italic; "><br class="Apple-interchange-newline">§ Gabriele Sclauzero, EPFL SB ITP CSEA</span></div><div><font class="Apple-style-span" color="#7E7E7E"><i> PH H2 462, Station 3, CH-1015 Lausanne</i></font></div></span>
</div>
<br></body></html>