小木虫 --- 500万硕博科研人员喜爱的学术科研平台
&&查看话题
Problem - nscf - c_bands: eigenvalues not converged
采用pwscf计算费米面时，前面若干步顺利通过，但在band_FS.x步骤时报错
./bands_FS.x bands_fs.out
forrtl: severe (59): list-directed I/O syntax error, unit 5, file stdin
Image& && && && &&&PC& && &&&Routine& && && && &Line& && &&&Source
bands_FS.x& && && &0809DB77&&Unknown& && && && && &Unknown&&Unknown
bands_FS.x& && && &&&Unknown& && && && && &Unknown&&Unknown
bands_FS.x& && && &08071A96&&Unknown& && && && && &Unknown&&Unknown
bands_FS.x& && && &08050B4C&&Unknown& && && && && &Unknown&&Unknown
bands_FS.x& && && &&&Unknown& && && && && &Unknown&&Unknown
bands_FS.x& && && &08061C93&&Unknown& && && && && &Unknown&&Unknown
请高手赐教。【真不希望该问题无解】
类似的问题还出现在：
http://qe-forge.org/pipermail/pw_forum/2011-October/096758.html；
http://www.democritos.it/pipermail/pw_forum/2010-July/017605.html；
以及小木虫上
http://emuch.net/html/7669.html
附件2：band_FS.x程序需要的input_FS如下：
& & 9.0728
&&1...000000
&&0...000000
&&0...628094
附3：非自洽计算所用的输入和输出结果文件如下：
3. 设定nbnd=16
& &calculation='nscf',
& &prefix='Sc'
& &pseudo_dir = '/home/users/spclibtt/home/yexq/shanghai-opt/pw-work/pseudo/',
& &outdir='/home/users/spclibtt/home/yexq/shanghai-opt/pw-work/environment/tmp/'
& &ibrav = 0,
& &nat=&&2,
& &ntyp= 1,
& &ecutwfc = 90,
& &ecutrho = 400.00,
& &nbnd=16
&electrons
& &diagonalization='cg'
& &mixing_beta = 0.7
& &conv_thr =&&1.0d-10
#cell_dynamics= 'damp-pr'
ATOMIC_SPECIES
Sc& & 44.9559& & Sc.pbe-nsp-van.UPF
CELL_PARAMETERS {bohr}
ATOMIC_POSITIONS {crystal}
Sc 0.4 0.0000
Sc 0.6 0.0000
& & 0.000000& & 0.000000& & 0..00
& & 0.000000& & 0.000000& & 0..00
& & 0.000000& & 0.000000& & 0..00
& & 0.000000& & 0.000000& & 0..00
& & 0.000000& & 0.000000& & 0..00
& & 0.000000& & 0.000000& & 0..00
& &&&Program PWSCF v.5.0.2 (svn rev. 9392) starts on 13Mar2013 at 11:17:57
& &&&This program is part of the open-source Quantum ESPRESSO suite
& &&&for quantum sim please cite
& && && &"P. Giannozzi et al., J. Phys.:Condens. Matter 21 09);
& && && & URL http://www.quantum-espresso.org",
& &&&in publications or presentations arising from this work. More details at
& &&&http://www.quantum-espresso.org/quote.php
& &&&Parallel version (MPI), running on& & 16 processors
& &&&R & G space division:&&proc/nbgrp/npool/nimage =& && &16
& &&&Current dimensions of program PWSCF are:
& &&&Max number of different atomic species (ntypx) = 10
& &&&Max number of k-points (npk) =&&40000
& &&&Max angular momentum in pseudopotentials (lmaxx) =&&3
& &&&Waiting for input...
& &&&Reading input from standard input
& &&&Atomic positions and unit cell read from directory:
& &&&/home/users/spclibtt/home/yexq/shanghai-opt/pw-work/environment/tmp/Sc.save/
& &&&Subspace diagonalization in iterative solution of the eigenvalue problem:
& &&&a serial algorithm will be used
& &&&Parallelization info
& &&&--------------------
& &&&sticks:& &dense&&smooth& &&&PW& &&&G-vecs:& & dense& &smooth& && &PW
& &&&Min& && && & 66& && &59& &&&22& && && && && &&&2846& &&&2423& &&&554
& &&&Max& && && & 67& && &60& &&&23& && && && && &&&2855& &&&2438& &&&561
& &&&Sum& && &&&1069& &&&955& & 361& && && && && & 45639& & 38911& & 8929
& &&&bravais-lattice index& &&&=& && && && &0
& &&&lattice parameter (alat)&&=& && & 6.2548&&a.u.
& &&&unit-cell volume& && && & =& &&&337.4042 (a.u.)^3
& &&&number of atoms/cell& && &=& && && && &2
& &&&number of atomic types& & =& && && && &1
& &&&number of electrons& && & =& && &&&22.00
& &&&number of Kohn-Sham states=& && && &&&16
& &&&kinetic-energy cutoff& &&&=& && &90.0000&&Ry
& &&&charge density cutoff& &&&=& &&&400.0000&&Ry&&#应该用到900吗？
& &&&Exchange-correlation& && &=&&SLA&&PW& &PBE&&PBE ( 1 4 3 4 0)
& &&&EXX-fraction& && && && &&&=& && &&&0.00
& &&&celldm(1)=& &6.254820&&celldm(2)=& &0.000000&&celldm(3)=& &0.000000
& &&&celldm(4)=& &0.000000&&celldm(5)=& &0.000000&&celldm(6)=& &0.000000
& &&&crystal axes: (cart. coord. in units of alat)
& && && && && &a(1) = (& &1.000000& &0.000000& &0.000000 )&&
& && && && && &a(2) = (&&-0.500000& &0.866025& &0.000000 )&&
& && && && && &a(3) = (& &0.000000& &0.000000& &1.592118 )&&
& &&&reciprocal axes: (cart. coord. in units 2 pi/alat)
& && && && && &b(1) = (&&1...000000 )&&
& && && && && &b(2) = (&&0...000000 )&&
& && && && && &b(3) = (&&0...628094 )&&
& &&&PseudoPot. # 1 for Sc read from file:
& &&&/home/users/spclibtt/home/yexq/shanghai-opt/pw-work/pseudo/Sc.pbe-nsp-van.UPF
& &&&MD5 check sum: ff19847cdcbe4b0e9f43ec
& &&&Pseudo is Ultrasoft + core correction, Zval = 11.0
& &&&Generated by new atomic code, or converted to UPF format
& &&&Using radial grid of&&849 points,&&6 beta functions with:
& && && && && & l(1) =& &0
& && && && && & l(2) =& &0
& && && && && & l(3) =& &1
& && && && && & l(4) =& &1
& && && && && & l(5) =& &2
& && && && && & l(6) =& &2
& &&&Q(r) pseudized with&&6 coefficients,&&rinner =& & 1.200& &1.200& &1.200
& && && && && && && && && && && && && && && && && && & 1.200& &1.200
& &&&atomic species& &valence& & mass& &&&pseudopotential
& && &&&Sc& && && && &11.00& & 44.95590& &&&Sc( 1.00)
& &&&24 Sym. Ops., with inversion, found (12 have fractional translation)
& &Cartesian axes
& &&&site n.& &&&atom& && && && && && &positions (alat units)
& && && &1& && && &&&Sc&&tau(& &1) = (& &0.0000000& &0.5773503& &0.3980295&&)
& && && &2& && && &&&Sc&&tau(& &2) = (& &0.5000000& &0.2886751& &1.1940885&&)
& &&&number of k points=&&2601
& &&&Number of k-points >= 100: set verbosity='high' to print them.
& &&&Dense&&grid:& & 45639 G-vectors& &&&FFT dimensions: (&&40,&&40,&&64)
& &&&Smooth grid:& & 38911 G-vectors& &&&FFT dimensions: (&&40,&&40,&&64)
& &&&Largest allocated arrays& &&&est. size (Mb)& &&&dimensions
& && &&&Kohn-Sham Wavefunctions& && && &0.08 Mb& &&&(& & 318,& &16)
& && &&&NL pseudopotentials& && && && & 0.17 Mb& &&&(& & 318,& &36)
& && &&&Each V/rho on FFT grid& && && & 0.10 Mb& &&&(& &6400)
& && &&&Each G-vector array& && && && & 0.02 Mb& &&&(& &2855)
& && &&&G-vector shells& && && && && &&&0.01 Mb& &&&(& &1357)
& &&&Largest temporary arrays& &&&est. size (Mb)& &&&dimensions
& && &&&Each subspace H/S matrix& && &&&0.00 Mb& &&&(&&16,&&16)
& && &&&Each
matrix& && &0.01 Mb& &&&(& &&&36,& &16)
& &&&The potential is recalculated from file :
& &&&/home/users/spclibtt/home/yexq/shanghai-opt/pw-work/environment/tmp/Sc.save/charge-density.dat
& &&&Starting wfc are& &20 randomized atomic wfcs
& &&&total cpu time spent up to now is& && &&&1.3 secs
& &&&per-process dynamical memory:& & 13.6 Mb
& &&&Band Structure Calculation& &#【启动了Band Structure Calculation ，但特征值不收敛】
& &&&CG style diagonalization
& &&&c_bands:&&1 eigenvalues not converged&&#这是什么错误
& &&&c_bands:&&1 eigenvalues not converged
& &&&c_bands:&&2 eigenvalues not converged
& &&&c_bands:&&3 eigenvalues not converged
&&.........#【略去n个...eigenvalues not converged】
& &&&c_bands:&&2 eigenvalues not converged
& &&&c_bands:&&3 eigenvalues not converged
& &&&c_bands:&&2 eigenvalues not converged
& &&&c_bands:&&2 eigenvalues not converged
& &&&c_bands:&&1 eigenvalues not converged
& &&&c_bands:&&1 eigenvalues not converged
& &&&ethr =&&4.55E-13,&&avg # of iterations = 57.3
& &&&total cpu time spent up to now is& &&&3955.9 secs
& &&&End of band structure calculation
& &&&Number of k-points >= 100: set verbosity='high' to print the bands.
& &&&highest occupied, lowest unoccupied level (ev):& &&&9.9643& & 8.6137
& &&&Writing output data file Sc.save
& &&&init_run& &&&:& && &0.50s CPU& && &0.67s WALL (& && & 1 calls)
& &&&electrons& & :& &3824.17s CPU& &3954.63s WALL (& && & 1 calls)
& &&&Called by init_run:
& &&&wfcinit& && &:& && &0.00s CPU& && &0.00s WALL (& && & 1 calls)
& &&&potinit& && &:& && &0.11s CPU& && &0.11s WALL (& && & 1 calls)
& &&&Called by electrons:
& &&&c_bands& && &:& &3824.16s CPU& &3954.62s WALL (& && & 1 calls)
& &&&v_of_rho& &&&:& && &0.08s CPU& && &0.09s WALL (& && & 1 calls)
& &&&newd& && && &:& && &0.04s CPU& && &0.04s WALL (& && & 1 calls)
& &&&Called by c_bands:
& &&&init_us_2& & :& && &1.08s CPU& && &1.03s WALL (& & 2601 calls)
& &&&ccgdiagg& &&&:& &3501.84s CPU& &3601.94s WALL (& &12559 calls)
& &&&wfcrot& && & :& & 317.80s CPU& & 319.55s WALL (& &12559 calls)
& &&&Called by *cgdiagg:
& &&&h_psi& && &&&:& &2831.31s CPU& &2893.42s WALL ( 2168039 calls)
& &&&s_psi& && &&&:& & 201.42s CPU& & 204.42s WALL ( 4323519 calls)
& &&&cdiaghg& && &:& && &4.95s CPU& && &5.16s WALL (& &12559 calls)
& &&&Called by h_psi:
& &&&add_vuspsi& &:& & 115.68s CPU& & 118.56s WALL ( 2168039 calls)
& &&&General routines
& &&&calbec& && & :& & 507.36s CPU& & 523.06s WALL ( 4323519 calls)
& &&&fft& && && & :& && &0.08s CPU& && &0.09s WALL (& && &13 calls)
& &&&ffts& && && &:& && &0.00s CPU& && &0.00s WALL (& && & 1 calls)
& &&&fftw& && && &:& &2282.87s CPU& &2296.81s WALL ( 4733656 calls)
& &&&interpolate&&:& && &0.00s CPU& && &0.00s WALL (& && & 1 calls)
& &&&davcio& && & :& && &0.14s CPU& &&&10.04s WALL (& & 2601 calls)
& &&&Parallel routines
& &&&fft_scatter&&:& &1293.12s CPU& &1310.33s WALL ( 4733670 calls)
& &&&PWSCF& && &&&:& &&&1h 3m CPU& && &&&1h 6m WALL
& &This run was terminated on:&&12:24:32&&13Mar2013& && && && &
=------------------------------------------------------------------------------=
& &JOB DONE.
=------------------------------------------------------------------------------=
可能是激发态确实不太好算吧。
vasp所用的赝势跟pwscf不一样，平衡位置肯定也不一样，vasp优化好了，在pwscf里面一般是有受力的。所以，没办法，哈哈
不好意思，插到楼主的帖子里问个问题，您好，昨天我的pwscf安装帖子您也帮助了很多，然后我还有几个问题，现在我在PC上跑的时候都是按照example在学习使用，但是里面的膺势都是人家选好的，我想知道我自己算的时候我该怎么选呢，我在网上找了下，人家都说要测试，要怎么测试呢，是不是从官网下了， 然后一个一个测试，算自洽测试么？还有，我怎么看这个膺势是LDA还是GGA呢，直接打开膺势文件看么？不好意思，我们组用的是VASP，现在导师突然让我学习pwscf，我是初学者，很多问题不懂，谢谢了
如何测赝势，可以参见castep的帖子，过程是一样的：
http://emuch.net/bbs/viewthread.php?tid=5153783
LDA还是GGA，参见此帖：
http://emuch.net/bbs/viewthread.php?tid=5464861&fpage=1
这是我仿照侯博的Cu的计算例子文件写的脚步，包括ESPRESSO 5.0.2 pp/examples/example02或老版的example08计算费米面的例子文件中非自洽计算部分均无需设定smearing。Sc是金属。
& && & Sc.pbe-sp-hgh.UPF是模守恒赝势，之前用的是超软赝势Sc.pbe-nsp-van.UPF：ecutwfc = 90,& &ecutrho = 960.00。结果同上，仍然是不收敛。
& && &赝势网址如下：
&&http://www.quantum-espresso.org/?page_id=887&el_id=21&unp_id=&fun_id=&colum_k=&origin_id=
#######################【自洽计算，产生$name.scf.out_$a文件，即Sc.scf.out_0文件】################################ self-consistent calculation -Non-Spin-Polarised case&&#【为什么要进行自洽计算，目的是.....】
for a in 0
cat > $name.scf.in_$a << EOF
& & calculation = 'scf'
& & restart_mode='from_scratch',
& & prefix='$name',
& & pseudo_dir = '$PSEUDO_DIR/',
& & outdir='$TMP_DIR/'
& & tstress=.t.,
& & tprnfor=.t.
&&ibrav=0,
&&ecutwfc=90.0, ecutrho = 360.0
&&occupations='smearing',
&&smearing='methfessel-paxton',
&&degauss=0.02
&electrons
conv_thr = 1.0d-8&&
mixing_beta = 0.7
#cell_dynamics= 'damp-pr'
# press=$a
ATOMIC_SPECIES
Sc& & 44.9559& & Sc.pbe-nsp-van.UPF
CELL_PARAMETERS {bohr}
ATOMIC_POSITIONS {crystal}
Sc 0.4 0.0000
Sc 0.6 0.0000
K_POINTS {automatic}
8 8 4 0 0 0
$PW_ROOT/pw.x&&$name.scf.out_$a
#########【准备并用kvecs_FS.x程序产生计算费米面要用到的密集网格k点，各k点坐标存储于文件kves_$Sysname中】#############
#kves_$Sysname的定义见keves_FS.f源程序
# prepare input file $name.fs.in
Sysname='$name'
Calc_Type='FS'
nabc=' 16 16&&8'& &&&#【k点分割数与一般的计算相同吗】
E_Fermi=`grep Fermi $name.scf.out_$a| cut-c 26-36`& & #【$name.scf.out_$a，由前一步自洽计算产生】
a1=`grep 'b(1)' $name.scf.out_$a| cut -c 24-54`
a2=`grep 'b(2)' $name.scf.out_$a| cut -c 24-54`
a3=`grep 'b(3)' $name.scf.out_$a| cut -c 24-54`
cat > kvecs_FS.in <<EOF
$a1& && && && && && && && && & #$a1具有xx,xx,xx的数字格式，相当于一维向量，与nabc格式类似，见上
kvecs_Fs.x kvecs_FS.out
#########################【非自洽计算，产生文件$name.nscf.out】########################
# 进行一次非自洽计算，计算上一步产生的各k点的本征值；
cat > $name.nscf.in << EOF
& &calculation='nscf',
& &prefix='$name'
& &pseudo_dir = '$PSEUDO_DIR/',
& &outdir='$TMP_DIR/'
& & ibrav = 0,
& &nat=&&2,
& &ntyp= 1,
& &ecutwfc = 90,
& &ecutrho = 400.00,
#&&occupations='smearing',& &#【非自洽计算这些参数均无必要】
&&# smearing='methfessel-paxton',
&&# degauss=0.02
& &nbnd=8&&#【一定要定义nbnd吗？默认值是.....】
&electrons
&&diagonalization='cg'&&
& &mixing_beta = 0.7
& &conv_thr =&&1.0d-10& && && & #【1.0d-10与1.0e-10含义一样吗】
#【diagonalization='cg' 英文注释：conjugate-gradient-like band-by-band diagonalization. Typically slower than 'david' but it uses less memory& &&&and is more robust (it seldom fails)】
#cell_dynamics= 'damp-pr'
# press=$a
ATOMIC_SPECIES
Sc& & 44.9559& & Sc.pbe-nsp-van.UPF
CELL_PARAMETERS {bohr}
ATOMIC_POSITIONS {crystal}
Sc 0.4 0.0000
Sc 0.6 0.0000
K_POINTS& & #【K_POINTS已在上一步产生】
cat kvecs_$Sysname >>$name.nscf.in 【将上一步产生的各k点坐标，附加到$name.nscf.in 文件后面】
$PW_ROOT/pw.x$name.nscf.out
##########【采用bands_FS.x程序将上一步得到的k点及相应本征值转化成xcrysden能读取的bxsf格式】##########################
#【下面是即将见证奇迹的时刻-采用bands_FS.x将计算的k点以及相应的本征值，转换成xcrysden软件的bxsf格式，以便采用xcrysden来画图。 】
# prepare input data (input_FS, Bands.out)for bands_FS
mv $name.nscf.out&&Bands_NSP.out
& & #【注意这利用了非自洽计算产生的$name.nscf.out文件，该步的意义可能是bands_FS.x陈旭默认读取文件名为 Bands_NSP.out?】
cat > input_FS <<EOF
$n_start $n_last
$PW_ROOT/bands_FS.x bands_fs.out& & #【这一操作将产生Bands_FS.bxsf文件】
mv Bands_FS.bxsf&&$name.fs_NSP.bxsf& & #【这一操作意义不大，是为了记忆方便】
$ECHO "&&Fermi surface plot: use 'xcrysden --bxsf $name.fs_NSP.bxsf' to plot ...\c"
$ECHO" done"
#######################################################
ESPRESSO 5.0.2 pp/examples/example02解读
您好，谢谢，也即是说测试赝势主要是看跑出来的结果和文献进行比对，相差不大即可；我昨天也找了下资料，总结了下，也就是说一般LDA包括PAW ,PZ ;GGA包括PBE、PW91、B88-P86、BLYP，然后不能混用，所有的元素要么都用PBE的，要么都用PW91的，对吧，是这样的吧？
贴出来看看，哪里不一样，我学习下。
测试赝势主要指测试其收敛性：ecutrho、degauss、k mesh。后面有一点小错误，在这里纠正一下：
PAW是与模守恒和超软相提并论的，它既可能是LDA也可能是GGA。
谢谢，我看VASP里我们的pot文件夹就是PAW，那个我学习pwscf第一个学习的也是画fermi surface，我可以画出Ni和Cu的，但是我在重复楼主的Sc的时候，用xcrysden画图的时候报错了，我现在有几个问题，1：nbnd=?到底怎么确定，看到有人回复为电子占据能带数的1.2倍以上，这个电子能带数又怎么知道，我是这样判断的，我打开自洽计算的out文件，在每一个K点下面都会有一些band的能量，我和fermi energy比较，最后相差太大的舍去，然后确定nbnd=多少；2，n_start=?n_last=？又是怎么确定，我还是看out文件，看fermi energy在那些能带附件，然后确定他们的值，不知道对不对，我这样判断；3，产生密集K点应该试差值法产生K点阿，为什么example里产生密集K点直接写16*16*16就可以了么，难道直接是把自洽计算里K的加大就可以了么，谢谢，问题有点多了，谢谢您的帮助！
1.nbnd的确定
首先电子是费米子，那么一条能带上只能呆两个电子，这个没问题吧？
然后看你的每种元素的赝势文件，里面给出一定数量的外层电子，将每个原子的外层电子相加，然后除以二，是不是就是总的有电子占据的能带数目？
其实对于绝缘体来说，能带数目取外层电子总数的一半就够了，除非你想算能隙，所以nbnd=体系外层电子数的一半
对于金属来说，最高能量的电子穿过好几条能带，它们很容易被激发，那么所选取的能带总数要大于一半的外层电子数，默认为电子数一半的1.2倍。你也可以取多一点，无妨。
2.n_start 和n_last
你的理解是对的，就是穿过费米面的能带。
3.k点的产生
我不懂什么是试差值法，故无法给出回答，抱歉！
（1）应将calculation='nscf'中的nscf改为'band'；
新版PWscf的要求。
（2）您可以尝试Ca的费米面计算，我算完后不能重现文献中的结果，非常诡异。
我不是高版本的，我是4.3.2所以应该nscf，example也是这样的，不过我遇到了新的我问题，麻烦帮我看看：
您好，我在测试ecut和K点的时候遇到了新的问题，我是这样做的，分别设ecut=30 40 50 60 70 80 90 100 110 120,然后在每个ecut值下面慢慢增加K点，但是奇怪的是，在ecut=50 60和70的时候自洽运算会出现对称性操作丢失的警告：
& & warning: symmetry operation #&&3 not allowed.& &fractional translation:
& && & 0...5000000&&in crystal coordinates
& &&&warning: symmetry operation #&&4 not allowed.& &fractional translation:
& && &-0...5000000&&in crystal coordinates
& &&&warning: symmetry operation #&&7 not allowed.& &fractional translation:
& && &-0...5000000&&in crystal coordinates
& &&&warning: symmetry operation #&&8 not allowed.& &fractional translation:
& && & 0...5000000&&in crystal coordinates
& &&&warning: symmetry operation # 11 not allowed.& &fractional translation:
& && &-0...5000000&&in crystal coordinates
& &&&warning: symmetry operation # 12 not allowed.& &fractional translation:
& && & 0...5000000&&in crystal coordinates
& &&&warning: symmetry operation # 15 not allowed.& &fractional translation:
& && & 0...5000000&&in crystal coordinates
& &&&warning: symmetry operation # 16 not allowed.& &fractional translation:
& && &-0...5000000&&in crystal coordinates
我看到小木虫上有说通过增加ecut可以解决，可是为什么30 和40又能跑呢？这个问题怎么解决，我同时贴出我的IN文件：
& && && && && &&&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' ,
& && && && && && && &&&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.& & 0.& & 0.
& &&&0.& & 8.& & 0.
& &&&0.& & 0.& & 5.
ATOMIC_SPECIES
& &Cr& &52.00000&&Cr.pz-hgh.UPF
& & O& &16.00000&&O.pz-hgh.UPF
ATOMIC_POSITIONS crystal
& &Cr& && &0.& & 0.& & 0.& &
& &Cr& && &0.& & 0.& & 0.& &
& & O& && &0.& & 0.& & 0.& &
& & O& && &0.& & 0.& & 0.& &
& & O& && &0.& & 0.& & 0.& &
& & O& && &0.& & 0.& & 0.& &
K_POINTS automatic
&&6 6 9& &0 0 0
2，第二个问题是：我发现计算结果的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&&)
麻烦帮我看下，谢谢！
研究生必备与500万研究生在线互动！
扫描下载送金币CF CF外挂源码请大家自己研究使用谢谢呵呵 过夏非法 Game Hook Crack 编程 238万源代码下载-
&文件名称: CF
& & & & &&]
&&所属分类:
&&开发工具: EasyLanguage易语言
&&文件大小: 3617 KB
&&上传时间:
&&下载次数: 236
&&提 供 者:
&详细说明：CF外挂源码请大家自己研究使用谢谢呵呵自己过夏非法-CFwaiguayuanmaqingdajzijyanjyanj
文件列表(点击判断是否您需要的文件，如果是垃圾请在下面评价投诉):
&&恶魔天尊网游模块1.0.ec&&窗口化游戏.ec&&单一实例.ec&&窗口句柄查询.ec&&Super-EC.ec&&枫式热键FS1.2.ec&&超级伪装者2.0.exe&&测试拿房.bmp&&皮肤模块.ec&&远洋外挂源码.e
&[]:纯粹是垃圾&[]:纯粹是垃圾&[]:很好，推荐下载&[]:很好，推荐下载
&近期下载过的用户:
&相关搜索:
&&&&&&&&&&
&输入关键字，在本站238万海量源码库中尽情搜索：
&[] - QQ超级农民农场外挂源码
需要的朋友用下
&[] - 做游戏外挂必备，快速搜寻角色各项相关信息。此为易语言源码版。
&[] - CF透视挂源码。。。。Delphi的.....
&[] - 穿越火线辅助器月9号编辑穿越火线外挂易语言源码
&[] - CF早出3秒,源.
比较简单的一个源码
&[] - 易语言编写的一个非常简单的加密解密程序,可以用在网络通信加密里面,速度比较快,可以根据需要修改一下算法,实现更好的加密
&[] - 游戏外挂客户端服务端充值系统（开源）`学习研究
&[] - 一篇很好的有关提取关键帧的论文，介绍了多种提取方法，值得一看。
&[] - 武林外传脱机外挂，易语言全部源码，可后台挂机
&[] - Autoit Script 写的跑跑卡丁车隐身外挂