<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv=Content-Type content="text/html; charset=us-ascii">
<meta name=Generator content="Microsoft Word 12 (filtered medium)">
<style>
<!--
/* Font Definitions */
@font-face
{font-family:PMingLiU;
panose-1:2 2 3 0 0 0 0 0 0 0;}
@font-face
{font-family:PMingLiU;
panose-1:2 2 3 0 0 0 0 0 0 0;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:PMingLiU;
panose-1:2 2 3 0 0 0 0 0 0 0;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:12.0pt;
font-family:"Calibri","sans-serif";}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
span.EmailStyle17
{mso-style-type:personal-compose;
font-family:"Calibri","sans-serif";
color:windowtext;}
span.apple-style-span
{mso-style-name:apple-style-span;}
.MsoChpDefault
{mso-style-type:export-only;
font-size:10.0pt;}
/* Page Definitions */
@page Section1
{size:612.0pt 792.0pt;
margin:72.0pt 90.0pt 72.0pt 90.0pt;}
div.Section1
{page:Section1;}
-->
</style>
<!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang=ZH-TW link=blue vlink=purple style='text-justify-trim:punctuation'>
<div class=Section1>
<p class=MsoNormal><span lang=EN-US>Dear all,<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US><o:p> </o:p></span></p>
<p class=MsoNormal><span lang=EN-US>I’ve got a problem when computing the
band structure for mono-layer CaC6.<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>The paper I followed is here (Arxiv
preprint arXiv:1105.3736, 2011 - arxiv.org).<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>They got their results from QE, too.<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US><o:p> </o:p></span></p>
<p class=MsoNormal><span lang=EN-US>The problem is not about QE running
process. It’s about the result I got from QE simulation.<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US><o:p> </o:p></span></p>
<p class=MsoNormal><span lang=EN-US>If I set the parameters (atomic positions
and celldm) as the same as paper.<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>The band structure I got was different from
paper.<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>Some bands are fine, but some are distorted.<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>For example, part of a band with flat shape
be tilted.<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>The link for my result is here. (<a
href="http://imgur.com/fUSna">http://imgur.com/fUSna</a>)<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US><o:p> </o:p></span></p>
<p class=MsoNormal><span lang=EN-US>It looks like something wrong with my
atomic structure.<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>But it isn’t. I had examined the
atomic structure and k-path for many times.<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US><o:p> </o:p></span></p>
<p class=MsoNormal><span lang=EN-US>To check the validity of atomic positions. I
also tried the same structure with VASP.<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>The result from VASP could fit the paper
pretty well.<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US><o:p> </o:p></span></p>
<p class=MsoNormal><span lang=EN-US>To check the validity of input files for
QE. (parameters and k-path)<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>I tried to plot the band structure of
mono-layer LiC6 with QE.<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>The band structure of LiC6 could fit the
paper, too.<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US><o:p> </o:p></span></p>
<p class=MsoNormal><span lang=EN-US>I had tried to relax the structure or use
other pseudo potential files, but still none better.<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US><o:p> </o:p></span></p>
<p class=MsoNormal><span lang=EN-US>Why I use the correct structure and correct
input files, but couldn’t get a correct result?<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>Have there been any wrong in my input setting?<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US><o:p> </o:p></span></p>
<p class=MsoNormal><span lang=EN-US>Could someone please give me any advices or
clues?<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>Any suggestions will be greatly
appreciated.<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US><o:p> </o:p></span></p>
<p class=MsoNormal><span lang=EN-US>Best regards,<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US><o:p> </o:p></span></p>
<p class=MsoNormal><span lang=EN-US>J. C. Huang IAMS, Taiwan<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US><o:p> </o:p></span></p>
<p class=MsoNormal><span lang=EN-US>================= the following is my input
files =================<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US><o:p> </o:p></span></p>
<p class=MsoNormal><span lang=EN-US>Scf.in<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>&control<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> calculation ='scf'<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> restart_mode
='from_scratch'<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> pseudo_dir
='/lustre/lwork/dyyang/bin/espresso-4.0.3/pseudo/'<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> outdir ='./'<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> prefix ='AuC'<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> disk_io='low',<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> wf_collect=.TRUE.,<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> /<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> &system<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> ibrav = 4,<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> celldm(1) = 8.050232711,<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> celldm(3) = 4.694835681,<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> nat = 7,<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> ntyp = 2,<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> ecutwfc = 45,<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> ecutrho = 400,<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> occupations =
"smearing",<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> smearing =
'gaussian'<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>
degauss = 0.002D0<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> /<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> &electrons<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> diagonalization='david'<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> mixing_mode = 'plain'<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> conv_thr = 1.0e-08<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> mixing_beta = 0.7<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> /<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> &ions<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> bfgs_ndim = 3,<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>ATOMIC_SPECIES<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> Ca 40.08 Ca.pz-n-vbc.UPF<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> C 12.0107 C.pz-vbc.UPF<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>ATOMIC_POSITIONS {crystal}<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> Ca
0.66667 0.33333 0.61200<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> C
0.33333 0.00000 0.50000<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> C
0.33333 0.33333 0.50000<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> C
0.00000 0.33333 0.50000<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> C
0.00000 0.66667 0.50000<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> C
0.66667 0.66667 0.50000<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> C
0.66667 0.00000 0.50000<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>K_POINTS (automatic)<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>16 16 1 0 0 0 <o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US><o:p> </o:p></span></p>
<p class=MsoNormal><span lang=EN-US>Band1.in<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US># self-consistent calculation<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> &control<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> calculation='bands',<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>
restart_mode='from_scratch',<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> prefix='AuC',<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> pseudo_dir =
'/lustre/lwork/dyyang/bin/espresso-4.0.3/pseudo/',<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> outdir='./'<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> tstress = .true.<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> tprnfor = .true.<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> forc_conv_thr = 1.0e-3<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> /<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> &system<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> ibrav = 4,<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> celldm(1) = 8.050232711,<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> celldm(3) = 4.694835681,<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> nat=7, ntyp=2<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> ecutwfc =45.0,<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> ecutrho =400,<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> occupations =
"smearing",<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>
smearing = "gaussian",<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>
degauss = 0.002D0,<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> nbnd= 20<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>/<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> &electrons<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> diagonalization='david'<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> mixing_mode = 'plain'<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> mixing_beta = 0.7<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> conv_thr = 1.0e-8<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>/<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> &ions<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> bfgs_ndim = 3,<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>ATOMIC_SPECIES<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> Ca 40.08 Ca.pz-n-vbc.UPF<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> C 12.0107 C.pz-vbc.UPF<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>ATOMIC_POSITIONS {crystal}<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> Ca
0.66667 0.33333 0.61200<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> C
0.33333 0.00000 0.50000<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> C
0.33333 0.33333 0.50000<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> C
0.00000 0.33333 0.50000<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> C
0.00000 0.66667 0.50000<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> C
0.66667 0.66667 0.50000<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> C
0.66667 0.00000 0.50000<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>K_POINTS {crystal}<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>49<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> 0.00000000 0.00000000
0.00000000 1<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> 0.03125000 0.00000000
0.00000000 2<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> 0.06250000 0.00000000
0.00000000 3<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> 0.09375000 0.00000000
0.00000000 4<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> 0.12500000 0.00000000
0.00000000 5<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> 0.15625000 0.00000000
0.00000000 6<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> 0.18750000 0.00000000
0.00000000 7<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> 0.21875000 0.00000000
0.00000000 8<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> 0.25000000 0.00000000
0.00000000 9<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> 0.28125000 0.00000000
0.00000000 10<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> 0.31250000 0.00000000
0.00000000 11<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> 0.34375000 0.00000000
0.00000000 12<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> 0.37500000 0.00000000
0.00000000 13<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> 0.40625000 0.00000000
0.00000000 14<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> 0.43750000 0.00000000
0.00000000 15<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> 0.46875000 0.00000000
0.00000000 16<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> 0.50000000 0.00000000
0.00000000 17<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> 0.48958333 0.02083333
0.00000000 18<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> 0.47916666 0.04166666
0.00000000 19<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> 0.46874999 0.06249999
0.00000000 20<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> 0.45833333 0.08333332
0.00000000 21<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> 0.44791666 0.10416666
0.00000000 22<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> 0.43749999 0.12499999
0.00000000 23<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> 0.42708332 0.14583332
0.00000000 24<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> 0.41666665 0.16666665
0.00000000 25<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> 0.40624998 0.18749998
0.00000000 26<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> 0.39583331 0.20833331
0.00000000 27<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> 0.38541664 0.22916664
0.00000000 28<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> 0.37499997 0.24999997
0.00000000 29<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> 0.36458331 0.27083331
0.00000000 30<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> 0.35416664 0.29166664
0.00000000 31<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> 0.34374997 0.31249997
0.00000000 32<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> 0.33333330 0.33333330
0.00000000 33<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> 0.31249997 0.31249997
0.00000000 34<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> 0.29166664 0.29166664
0.00000000 35<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> 0.27083331 0.27083331
0.00000000 36<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> 0.24999997 0.24999997
0.00000000 37<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> 0.22916664 0.22916664 0.00000000
38<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> 0.20833331 0.20833331
0.00000000 39<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> 0.18749998 0.18749998
0.00000000 40<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> 0.16666665 0.16666665
0.00000000 41<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> 0.14583332 0.14583332
0.00000000 42<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> 0.12499999 0.12499999
0.00000000 43<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> 0.10416666 0.10416666
0.00000000 44<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> 0.08333333 0.08333333
0.00000000 45<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> 0.06249999 0.06249999
0.00000000 46<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> 0.04166666 0.04166666
0.00000000 47<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> 0.02083333 0.02083333
0.00000000 48<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> 0.00000000 0.00000000
0.00000000 49<o:p></o:p></span></p>
</div>
</body>
</html>