【求助】求助 castep计算
我刚接触MS计算,用castep中进行几何优化和能带计算时,出现错误,最后给出:
Warning: electronic minimisation did not converge when finding ground state.
不知道要调试哪个参数,请高手指导,谢谢。
如能根据文件信息给出相应的指导解释信息,则不甚感激。
具体文档为
Pseudo atomic calculation performed for B 2s2 2p1
Converged in 16 iterations to a total energy of -70.4905 eV
Pseudo atomic calculation performed for O 2s2 2p4
Converged in 25 iterations to a total energy of -429.5309 eV
Pseudo atomic calculation performed for Y 4d1 5s2
Converged in 17 iterations to a total energy of -188.5056 eV
Calculation parallelised over 1 nodes.
K-points are distributed over 1 groups, each containing 1 nodes.
************************************ Title ************************************
CASTEP calculation from Materials Studio
***************************** General Parameters ******************************
output verbosity : normal (1)
write checkpoint data to : YBO3.check
type of calculation : geometry optimization
stress calculation : on
density difference calculation : off
electron localisation func (ELF) calculation : off
unlimited duration calculation
timing information : on
memory usage estimate : on
write final potential to formatted file : off
write final density to formatted file : off
output length unit : A
output mass unit : amu
output time unit : ps
output charge unit : e
output energy unit : eV
output force unit : eV/A
output velocity unit : A/ps
output pressure unit : GPa
output inv_length unit : 1/A
output frequency unit : cm-1
output force constant unit : eV/A**2
output volume unit : A**3
output IR intensity unit : (D/A)**2/amu
output dipole unit : D
output efield unit : eV/A/e
wavefunctions paging : none
random number generator seed : randomised (191520921)
data distribution : optimal for this architecture
optimization strategy : balance speed and memory
*********************** Exchange-Correlation Parameters ***********************
using functional : Perdew Burke Ernzerhof
Divergence correction : off
************************* Pseudopotential Parameters **************************
pseudopotential representation : reciprocal space
<beta|phi> representation : reciprocal space
**************************** Basis Set Parameters *****************************
basis set accuracy : FINE
plane wave basis set cut-off : 340.0000 eV
size of standard grid : 2.0000
largest prime factor in FFT : 5
finite basis set correction : automatic
number of sample energies : 3
sample spacing : 5.0000 eV
**************************** Electronic Parameters ****************************
number of electrons : 135.0
net charge of system : 0.000
net spin of system : 1.000
number of up spins : 68.00
number of down spins : 67.00
treating system as spin-polarized
number of bands : 68
********************* Electronic Minimization Parameters **********************
Method: Treating system as non-metallic,
and number of SD steps : 1
and number of CG steps : 4
total energy / atom convergence tol. : 0.1000E-05 eV
eigen-energy convergence tolerance : 0.4265E-06 eV
max force / atom convergence tol. : ignored
convergence tolerance window : 3 cycles
max. number of SCF cycles : 100
************************** Density Mixing Parameters **************************
density-mixing scheme : Pulay
max. length of mixing history : 20
charge density mixing amplitude : 0.5000
spin density mixing amplitude : 2.000
charge density mixing g-vector : 1.500 1/A
spin density mixing g-vector : 1.500 1/A
********************** Geometry Optimization Parameters ***********************
optimization method : BFGS
variable cell method : fixed basis quality
max. number of steps : 100
estimated bulk modulus : 500.0 GPa
estimated <frequency> : 1668. cm-1
line minimiser tolerance : 0.4000
with spin fully able to relax for all steps
total energy convergence tolerance : 0.1000E-04 eV/atom
max ionic |force| tolerance : 0.3000E-01 eV/A
max ionic |displacement| tolerance : 0.1000E-02 A
max |stress component| tolerance : 0.5000E-01 GPa
convergence tolerance window : 2 steps
backup results every : 5 steps
*******************************************************************************
-------------------------------
Unit Cell
-------------------------------
Real Lattice(A) Reciprocal Lattice(1/A)
3.2701119 -1.8880000 0.0000000 1.9213976 0.0000000 0.0000000
0.0000000 3.7760000 0.0000000 0.9606988 1.6639791 0.0000000
0.0000000 0.0000000 8.8059998 0.0000000 0.0000000 0.7135119
Lattice parameters(A) Cell Angles
a = 3.776000 alpha = 90.000000
b = 3.776000 beta = 90.000000
c = 8.806000 gamma = 120.000000
Current cell volume = 108.735981 A**3
-------------------------------
Cell Contents
-------------------------------
Total number of ions in cell = 29
Total number of species in cell = 3
Max number of any one species = 16
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
x Element Atom Fractional coordinates of atoms x
x Number u v w x
x----------------------------------------------------------x
x B 1 0.586000 0.435000 0.250000 x
x B 2 -0.435000 0.151000 0.250000 x
x B 3 -0.151000 -0.586000 0.250000 x
x B 4 -0.586000 -0.435000 0.250000 x
x B 5 0.435000 -0.151000 0.250000 x
x B 6 0.151000 0.586000 0.250000 x
x B 7 -0.586000 -0.435000 -0.250000 x
x B 8 0.435000 -0.151000 -0.250000 x
x B 9 0.151000 0.586000 -0.250000 x
x B 10 0.586000 0.435000 -0.250000 x
x B 11 -0.435000 0.151000 -0.250000 x
x B 12 -0.151000 -0.586000 -0.250000 x
x O 1 0.774000 0.881000 0.250000 x
x O 2 -0.881000 -0.107000 0.250000 x
x O 3 0.107000 -0.774000 0.250000 x
x O 4 -0.774000 -0.881000 0.250000 x
x O 5 0.881000 0.107000 0.250000 x
x O 6 -0.107000 0.774000 0.250000 x
x O 7 -0.774000 -0.881000 -0.250000 x
x O 8 0.881000 0.107000 -0.250000 x
x O 9 -0.107000 0.774000 -0.250000 x
x O 10 0.774000 0.881000 -0.250000 x
x O 11 -0.881000 -0.107000 -0.250000 x
x O 12 0.107000 -0.774000 -0.250000 x
x O 13 0.666667 0.333333 0.110300 x
x O 14 -0.666667 -0.333333 0.110300 x
x O 15 -0.666667 -0.333333 -0.110300 x
x O 16 0.666667 0.333333 -0.110300 x
x Y 1 0.000000 0.000000 0.000000 x
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
No user defined ionic velocities
-------------------------------
Details of Species
-------------------------------
Mass of species in AMU
B 10.8100004
O 15.9989996
Y 88.9059982
Electric Quadrupole Moment (Barn)
B 0.0405900 Isotope 11
O -0.0255800 Isotope 17
Y 1.0000000 No Isotope Defined
Files used for pseudopotentials:
B B_00PBE.usp
O O_00PBE.usp
Y Y_00PBE.uspcc
-------------------------------
k-Points For BZ Sampling
-------------------------------
MP grid size for SCF calculation is 8 8 4
Number of kpoints used = 24
返回小木虫查看更多
京公网安备 11010802022153号
把截断能量值改了
但是把截断能改小,会影响结果啊
恩,对你优化好之后的几何构型再次进行优化,应该就可以了,而且时间应该会很短。或者你降低收敛标准试试看,