[Pw_forum] hydrogen adsorption on graphene

Eyvaz Isaev eyvaz_isaev at yahoo.com
Thu Jul 31 12:33:49 CEST 2008


Hi,

Please read carefully Troubleshootings in the UserGuide (page 55 or around). This is very very very ... famous question.

Bests,
Eyvaz.

 -------------------------------------------------------------------
Prof. Eyvaz Isaev, 
Theoretical Physics Department, Moscow State Institute of Steel & Alloys, Russia, 
Department of Physics, Chemistry, and Biology (IFM), Linkoping University, Sweden 
Condensed Matter Theory Group, Uppsala University, Sweden 
Eyvaz.Isaev at fysik.uu.se, isaev at ifm.liu.se, eyvaz_isaev at yahoo.com



----- Original Message ----
From: "s305igun at mail.chem.itb.ac.id" <s305igun at mail.chem.itb.ac.id>
To: wangxinquan at tju.edu.cn; PWSCF Forum <pw_forum at pwscf.org>
Sent: Thursday, July 31, 2008 11:55:55 AM
Subject: [Pw_forum] hydrogen adsorption on graphene

Dear PW_Forum

I need your help, I was calculated hydrogen adsorption on graphene layer
and this is my trouble output:

starting charge  281.99714, renormalised to  282.00000
     Starting wfc are  282 atomic wfcs

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
     from cdiaghg : error #       296
     info =/= 0
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

     stopping ...

Would you like to explain to me about info =/= 0?

thanks very mach. and this is my input list:

#!/bin/bash

. ../environment_variables

PW_COMMAND=$PW_ROOT/bin/pw.x
PP_COMMAND=$PW_ROOT/bin/pp.x

##### Generate input files for scf calculation

cat <<EOF > h2-graphite.scf.in
&CONTROL
  calculation   = 'scf'
  title         = 'h2 / carbon layers slab - for demonstrative purposes only'
  restart_mode  = 'from_scratch'
  prefix        = 'h2-graphite'
  outdir        = '$TMP_DIR/'
  pseudo_dir    = '$PSEUDO_DIR/'
/
&SYSTEM
    ibrav     = 0,
    celldm(1) = 5.3d0
    nat       = 72
    ntyp      = 2
    ecutwfc   = 16.0d0
    ecutrho   = 160.0d0
    occupations = 'smearing'
    degauss     = 0.03D0
    smearing    = 'mp'
/
&ELECTRONS
    conv_thr    = 1.0d-5
    mixing_mode = 'plain'
    mixing_beta = 0.3d0
