<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Verdana
}
--></style>
</head>
<body class='hmmessage'>
Dear Linh,<BR>
&nbsp;<BR>
Thanks a&nbsp;lot!<BR>
&nbsp;<BR>
&nbsp;Could you tell me how to get the <SPAN lang=EN>NORM-CONSERVING pseudopotential of Cd (Cadmium) which can be used in hybrid-functionals calculation?&nbsp;</SPAN><BR>
<SPAN lang=EN>&nbsp;I only find the ultrosoft&nbsp;</SPAN>&nbsp;pseudopotential of Cd which&nbsp;has not been&nbsp;implemented in the hybrid-functionals calculation.<BR>&nbsp;<BR>
Thank you again!&nbsp;<BR>
&nbsp;<BR>
&nbsp;Rgds<BR>Jianchun Wu<BR>Departmant of Physics<BR>Soochow Univ. China<BR>215006<BR>
&nbsp;<BR>&gt; Date: Tue, 13 Apr 2010 10:25:09 +0200<BR>&gt; From: nnlinh@sissa.it<BR>&gt; To: pw_forum@pwscf.org<BR>&gt; Subject: Re: [Pw_forum] How to calculate DOS by Hybrid functionals.<BR>&gt; <BR>&gt; <BR>&gt; I rerun your calculations basing on your input files with a modify for <BR>&gt; pseudo directory, and they work well.<BR>&gt; <BR>&gt; So, pls check carefully the input file before starting calculations.<BR>&gt; <BR>&gt; Good luck,<BR>&gt; <BR>&gt; Linh<BR>&gt; <BR>&gt; wujianchun wrote:<BR>&gt; &gt; Dear All,<BR>&gt; &gt;<BR>&gt; &gt; I want to get the DOS of Si by use of /Hybrid functionals. /The first <BR>&gt; &gt; step scf calculation can be run normally and the result is same as the <BR>&gt; &gt; EXX_example. but the nscf calculation reports error.<BR>&gt; &gt;<BR>&gt; &gt; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%<BR>&gt; &gt; from read_namelists : error # 1<BR>&gt; &gt; reading namelist control<BR>&gt; &gt; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%<BR>&gt; &gt;<BR>&gt; &gt; Below is my input file. Could someone help me to deal with DOS <BR>&gt; &gt; calculation?<BR>&gt; &gt;<BR>&gt; &gt; Thanks in advance!<BR>&gt; &gt;<BR>&gt; &gt; Rgds<BR>&gt; &gt; Jianchun Wu<BR>&gt; &gt;<BR>&gt; &gt; Departmant of Physics<BR>&gt; &gt; Soochow Univ. China<BR>&gt; &gt; 215006<BR>&gt; &gt;<BR>&gt; &gt; _____________________________________________________________________________<BR>&gt; &gt; cat &gt; CdS-0-o.pw.inp &lt;&lt; EOF<BR>&gt; &gt; &amp;CONTROL<BR>&gt; &gt; title = 'Si-0',<BR>&gt; &gt; calculation = 'scf',<BR>&gt; &gt; prefix = 'Si-0.pw',<BR>&gt; &gt; outdir = './Si-0', ! Need to create this directory<BR>&gt; &gt; wf_collect=.true., !<BR>&gt; &gt; restart_mode='from_scratch', ! Change to restart if you want to restart<BR>&gt; &gt; max_seconds=255000 ! Just under 259200 == 72 hours<BR>&gt; &gt; /<BR>&gt; &gt; &amp;system<BR>&gt; &gt; ibrav= 2, celldm(1) =10.20, nat= 2, ntyp= 1,<BR>&gt; &gt; ecutwfc =12.0, nbnd = 8,<BR>&gt; &gt; input_dft='pbe0', nqx1 = 1, nqx2 = 1, nqx3 = 1,<BR>&gt; &gt; /<BR>&gt; &gt; &amp;electrons<BR>&gt; &gt; mixing_beta = 0.7<BR>&gt; &gt; /<BR>&gt; &gt; ATOMIC_SPECIES<BR>&gt; &gt; Si 28.086 Si.pz-vbc.UPF<BR>&gt; &gt; ATOMIC_POSITIONS<BR>&gt; &gt; Si 0.00 0.00 0.00<BR>&gt; &gt; Si 0.25 0.25 0.25<BR>&gt; &gt; K_POINTS<BR>&gt; &gt; 10<BR>&gt; &gt; 0.1250000 0.1250000 0.1250000 1.00<BR>&gt; &gt; 0.1250000 0.1250000 0.3750000 3.00<BR>&gt; &gt; 0.1250000 0.1250000 0.6250000 3.00<BR>&gt; &gt; 0.1250000 0.1250000 0.8750000 3.00<BR>&gt; &gt; 0.1250000 0.3750000 0.3750000 3.00<BR>&gt; &gt; 0.1250000 0.3750000 0.6250000 6.00<BR>&gt; &gt; 0.1250000 0.3750000 0.8750000 6.00<BR>&gt; &gt; 0.1250000 0.6250000 0.6250000 3.00<BR>&gt; &gt; 0.3750000 0.3750000 0.3750000 1.00<BR>&gt; &gt; 0.3750000 0.3750000 0.6250000 3.00<BR>&gt; &gt; EOF<BR>&gt; &gt; cd $QSUB_WORKDIR<BR>&gt; &gt; tot_cores=$(($QSUB_VNODES*$QSUB_CPUS))<BR>&gt; &gt; time mpirun -perhost $QSUB_CPUS -n $tot_cores <BR>&gt; &gt; /apps/Espresso/espresso-4.1.2-i11scaEXX/bin/pw.x -in CdS-0-o.pw.inp &gt;&gt; <BR>&gt; &gt; CdS-0-o.pw.out<BR>&gt; &gt; #Create this file when done<BR>&gt; &gt; touch CdS-0-o.pw.done<BR>&gt; &gt;<BR>&gt; &gt; #DOS calculation for VCd<BR>&gt; &gt; cat &gt; ni.dos.in &lt;&lt; EOF<BR>&gt; &gt; &amp;CONTROL<BR>&gt; &gt; title = 'Si-0',<BR>&gt; &gt; calculation = 'nscf',<BR>&gt; &gt; prefix = 'Si-0.pw',<BR>&gt; &gt; outdir = './Si-0', ! Need to create this directory<BR>&gt; &gt; wf_collect=.true., ! Collect data all in one place<BR>&gt; &gt; restart_mode='from_scratch', ! Change to restart if you want to restart<BR>&gt; &gt; max_seconds=255000 ! Just under 259200 == 72 hours<BR>&gt; &gt; /<BR>&gt; &gt; /<BR>&gt; &gt; &amp;system<BR>&gt; &gt; ibrav= 2, celldm(1) =10.20, nat= 2, ntyp= 1,<BR>&gt; &gt; ecutwfc =12.0, nbnd = 8,<BR>&gt; &gt; input_dft='pbe0', nqx1 = 1, nqx2 = 1, nqx3 = 1,<BR>&gt; &gt; /<BR>&gt; &gt; &amp;electrons<BR>&gt; &gt; mixing_beta = 0.7<BR>&gt; &gt; /<BR>&gt; &gt; ATOMIC_SPECIES<BR>&gt; &gt; Si 28.086 Si.pz-vbc.UPF<BR>&gt; &gt; ATOMIC_POSITIONS<BR>&gt; &gt; Si 0.00 0.00 0.00<BR>&gt; &gt; Si 0.25 0.25 0.25<BR>&gt; &gt; K_POINTS<BR>&gt; &gt; 10<BR>&gt; &gt; 0.1250000 0.1250000 0.1250000 1.00<BR>&gt; &gt; 0.1250000 0.1250000 0.3750000 3.00<BR>&gt; &gt; 0.1250000 0.1250000 0.6250000 3.00<BR>&gt; &gt; 0.1250000 0.1250000 0.8750000 3.00<BR>&gt; &gt; 0.1250000 0.3750000 0.3750000 3.00<BR>&gt; &gt; 0.1250000 0.3750000 0.6250000 6.00<BR>&gt; &gt; 0.1250000 0.3750000 0.8750000 6.00<BR>&gt; &gt; 0.1250000 0.6250000 0.6250000 3.00<BR>&gt; &gt; 0.3750000 0.3750000 0.3750000 1.00<BR>&gt; &gt; 0.3750000 0.3750000 0.6250000 3.00<BR>&gt; &gt; EOF<BR>&gt; &gt;<BR>&gt; &gt; cat &gt; ni.dos2.in &lt;&lt; EOF<BR>&gt; &gt; &amp;inputpp<BR>&gt; &gt; prefix = 'Si-0.pw',<BR>&gt; &gt; outdir = './Si-0',<BR>&gt; &gt; fildos='Si.dos',<BR>&gt; &gt; Emin=-8.0, Emax=12.0, DeltaE=0.01<BR>&gt; &gt; /<BR>&gt; &gt; EOF<BR>&gt; &gt;<BR>&gt; &gt; cd $QSUB_WORKDIR<BR>&gt; &gt; tot_cores=$(($QSUB_VNODES*$QSUB_CPUS))<BR>&gt; &gt; time mpirun -perhost $QSUB_CPUS -n $tot_cores <BR>&gt; &gt; /apps/Espresso/espresso-4.1.2-i11scaEXX/bin/pw.x &lt; ni.dos.in &gt; ni.dos.out<BR>&gt; &gt; time mpirun -perhost $QSUB_CPUS -n $tot_cores <BR>&gt; &gt; /apps/Espresso/espresso-4.1.2-i11scaEXX/bin/dos.x &lt; ni.dos2.in &gt; <BR>&gt; &gt; ni.dos2.out<BR>&gt; &gt;<BR>&gt; &gt; touch dos.done<BR>&gt; &gt;<BR>&gt; &gt;<BR>&gt; &gt; ------------------------------------------------------------------------<BR>&gt; &gt; ??+??+?? ??????,??MSN????! ????! &lt;http://3g.msn.cn/&gt;<BR>&gt; &gt; ------------------------------------------------------------------------<BR>&gt; &gt;<BR>&gt; &gt; _______________________________________________<BR>&gt; &gt; Pw_forum mailing list<BR>&gt; &gt; Pw_forum@pwscf.org<BR>&gt; &gt; http://www.democritos.it/mailman/listinfo/pw_forum<BR>&gt; &gt; <BR>&gt; <BR>&gt; <BR>&gt; -- <BR>&gt; -----------------------------------------------------<BR>&gt; Nguyen Ngoc Linh, PhD Student<BR>&gt; c/o: SISSA &amp; CNR-INFM Democritos,<BR>&gt; via Beirut 2-4, 34014 Trieste (Italy)<BR>&gt; email: nnlinh@sissa.it<BR>&gt; phone: +39 04 03787 319<BR>&gt; skype: ngoclinh84phys<BR>&gt; -----------------------------------------------------<BR>&gt; "The physics is theoretical but the fun is real"<BR>&gt; <BR>&gt; <BR>&gt; _______________________________________________<BR>&gt; Pw_forum mailing list<BR>&gt; Pw_forum@pwscf.org<BR>&gt; http://www.democritos.it/mailman/listinfo/pw_forum<BR>                                               <br /><hr />搜索本应是彩色的,快来体验新一代搜索引擎-必应,精美图片每天换哦! <a href='http://cn.bing.com/?form=CRMADS ' target='_new'>立即试用!</a></body>
</html>