<DIV>
<DIV>Dear all,</DIV>
<DIV>&nbsp;&nbsp;&nbsp; Recently, I&nbsp;have been&nbsp;performing calulation of frequency and zero point energy of CO2 molecule. However, unfortunately, I have obtained&nbsp;negative frequency. According to the manual, there are several possible reseaons: <BR>&nbsp;&nbsp;&nbsp; (1) wrong data file read.<BR>&nbsp;&nbsp;&nbsp; (2) wrong atomic masses given in input will yield wrong frequencies (but<BR>the content of file fildyn should be valid, since the force constants, not<BR>the dynamical matrix, are written to file).<BR>&nbsp;&nbsp;&nbsp; (3) convergence threshold for either SCF (conv thr) or phonon calculation<BR>(tr2 ph) too large: try to reduce them.<BR>&nbsp;&nbsp;&nbsp; (4) maybe your system does have negative or strange phonon frequencies,<BR>with the approximations you used. A negative frequency signals a<BR>mechanical instability of the chosen structure. Check that the structure<BR>is reasonable, and check the following parameters:<BR>&nbsp;&nbsp;&nbsp; – The cutoff for wavefunctions, ecutwfc<BR>&nbsp;&nbsp;&nbsp; – For US PP: the cutoff for the charge density, ecutrho<BR>&nbsp;&nbsp;&nbsp; – The k-point grid, especially for metallic systems!<BR></DIV>
<DIV>My calculated process as follows: </DIV>
<DIV>(1) Optimizing CO2 structure(input file as follows);</DIV>
<DIV>&nbsp;&amp;CONTROL<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; calculation = 'relax' ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; restart_mode = 'from_scratch' ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; outdir = '/home/olh/tmp/' ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; pseudo_dir = '/home/olh/pseudo/' ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; prefix = 'CO2' ,<BR>&nbsp;/<BR>&nbsp;&amp;SYSTEM<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ibrav = 1,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; celldm(1) = 15.11781502,<BR>&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 = 3,<BR>&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>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ecutwfc = 26 ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ecutrho = 260 ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; occupations = 'smearing' ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; degauss = 0.02 ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; smearing = 'methfessel-paxton' ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; nspin = 2 ,<BR>&nbsp;&nbsp; starting_magnetization(1) = 0.5,<BR>&nbsp;&nbsp; starting_magnetization(2) = 0.5,<BR>&nbsp;/<BR>&nbsp;&amp;ELECTRONS<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; conv_thr = 1.D-12 ,</DIV>
<DIV>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; mixing_beta = 0.4D0 ,<BR>&nbsp;/<BR>&nbsp;&amp;IONS<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ion_dynamics = 'bfgs' ,<BR>&nbsp;/<BR>ATOMIC_SPECIES<BR>&nbsp;&nbsp;&nbsp; C&nbsp;&nbsp; 12.01100&nbsp; C.pbe-rrkjus.UPF<BR>&nbsp;&nbsp;&nbsp; O&nbsp;&nbsp; 15.99940&nbsp; O.pbe-rrkjus.UPF<BR>ATOMIC_POSITIONS angstrom<BR>&nbsp;&nbsp;&nbsp; C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 4.000000000&nbsp;&nbsp;&nbsp; 4.000000000&nbsp;&nbsp;&nbsp; 4.000000000&nbsp;&nbsp;&nbsp; 1&nbsp; 1&nbsp; 1<BR>&nbsp;&nbsp;&nbsp; O&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 5.160000000&nbsp;&nbsp;&nbsp; 4.000000000&nbsp;&nbsp;&nbsp; 4.000000000&nbsp;&nbsp;&nbsp; 1&nbsp; 1&nbsp; 1<BR>&nbsp;&nbsp;&nbsp; O&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 2.840000000&nbsp;&nbsp;&nbsp; 4.000000000&nbsp;&nbsp;&nbsp; 4.000000000&nbsp;&nbsp;&nbsp; 1&nbsp; 1&nbsp; 1<BR>K_POINTS gamma</DIV>
<DIV>&nbsp;</DIV>
<DIV>(2) Performing scf calculations:</DIV>
<DIV>&nbsp;&amp;CONTROL<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; calculation = 'scf' ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; restart_mode = 'from_scratch' ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; outdir = '/home/olh/tmp/' ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; pseudo_dir = '/home/olh/pseudo/' ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; prefix = 'CO2' ,<BR>&nbsp;/<BR>&nbsp;&amp;SYSTEM<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ibrav = 1,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; celldm(1) = 15.11781502,<BR>&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 = 3,<BR>&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>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ecutwfc = 26 ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ecutrho = 260 ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; occupations = 'smearing' ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; degauss = 0.02 ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; smearing = 'methfessel-paxton' ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; nspin = 2 ,<BR>&nbsp;&nbsp; starting_magnetization(1) = 0.5,<BR>&nbsp;&nbsp; starting_magnetization(2) = 0.5,<BR>&nbsp;/<BR>&nbsp;&amp;ELECTRONS<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; conv_thr = 1.D-12 ,</DIV>
<DIV>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; mixing_beta = 0.4D0 ,<BR>&nbsp;/<BR>ATOMIC_SPECIES<BR>&nbsp;&nbsp;&nbsp; C&nbsp;&nbsp; 12.01100&nbsp; C.pbe-rrkjus.UPF<BR>&nbsp;&nbsp;&nbsp; O&nbsp;&nbsp; 15.99940&nbsp; O.pbe-rrkjus.UPF<BR>ATOMIC_POSITIONS angstrom<BR>&nbsp;&nbsp;&nbsp; C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 4.000000000&nbsp;&nbsp; 4.000000000&nbsp;&nbsp; 4.000000000&nbsp;&nbsp;&nbsp; 0&nbsp; 0&nbsp; 0<BR>&nbsp;&nbsp;&nbsp; O&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 5.172018724&nbsp;&nbsp; 4.000000000&nbsp;&nbsp; 4.000000000&nbsp;&nbsp;&nbsp; 0&nbsp; 0&nbsp; 0<BR>&nbsp;&nbsp;&nbsp; O&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 2.827981276&nbsp;&nbsp; 4.000000000&nbsp;&nbsp; 4.000000000&nbsp;&nbsp;&nbsp; 0&nbsp; 0&nbsp; 0<BR>K_POINTS automatic<BR>&nbsp; 16 16 16&nbsp; 1 1 1</DIV>
<DIV>&nbsp;</DIV>
<DIV>(3) Performing frequency calculation using PHONOS codes:</DIV>
<DIV>phonons of CO2 at Gamma<BR>&nbsp;&amp;inputph<BR>&nbsp; tr2_ph=1.0d-20,<BR>&nbsp; prefix='CO2',<BR>&nbsp; epsil=.false.,<BR>&nbsp; amass(1)=12.0110,<BR>&nbsp; amass(2)=15.9994,<BR>&nbsp; outdir='/home/olh/tmp/',<BR>&nbsp; fildyn='CO2.dynG',<BR>&nbsp;/<BR>0.0 0.0 0.0</DIV>
<DIV>&nbsp;</DIV>
<DIV>Finally, I obtained the frequency&nbsp;of CO2 molecule, output file as follows:</DIV>
<DIV>&nbsp;&nbsp;&nbsp; Convergence has been achieved</DIV>
<DIV>&nbsp;&nbsp;&nbsp; Number of q in the star =&nbsp;&nbsp;&nbsp; 1<BR>&nbsp;&nbsp;&nbsp; List of q in the star:<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp; 0.000000000&nbsp;&nbsp; 0.000000000&nbsp;&nbsp; 0.000000000</DIV>
<DIV>&nbsp;&nbsp;&nbsp; Diagonalizing the dynamical matrix</DIV>
<DIV>&nbsp;&nbsp;&nbsp; q = (&nbsp;&nbsp;&nbsp; 0.000000000&nbsp;&nbsp; 0.000000000&nbsp;&nbsp; 0.000000000 )</DIV>
<DIV>**************************************************************************<BR>&nbsp;&nbsp;&nbsp; omega( 1) =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -3.661345 [THz] =&nbsp;&nbsp;&nbsp; -122.130124 [cm-1]<BR>&nbsp;&nbsp;&nbsp; omega( 2) =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -3.661345 [THz] =&nbsp;&nbsp;&nbsp; -122.130124 [cm-1]<BR>&nbsp;&nbsp;&nbsp; omega( 3) =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 3.659957 [THz] =&nbsp;&nbsp;&nbsp;&nbsp; 122.083850 [cm-1]<BR>&nbsp;&nbsp;&nbsp; omega( 4) =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 7.879986 [THz] =&nbsp;&nbsp;&nbsp;&nbsp; 262.849774 [cm-1]<BR>&nbsp;&nbsp;&nbsp; omega( 5) =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 7.879986 [THz] =&nbsp;&nbsp;&nbsp;&nbsp; 262.849774 [cm-1]<BR>&nbsp;&nbsp;&nbsp; omega( 6) =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 19.746862 [THz] =&nbsp;&nbsp;&nbsp;&nbsp; 658.688787 [cm-1]<BR>&nbsp;&nbsp;&nbsp; omega( 7) =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 19.746862 [THz] =&nbsp;&nbsp;&nbsp;&nbsp; 658.688787 [cm-1]<BR>&nbsp;&nbsp;&nbsp; omega( 8) =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 38.945041 [THz] =&nbsp;&nbsp;&nbsp; 1299.075383 [cm-1]<BR>&nbsp;&nbsp;&nbsp; omega( 9) =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 70.000609 [THz] =&nbsp;&nbsp;&nbsp; 2334.984491 [cm-1]<BR>**************************************************************************<BR>&nbsp;&nbsp;&nbsp;&nbsp; Mode symmetry, D_4h(4/mmm) point group:</DIV>
<DIV>&nbsp;&nbsp;&nbsp;&nbsp; omega(&nbsp; 1 -&nbsp; 2) =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -122.1&nbsp; [cm-1]&nbsp;&nbsp; --&gt; E_g&nbsp; X_5&nbsp; M_5&nbsp;&nbsp; R<BR>&nbsp;&nbsp;&nbsp;&nbsp; omega(&nbsp; 3 -&nbsp; 3) =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 122.1&nbsp; [cm-1]&nbsp;&nbsp; --&gt; A_2u X_4' M_4'&nbsp; I<BR>&nbsp;&nbsp;&nbsp;&nbsp; omega(&nbsp; 4 -&nbsp; 5) =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 262.8&nbsp; [cm-1]&nbsp;&nbsp; --&gt; E_u&nbsp; X_5' M_5'&nbsp; I<BR>&nbsp;&nbsp;&nbsp;&nbsp; omega(&nbsp; 6 -&nbsp; 7) =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 658.7&nbsp; [cm-1]&nbsp;&nbsp; --&gt; E_u&nbsp; X_5' M_5'&nbsp; I<BR>&nbsp;&nbsp;&nbsp;&nbsp; omega(&nbsp; 8 -&nbsp; 8) =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1299.1&nbsp; [cm-1]&nbsp;&nbsp; --&gt; A_1g X_1&nbsp; M_1&nbsp;&nbsp; R<BR>&nbsp;&nbsp;&nbsp;&nbsp; omega(&nbsp; 9 -&nbsp; 9) =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 2335.0&nbsp; [cm-1]&nbsp;&nbsp; --&gt; A_2u X_4' M_4'&nbsp; I</DIV>
<DIV>&nbsp;**************************************************************************</DIV>
<DIV>&nbsp;&nbsp;&nbsp;&nbsp; PHONON&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; :&nbsp;&nbsp;&nbsp;&nbsp; 1h48m CPU time,&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 2h25m wall time<BR>......</DIV>
<DIV>There are still negative frequency in output file by&nbsp;increasing K_POINTS and reducing tr2_ph, anyone who help me will be appreciated. </DIV>
<DIV>In addition, I intend to calculate the zero point energy of system by PHNONS code, but it does not give the zero point energy and other thermodynamics properties in output file. How could I do? I need your help urgently.</DIV>
<DIV>&nbsp;</DIV>
<DIV>Thanks in advance </DIV>
<DIV>Best regards </DIV>
<DIV>Lihui </DIV><BR>--<BR></DIV>
<DIV>
<DIV>
<DIV><FONT face="Arial" size="3">======================================</FONT></DIV>
<DIV><FONT face="Arial" size="3">Lihui&nbsp;Ou</FONT></DIV>
<DIV><FONT face="Arial" size="3">PH.D Candidate&nbsp;in Electrochemistry&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </FONT></DIV>
<DIV><FONT face="Arial" size="3">College of Chemistry and Molecular Science&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </FONT></DIV>
<DIV><FONT face="Arial" size="3">Wuhan University,<EM>430072</EM>,Hubei Province,China&nbsp;</FONT></DIV>
<DIV><EM><FONT face="Arial" size="3">E-mail:</FONT><A href="mailto:oulihui666@126.com" target="_blank"><FONT face="Arial" size="3">oulihui666@126.com</FONT></A></EM></DIV>
<DIV><FONT face="Arial" size="3">======================================</FONT></DIV></DIV></DIV><br><br><span title="neteasefooter"/><hr/>
<a href="http://www.yeah.net/?from=footer">没有广告的终身免费邮箱,www.yeah.net</a>
</span>