/
ATOMIC_SPECIES
C     12.0107  C.pbe-rrkjus.UPF
H      1.0079  H.pbe-rrkjus.UPF
ATOMIC_POSITIONS angstrom
H       8.741687    2.884450    1.000000
H       8.741715    2.884450    1.500000
C       1.249543   -0.711074    0.000000
C       0.000947   -0.016393    0.000000
C       0.000844    1.443145    0.000000
C       1.249381    2.159954    0.000000
C       1.249409    3.608946    0.000000
C       0.000843    4.325668    0.000000
C       0.000939    5.785283    0.000000
C       1.249515    6.479991    0.000000
C       1.249981   -1.825631    0.000000
C       1.249981    7.594413    0.000000
C       3.746978   -0.711074    0.000000
C       2.498382   -0.016393    0.000000
C       2.498279    1.443145    0.000000
C       3.746816    2.159954    0.000000
C       3.746844    3.608946    0.000000
C       2.498278    4.325668    0.000000
C       2.498374    5.785283    0.000000
C       3.746950    6.479991    0.000000
C       3.747416   -1.825631    0.000000
C       3.747416    7.594413    0.000000
C       6.244413   -0.711074    0.000000
C       4.995817   -0.016393    0.000000
C       4.995714    1.443145    0.000000
C       6.244251    2.159954    0.000000
C       6.244279    3.608946    0.000000
C       4.995713    4.325668    0.000000
C       4.995809    5.785283    0.000000
C       6.244385    6.479991    0.000000
C       6.244851   -1.825631    0.000000
C       6.244851    7.594413    0.000000
C       8.741849   -0.711074    0.000000
C       7.493253   -0.016393    0.000000
C       7.493150    1.443145    0.000000
C       8.741687    2.159954    0.000000
C       8.741715    3.608946    0.000000
C       7.493149    4.325668    0.000000
C       7.493245    5.785283    0.000000
C       8.741821    6.479991    0.000000
C       8.742287   -1.825631    0.000000
C       8.742287    7.594413    0.000000
C      11.239284   -0.711074    0.000000
C       9.990688   -0.016393    0.000000
C       9.990585    1.443145    0.000000
C      11.239122    2.159954    0.000000
C      11.239150    3.608946    0.000000
C       9.990584    4.325668    0.000000
C       9.990680    5.785283    0.000000
C      11.239256    6.479991    0.000000
C      11.239722   -1.825631    0.000000
C      11.239722    7.594413    0.000000
C      13.736719   -0.711074    0.000000
C      12.488123   -0.016393    0.000000
C      12.488020    1.443145    0.000000
C      13.736557    2.159954    0.000000
C      13.736585    3.608946    0.000000
C      12.488019    4.325668    0.000000
C      12.488115    5.785283    0.000000
C      13.736691    6.479991    0.000000
C      13.737157   -1.825631    0.000000
C      13.737157    7.594413    0.000000
C      16.234154   -0.711074    0.000000
C      14.985558   -0.016393    0.000000
C      14.985455    1.443145    0.000000
C      16.233992    2.159954    0.000000
C      16.234020    3.608946    0.000000
C      14.985454    4.325668    0.000000
C      14.985550    5.785283    0.000000
C      16.234126    6.479991    0.000000
C      16.234592   -1.825631    0.000000
C      16.234592    7.594413    0.000000
K_POINTS { automatic }
     2   1   1  1  1  1
CELL_PARAMETERS
   1.000000    0.000000    0.000000
   0.000000    1.144031    0.000000
   0.000000    0.000000    0.858023
EOF

cat <<EOF > h2.scf.in
&CONTROL
  calculation   = 'scf'
  title         = 'h2 - for demonstrative purposes only'
  restart_mode  = 'from_scratch'
  prefix        = 'h2'
  outdir        = '$TMP_DIR'
  pseudo_dir    = '$PSEUDO_DIR'
/
&SYSTEM
    ibrav     = 0,
    celldm(1) = 5.3d0
    nat       = 2
    ntyp      = 1
    ecutwfc   = 16.0d0
    ecutrho   = 160.0d0
    occupations = 'smearing'
    degauss     = 0.03D0
    smearing    = 'mp'
/
&ELECTRONS
    conv_thr    = 1.0d-5
    mixing_mode = 'plain'
    mixing_beta = 0.3d0
/
ATOMIC_SPECIES
H      1.0079  H.pbe-rrkjus.UPF
ATOMIC_POSITIONS angstrom
H       8.741687    2.884450    1.000000
H       8.741715    2.884450    1.500000
K_POINTS { automatic }
     2   2   1  1  1  1
CELL_PARAMETERS
   1.000000    0.000000    0.000000
   0.000000    1.144031    0.000000
   0.000000    0.000000    0.858023
EOF

cat <<EOF > graphite.scf.in
&CONTROL
  calculation   = 'scf'
  title         = 'graphite layers slab - for demonstrative purposes only'
  restart_mode  = 'from_scratch'
  prefix        = 'graphite'
  outdir        = '$TMP_DIR'
  pseudo_dir    = '$PSEUDO_DIR'
/
&SYSTEM
    ibrav     = 0,
    celldm(1) = 5.3d0
    nat       = 70
    ntyp      = 1
    ecutwfc   = 16.0d0
    ecutrho   = 160.0d0
    occupations = 'smearing'
    degauss     = 0.03D0
    smearing    = 'mp'
/
&ELECTRONS
    conv_thr    = 1.0d-5
    mixing_mode = 'plain'
    mixing_beta = 0.3d0
