<table cellspacing="0" cellpadding="0" border="0" ><tr><td valign="top" style="font: inherit;"><div id="yiv1021237495"><table id="yiv1021237495bodyDrftID" class="yiv1021237495" border="0" cellpadding="0" cellspacing="0"><tbody><tr><td id="yiv1021237495drftMsgContent" style="font-style: inherit; font-variant: inherit; font-weight: inherit; line-height: inherit; font-size-adjust: inherit; font-stretch: inherit; -x-system-font: none; font-family: arial; font-size: 10pt;">Dear Giuseppe:<br>I see my output file:<br>bfgs converged in&nbsp;&nbsp; 7 scf cycles and&nbsp;&nbsp; 4 bfgs steps<br>
&nbsp;&nbsp;&nbsp;&nbsp; (criteria: energy &lt; 0.10E-04, force &lt; 0.10E-03, cell &lt; 0.50E+00)<br><br>the con<br><br>total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp; -338.60645412 Ry<br>&nbsp;&nbsp;&nbsp;&nbsp; Harris-Foulkes estimate&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp; -338.60645412 Ry<br>&nbsp;&nbsp;&nbsp;&nbsp; estimated scf accuracy&nbsp;&nbsp;&nbsp; &lt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 3.3E-12 Ry<br><br>&nbsp;&nbsp;&nbsp;&nbsp; The total energy is the sum of the following terms:<br><br>&nbsp;&nbsp;&nbsp;&nbsp; one-electron contribution =&nbsp;&nbsp;&nbsp;&nbsp; -91.38378899 Ry<br>&nbsp;&nbsp;&nbsp;&nbsp; hartree contribution&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 81.54028728 Ry<br>&nbsp;&nbsp;&nbsp;&nbsp; xc contribution&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp; -113.31772437 Ry<br>&nbsp;&nbsp;&nbsp;&nbsp;
 ewald contribution&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp; -215.44522804 Ry<br><br>&nbsp;&nbsp;&nbsp;&nbsp; convergence has been achieved in&nbsp;&nbsp; 8 iterations<br><br>&nbsp;&nbsp;&nbsp;&nbsp; Forces acting on atoms (Ry/au):<br><br>&nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp;&nbsp; 1 type&nbsp; 1&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp;&nbsp; 0.00000000&nbsp;&nbsp;&nbsp; 0.00000000&nbsp;&nbsp;&nbsp; 0.00004539<br>&nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp;&nbsp; 2 type&nbsp; 1&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp;&nbsp; 0.00000000&nbsp;&nbsp;&nbsp; 0.00000000&nbsp;&nbsp;&nbsp; 0.00004539<br>&nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp;&nbsp; 3 type&nbsp; 2&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp;&nbsp; 0.00000000&nbsp;&nbsp;&nbsp; 0.00000000&nbsp;&nbsp; -0.00004539<br>&nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp;&nbsp; 4 type&nbsp; 2&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp;&nbsp; 0.00000000&nbsp;&nbsp;&nbsp; 0.00000000&nbsp;&nbsp;
 -0.00004539<br><br>&nbsp;&nbsp;&nbsp;&nbsp; Total force =&nbsp;&nbsp;&nbsp;&nbsp; 0.000091&nbsp;&nbsp;&nbsp;&nbsp; Total SCF correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000001<br><br><br>&nbsp;&nbsp;&nbsp;&nbsp; entering subroutine stress ...<br><br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; total&nbsp;&nbsp; stress&nbsp; (Ry/bohr**3)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; (kbar)&nbsp;&nbsp;&nbsp;&nbsp; P=&nbsp;&nbsp;&nbsp; 0.04<br>&nbsp; -0.00000057&nbsp;&nbsp; 0.00000000&nbsp;&nbsp; 0.00000000&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -0.08&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.00&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.00<br>&nbsp;&nbsp; 0.00000000&nbsp; -0.00000057&nbsp;&nbsp; 0.00000000&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.00&nbsp;&nbsp;&nbsp;&nbsp; -0.08&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.00<br>&nbsp;&nbsp; 0.00000000&nbsp;&nbsp; 0.00000000&nbsp;&nbsp;
 0.00000198&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.00&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.00&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.29<br><br><br>&nbsp;&nbsp;&nbsp;&nbsp; bfgs converged in&nbsp;&nbsp; 7 scf cycles and&nbsp;&nbsp; 4 bfgs steps<br>&nbsp;&nbsp;&nbsp;&nbsp; (criteria: energy &lt; 0.10E-04, force &lt; 0.10E-03, cell &lt; 0.50E+00)<br><br>&nbsp;&nbsp;&nbsp;&nbsp; End of BFGS Geometry Optimization<br><br>&nbsp;&nbsp;&nbsp;&nbsp; Final enthalpy =&nbsp;&nbsp;&nbsp; -338.6064541238 Ry<br>Begin final coordinates<br>&nbsp;&nbsp;&nbsp;&nbsp; new unit-cell volume =&nbsp;&nbsp;&nbsp; 331.68137 a.u.^3 (&nbsp;&nbsp;&nbsp; 49.15011 Ang^3 )<br><br>CELL_PARAMETERS (alat=&nbsp; 6.22300000)<br>&nbsp;&nbsp; 0.994801477&nbsp;&nbsp; 0.000000000&nbsp;&nbsp; 0.000000000<br>&nbsp; -0.497400738&nbsp;&nbsp; 0.861523350&nbsp;&nbsp; 0.000000000<br>&nbsp;&nbsp; 0.000000000&nbsp;&nbsp; 0.000000000&nbsp;&nbsp; 1.605898909<br><br>ATOMIC_POSITIONS
 (crystal)<br>Zn&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.333333333&nbsp;&nbsp; 0.666666667&nbsp; -0.000036499<br>Zn&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.666666667&nbsp;&nbsp; 0.333333333&nbsp;&nbsp; 0.499963501<br>O&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.333333333&nbsp;&nbsp; 0.666666667&nbsp;&nbsp; 0.379036499<br>O&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.666666667&nbsp;&nbsp; 0.333333333&nbsp;&nbsp; 0.879036499<br>End final coordinates<br><br>Best <br>Regards<br>Yun Song,Kang<br>
