[Pw_forum] Plot fermi surface
oyxf328
oyxf328 at 126.com
Tue Apr 26 10:01:55 CEST 2016
Dear all
During the fermi surface plot of the PP/example02, in the last step of bands_FS.x , I found a
problem as follows, could you give some hints ?
There are no errors in the previous three steps.
First, I do a scf calculation using the pw.x <ni.scf_SP.in> ni.scf_SP.out, the input file ni.scf_SP.in is as follows:
&control
calculation='scf'
restart_mode='from_scratch',
prefix='ni',
pseudo_dir = '/data2/ouyang/software/espresso-5.1/espresso-5.1/pseudo/'
outdir= './tempdir/'
/
&system
ibrav=2, celldm(1) =6.48, nat=1, ntyp=1,
ecutwfc = 24.0, ecutrho = 288.0,
occupations='smearing',
smearing='methfessel-paxton',
degauss=0.02
nspin=2,
starting_magnetization(1)=0.8
/
&electrons
conv_thr = 1.0e-10
mixing_beta = 0.7
/
ATOMIC_SPECIES
Ni 58.69 Ni.pz-nd-rrkjus.UPF
ATOMIC_POSITIONS
Ni 0.0 0.0 0.0
K_POINTS {automatic}
8 8 8 0 0 0
The second step is generate the kpoints using the shell:
#
Sysname='ni'
nabc=' 16 16 16 '
n_start=3
n_last=6
#
E_Fermi=`grep Fermi ni.scf_SP.out | cut -c 26-36`
a1=`grep 'b(1)' ni.scf_SP.out | cut -c 24-54`
a2=`grep 'b(2)' ni.scf_SP.out | cut -c 24-54`
a3=`grep 'b(3)' ni.scf_SP.out | cut -c 24-54`
cat > kvecs_FS.in <<EOF
$a1
$a2
$a3
$nabc
$Sysname
EOF
and generate kpoints using this command: kvecs_FS.x <kvecs_FS.in> kvecs_FS.out
The third step is do a nscf calculation using this command pw.x <ni.fs_SP.in> ni.fs_SP.out
the input file ni.fs_SP.in is following:
&control
calculation='bands'
prefix='ni',
pseudo_dir = '/data2/ouyang/software/espresso-5.1/espresso-5.1/pseudo/'
outdir= './tempdir/'
verbosity='high'
/
&system
ibrav=2,
celldm(1) =6.48,
nat=1,
ntyp=1,
ecutwfc = 24.0,
ecutrho = 288.0,
nbnd=8
nspin=2,
starting_magnetization(1)=0.8
/
&electrons
diagonalization='cg'
conv_thr = 1.0e-10
mixing_beta = 0.7
/
ATOMIC_SPECIES
Ni 58.69 Ni.pz-nd-rrkjus.UPF
ATOMIC_POSITIONS
Ni 0.0 0.0 0.0
K_POINTS
4913
0.000000 0.000000 0.000000 1.00
-0.062500 0.062500 -0.062500 1.00
-0.125000 0.125000 -0.125000 1.00
-0.187500 0.187500 -0.187500 1.00
-0.250000 0.250000 -0.250000 1.00
here i omit other kpoints
Then mv ni.fs_SP.out Bands_SP.out and prepare the input_FS file as follows:
3 6
15.2788
ni
16 16 16
-1.000000 -1.000000 1.000000
1.000000 1.000000 1.000000
-1.000000 1.000000 -1.000000
the last step is using the bands_FS.x <Bans_SP.out> bands_fs.out to obtain the result ,but it always have this error:
Image PC Routine Line Source
bands_FS.x 0000000000491D2D Unknown Unknown Unknown
bands_FS.x 0000000000490835 Unknown Unknown Unknown
bands_FS.x 000000000044B819 Unknown Unknown Unknown
bands_FS.x 000000000041227A Unknown Unknown Unknown
bands_FS.x 0000000000411A70 Unknown Unknown Unknown
bands_FS.x 000000000040D5F2 Unknown Unknown Unknown
bands_FS.x 0000000000408261 MAIN__ 97 bands_FS.f90
bands_FS.x 0000000000407B4C Unknown Unknown Unknown
libc.so.6 0000003209A1D974 Unknown Unknown Unknown
bands_FS.x 0000000000407A59 Unknown Unknown Unknown
what's wrong ??? I did it following the example02 step by step.
I am looking forward to your kindly suggestions, thanks.
ouyang
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.quantum-espresso.org/pipermail/users/attachments/20160426/a13034de/attachment.html>
More information about the users
mailing list