/
ATOMIC_SPECIES
C     12.0107  C.pbe-rrkjus.UPF
ATOMIC_POSITIONS angstrom
C       1.249543   -0.711074    0.000000
C       0.000947   -0.016393    0.000000
C       0.000844    1.443145    0.000000
C       1.249381    2.159954    0.000000
C       1.249409    3.608946    0.000000
C       0.000843    4.325668    0.000000
C       0.000939    5.785283    0.000000
C       1.249515    6.479991    0.000000
C       1.249981   -1.825631    0.000000
C       1.249981    7.594413    0.000000
C       3.746978   -0.711074    0.000000
C       2.498382   -0.016393    0.000000
C       2.498279    1.443145    0.000000
C       3.746816    2.159954    0.000000
C       3.746844    3.608946    0.000000
C       2.498278    4.325668    0.000000
C       2.498374    5.785283    0.000000
C       3.746950    6.479991    0.000000
C       3.747416   -1.825631    0.000000
C       3.747416    7.594413    0.000000
C       6.244413   -0.711074    0.000000
C       4.995817   -0.016393    0.000000
C       4.995714    1.443145    0.000000
C       6.244251    2.159954    0.000000
C       6.244279    3.608946    0.000000
C       4.995713    4.325668    0.000000
C       4.995809    5.785283    0.000000
C       6.244385    6.479991    0.000000
C       6.244851   -1.825631    0.000000
C       6.244851    7.594413    0.000000
C       8.741849   -0.711074    0.000000
C       7.493253   -0.016393    0.000000
C       7.493150    1.443145    0.000000
C       8.741687    2.159954    0.000000
C       8.741715    3.608946    0.000000
C       7.493149    4.325668    0.000000
C       7.493245    5.785283    0.000000
C       8.741821    6.479991    0.000000
C       8.742287   -1.825631    0.000000
C       8.742287    7.594413    0.000000
C      11.239284   -0.711074    0.000000
C       9.990688   -0.016393    0.000000
C       9.990585    1.443145    0.000000
C      11.239122    2.159954    0.000000
C      11.239150    3.608946    0.000000
C       9.990584    4.325668    0.000000
C       9.990680    5.785283    0.000000
C      11.239256    6.479991    0.000000
C      11.239722   -1.825631    0.000000
C      11.239722    7.594413    0.000000
C      13.736719   -0.711074    0.000000
C      12.488123   -0.016393    0.000000
C      12.488020    1.443145    0.000000
C      13.736557    2.159954    0.000000
C      13.736585    3.608946    0.000000
C      12.488019    4.325668    0.000000
C      12.488115    5.785283    0.000000
C      13.736691    6.479991    0.000000
C      13.737157   -1.825631    0.000000
C      13.737157    7.594413    0.000000
C      16.234154   -0.711074    0.000000
C      14.985558   -0.016393    0.000000
C      14.985455    1.443145    0.000000
C      16.233992    2.159954    0.000000
C      16.234020    3.608946    0.000000
C      14.985454    4.325668    0.000000
C      14.985550    5.785283    0.000000
C      16.234126    6.479991    0.000000
C      16.234592   -1.825631    0.000000
C      16.234592    7.594413    0.000000
K_POINTS { automatic }
     2   1   1  1  1  1
CELL_PARAMETERS
   1.000000    0.000000    0.000000
   0.000000    1.144031    0.000000
   0.000000    0.000000    0.858023
EOF

##### Input files generated. Now execute!

for PREFIX in graphite h2 h2-graphite; do
echo Now doing $PREFIX ...
$PW_COMMAND < $PREFIX.scf.in > $PREFIX.scf.out
cat <<EOF > $PREFIX.ppin
&inputPP
   prefix   = '$PREFIX'
   outdir   = '$TMP_DIR'
   filplot  = '$PREFIX.charge'
   plot_num = 0
/
EOF
$PP_COMMAND < $PREFIX.ppin > $PREFIX.ppout
done
echo Done

oke everybody, thanks for your help

best regards

Rahmat Gunawan
Students of Chemistry Dept.
Lab of Theoretical and Computational Chemistry
Faculty of Mathematics and Natural Sciences
Bandung Institute of Technology
_______________________________________________
Pw_forum mailing list
Pw_forum at pwscf.org
http://www.democritos.it/mailman/listinfo/pw_forum



      


More information about the Pw_forum mailing list