|
|
我不是高版本的，我是4.3.2所以應(yīng)該nscf，example也是這樣的，不過我遇到了新的我問題，麻煩幫我看看：
您好，我在測試ecut和K點的時候遇到了新的問題，我是這樣做的，分別設(shè)ecut=30 40 50 60 70 80 90 100 110 120,然后在每個ecut值下面慢慢增加K點，但是奇怪的是，在ecut=50 60和70的時候自洽運算會出現(xiàn)對稱性操作丟失的警告：
warning: symmetry operation # 3 not allowed. fractional translation:
0.5000000 -0.5000000 0.5000000 in crystal coordinates
warning: symmetry operation # 4 not allowed. fractional translation:
-0.5000000 0.5000000 0.5000000 in crystal coordinates
warning: symmetry operation # 7 not allowed. fractional translation:
-0.5000000 0.5000000 -0.5000000 in crystal coordinates
warning: symmetry operation # 8 not allowed. fractional translation:
0.5000000 -0.5000000 -0.5000000 in crystal coordinates
warning: symmetry operation # 11 not allowed. fractional translation:
-0.5000000 0.5000000 -0.5000000 in crystal coordinates
warning: symmetry operation # 12 not allowed. fractional translation:
0.5000000 -0.5000000 -0.5000000 in crystal coordinates
warning: symmetry operation # 15 not allowed. fractional translation:
0.5000000 -0.5000000 0.5000000 in crystal coordinates
warning: symmetry operation # 16 not allowed. fractional translation:
-0.5000000 0.5000000 0.5000000 in crystal coordinates
我看到小木蟲上有說通過增加ecut可以解決，可是為什么30 和40又能跑呢？這個問題怎么解決，我同時貼出我的IN文件：
&CONTROL
calculation = 'scf' ,
restart_mode = 'from_scratch' ,
outdir = '/home/liushichen/work/pwscf/tmp/' ,
pseudo_dir = '/home/liushichen/tools/codes/espresso-4.3.2/pseudo/' ,
prefix = 'CrO2' ,
/
&SYSTEM
ibrav = 0,
nat = 6,
ntyp = 2,
ecutwfc = 50 ,
ecutrho = 200 ,
occupations = 'smearing' ,
degauss = 0.02 ,
smearing = 'gaussian' ,
nspin = 2 ,
starting_magnetization(1) = 0.5,
lda_plus_u = .true. ,
Hubbard_U(1) = 3,
/
&ELECTRONS
conv_thr = 1.0e-10 ,
mixing_beta = 0.7 ,
/
CELL_PARAMETERS cubic
8.350300000 0.000000000 0.000000000
0.000000000 8.350300000 0.000000000
0.000000000 0.000000000 5.502600000
ATOMIC_SPECIES
Cr 52.00000 Cr.pz-hgh.UPF
O 16.00000 O.pz-hgh.UPF
ATOMIC_POSITIONS crystal
Cr 0.000000000 0.000000000 0.000000000
Cr 0.500000000 0.500000000 0.500000000
O 0.303000000 0.303000000 0.000000000
O 0.697000000 0.697000000 0.000000000
O 0.803000000 0.197000000 0.500000000
O 0.197000000 0.803000000 0.500000000
K_POINTS automatic
6 6 9 0 0 0
2，第二個問題是：我發(fā)現(xiàn)計算結(jié)果的out文件中的原子位置的Z方向有移動，但是我沒有relax阿，我只是跑的自洽阿，這是為什么，貼出那部分文件：
Cartesian axes
site n. atom positions (alat units)
1 Cr tau( 1) = ( 0.0000000 0.0000000 0.0000000 )
2 Cr tau( 2) = ( 0.5000000 0.5000000 0.3294852 )
3 O tau( 3) = ( 0.3030000 0.3030000 0.0000000 )
4 O tau( 4) = ( 0.6970000 0.6970000 0.0000000 )
5 O tau( 5) = ( 0.8030000 0.1970000 0.3294852 )
6 O tau( 6) = ( 0.1970000 0.8030000 0.3294852 )
麻煩幫我看下，謝謝！ |
|
回復此樓
10+3