Department Physical Science and Technology of  Inner Mongolia University.<br><br>--- <b>11年10月14日，周五, Giuseppe Mattioli <i>&lt;giuseppe.mattioli@mlib.ism.cnr.it&gt;</i></b> 写道：<br><blockquote style="border-left: 2px solid rgb(16, 16, 255); margin-left: 5px; padding-left: 5px;"><br>发件人: Giuseppe Mattioli &lt;giuseppe.mattioli@mlib.ism.cnr.it&gt;<br>主题: Re: [Pw_forum] Zn NC_PP from qe-forge<br>收件人: "PWSCF Forum" &lt;pw_forum@pwscf.org&gt;<br>日期: 2011年10月14日,周五,下午9:12<br><br><div class="yiv1021237495plainMail"><br>Dear Yun Song<br>A 62Ry cutoff could be quite poor both in the case of your NC Zn (it contains <br>a quite deep 3d shell) and for your NC O.<br>You should reach concergence with respect to the wf cutoff before attempting <br>to perform phonon calculations.<br>HTH<br>Giuseppe<br><br><br>On Friday 14 October 2011 15:03:48 之王 改王 wrote:<br>&gt; Hi,<br>&gt; &nbsp; Recently,I generate Zn NC-PP
 using the PP
 library from qe-forge project.<br>&gt; I want to generate Zn.pz-nc.UPF using ld1.x.<br>&gt; My Zn.-nc.in is followed:<br>&gt; &nbsp;&amp;input<br>&gt; &nbsp;&nbsp; title='Zn',<br>&gt; &nbsp;&nbsp; zed=30.,<br>&gt; &nbsp;&nbsp; rel=0,<br>&gt; &nbsp;&nbsp; config='[Ar] 4s2.0 4p0.0 3d10.0',<br>&gt; &nbsp;&nbsp; iswitch=3,<br>&gt; &nbsp;&nbsp; dft='PZ'<br>&gt; &nbsp;/<br>&gt; &nbsp;&amp;inputp<br>&gt; &nbsp;&nbsp; pseudotype=2,<br>&gt; &nbsp;&nbsp; file_pseudopw='Zn.pz-nc.UPF',<br>&gt; &nbsp;&nbsp; author='TM',<br>&gt; &nbsp;&nbsp; lloc=0,<br>&gt; &nbsp;&nbsp; nlcc=.true.,<br>&gt; &nbsp;&nbsp; tm=.true.<br>&gt; &nbsp;/<br>&gt; 3<br>&gt; 4P&nbsp; 2&nbsp; 1&nbsp; 0.00&nbsp; 0.00&nbsp; 2.30&nbsp; 2.30&nbsp; 0.0<br>&gt; 3D&nbsp; 3&nbsp; 2 10.00&nbsp; 0.00&nbsp; 2.30&nbsp; 2.30&nbsp; 0.0<br>&gt; 4S&nbsp; 1&nbsp; 0&nbsp; 2.00&nbsp; 0.00&nbsp; 2.30&nbsp; 2.30&nbsp; 0.0<br>&gt;<br>&gt;<br>&gt; And I calculate ZnO vc-relax using Zn.pz-nc.UPF and O.pz-mt.UPF,<br>&gt;
 the theory is a=6.223 c/a=1.614, after calculating vc-relax is a=6.191<br>&gt; c/a=1.614 with ecutwfc=62Ry.<br>&gt;<br>&gt; However, No matter how I try to calculate,the phonon frequency is negative<br>&gt;<br>&gt; The scf and ph file is followed:<br>&gt; &amp;CONTROL<br>&gt; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; calculation = 'scf' ,<br>&gt; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; restart_mode = 'from_scratch' ,<br>&gt; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; outdir = '/root/ZnO/915-phonon' ,<br>&gt; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; pseudo_dir = '/root/pseudopotential' ,<br>&gt;
 &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; prefix = 'ZnO' ,<br>&gt; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; etot_conv_thr = 1.0D-5 ,<br>&gt; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; forc_conv_thr = 1.0D-4 ,<br>&gt; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; tstress = .true. ,<br>&gt; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; tprnfor = .true. ,<br>&gt; &nbsp;/<br>&gt; &nbsp;&amp;SYSTEM<br>&gt; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ibrav = 4,<br>&gt;
 &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; celldm(1) = 6.191,<br>&gt; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; celldm(3) = 1.614,<br>&gt; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; nat = 4,<br>&gt; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ntyp = 2,<br>&gt; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ecutwfc = 62 ,<br>&gt; &nbsp;/<br>&gt; &nbsp;&amp;ELECTRONS<br>&gt; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; conv_thr = 1.0d-10 ,<br>&gt;
 &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; mixing_mode = 'plain' ,<br>&gt; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; mixing_beta = 0.6 ,<br>&gt; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; diagonalization = 'david' ,<br>&gt; &nbsp;/<br>&gt; ATOMIC_SPECIES<br>&gt; &nbsp;&nbsp; Zn&nbsp;&nbsp; 65.40900&nbsp; Zn.pz-nc.UPF<br>&gt; &nbsp;&nbsp;&nbsp; O&nbsp;&nbsp; 15.99940&nbsp; O.pz-mt.UPF<br>&gt; ATOMIC_POSITIONS crystal<br>&gt; Zn&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.333333333&nbsp;&nbsp; 0.666666667&nbsp; -0.000036499&nbsp;&nbsp; 1 1 1<br>&gt; Zn&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.666666667&nbsp;&nbsp; 0.333333333&nbsp;&nbsp; 0.499963501&nbsp;&nbsp; 1 1 1<br>&gt; O&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.333333333&nbsp;&nbsp; 0.666666667&nbsp;&nbsp; 0.379036499&nbsp;&nbsp; 1 1 1<br>&gt;
 O&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.666666667&nbsp;&nbsp; 0.333333333&nbsp;&nbsp; 0.879036499&nbsp;&nbsp; 1 1 1<br>&gt; K_POINTS automatic<br>&gt; &nbsp; 4 4 4&nbsp;&nbsp; 1 1 1<br>&gt;<br>&gt;<br>&gt;<br>&gt; phonons of ZnO&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;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<br>&gt; &nbsp;&amp;INPUTPH<br>&gt; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; outdir = '/root/ZnO/915-phonon' ,<br>&gt; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; prefix = 'ZnO' ,<br>&gt; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; fildyn =
 'ZnO_IR.dyn' ,<br>&gt; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; fildvscf = 'ZnO.dv' ,<br>&gt; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ldisp = .true.,<br>&gt; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; nq1 = 4 ,<br>&gt; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; nq2 = 4 ,<br>&gt; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; nq3 = 4 ,<br>&gt; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; epsil =
 .false.,<br>&gt; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; elph = .false.,<br>&gt; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; fpol = .false. ,<br>&gt; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; recover = .false. ,<br>&gt; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; amass(1) = 65.4090,<br>&gt; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; amass(2) = 15.9994,<br>&gt; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; tr2_ph = 1.0d-12 ,<br>&gt; &nbsp;/<br>&gt;<br>&gt;
 Best<br>&gt; Regards<br>&gt; Yun Song,Kang<br>&gt;<br>&gt; Department Physical Science and Technology of&nbsp; Inner Mongolia University.<br><br><br><br>-- <br>********************************************************<br>- Article premier - Les hommes naissent et demeurent<br>libres et ègaux en droits. Les distinctions sociales<br>ne peuvent être fondèes que sur l'utilitè commune<br>- Article 2 - Le but de toute association politique<br>est la conservation des droits naturels et <br>imprescriptibles de l'homme. Ces droits sont la libertè,<br>la propriètè, la sùretè et la rèsistance à l'oppression.<br>********************************************************<br><br>&nbsp; &nbsp;Giuseppe Mattioli &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;<br>&nbsp; &nbsp;CNR - ISTITUTO DI STRUTTURA DELLA MATERIA &nbsp; <br>&nbsp; &nbsp;v. Salaria Km 29,300 - C.P. 10 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;
 &nbsp; &nbsp;<br>&nbsp; &nbsp;I 00015 - Monterotondo Stazione (RM) &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;<br>&nbsp; &nbsp;Tel + 39 06 90672836 - Fax +39 06 90672316 &nbsp; &nbsp;<br>&nbsp; &nbsp;E-mail: &lt;<a rel="nofollow">giuseppe.mattioli@ism.cnr.it</a>&gt;<br>_______________________________________________<br>Pw_forum mailing list<br><a rel="nofollow">Pw_forum@pwscf.org</a><br><a rel="nofollow" target="_blank" href="http://www.democritos.it/mailman/listinfo/pw_forum">http://www.democritos.it/mailman/listinfo/pw_forum</a><br></div></blockquote></td></tr></tbody></table></div></td></tr></table>