
<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 3.2//EN">

<HTML>

<HEAD>

<META HTTP-EQUIV="Content-Type" CONTENT="text/html; charset=iso-8859-1">

<META NAME="Generator" CONTENT="MS Exchange Server version 6.5.7638.1">

<TITLE>scf calculation for isolated hydrogen atoms</TITLE>

</HEAD>

<BODY>

<!-- Converted from text/plain format -->


<P><FONT SIZE=2>Hi there,<BR>

I tried to calculate the ground state energy for isolated hydrogen atom with the following input,however I got the fermi energy level at -3.5661eV, while it should be -13.6eV.<BR>

<BR>

&amp;CONTROL.<BR>

&nbsp;&nbsp;&nbsp;&nbsp; calculation ='scf'<BR>

&nbsp;&nbsp;&nbsp; restart_mode = 'from_scratch'<BR>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; outdir = '/sfiwork/dsharma/eigenvaluesSi/wf_H/'<BR>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; pseudo_dir = '/sfiwork/dsharma/QE_Pseudos/'<BR>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; prefix = 'H'<BR>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; tstress = .true.<BR>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; tprnfor = .true.<BR>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; wf_collect =.true.<BR>

&nbsp;/<BR>

&nbsp;&amp;SYSTEM<BR>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ibrav = 1<BR>

&nbsp;&nbsp;&nbsp; celldm(1) = 15<BR>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; nat = 1<BR>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ntyp = 1<BR>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ecutwfc = 40<BR>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ecutrho = 160.0<BR>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; nbnd = 2,<BR>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; nspin=2,<BR>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; starting_magnetization(1)=0.5<BR>

&nbsp;&nbsp;&nbsp; nosym=.true.,<BR>

&nbsp;&nbsp;&nbsp; occupations='smearing',<BR>

&nbsp;&nbsp;&nbsp; smearing='mp',<BR>

&nbsp;&nbsp;&nbsp; degauss=0.01,<BR>

<BR>

&nbsp;/<BR>

&nbsp;&amp;ELECTRONS<BR>

&nbsp;&nbsp;&nbsp; diagonalization ='cg'<BR>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; mixing_mode = 'local-TF'<BR>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; conv_thr = 1.0d-6<BR>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; mixing_beta = 0.7<BR>

/<BR>

ATOMIC_SPECIES<BR>

&nbsp;&nbsp; H&nbsp;&nbsp; 1.00794&nbsp; H.pz-vbc.UPF<BR>

&nbsp; ATOMIC_POSITIONS angstroms<BR>

&nbsp;&nbsp;&nbsp; H&nbsp;&nbsp; 0.0000&nbsp; 0.0000 0.0000<BR>

K_POINTS automatic<BR>

&nbsp;1 1 1&nbsp;&nbsp; 0 0 0<BR>

<BR>

Thanks<BR>

<BR>

Dimpy<BR>

UCC<BR>

Ireland<BR>

<BR>

</FONT>

</P>


</BODY>

</HTML>