当前位置: 首页 > 第一原理 >【求助】求助 castep计算

【求助】求助 castep计算

作者 cwp0918
来源: 小木虫 500 10 举报帖子
+关注

我刚接触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 返回小木虫查看更多

今日热帖
  • 精华评论
  • puchunying

    把截断能量值改了

  • cwp0918

    但是把截断能改小,会影响结果啊

  • hzy900211

    恩,对你优化好之后的几何构型再次进行优化,应该就可以了,而且时间应该会很短。或者你降低收敛标准试试看,

猜你喜欢
下载小木虫APP
与700万科研达人随时交流
  • 二维码
  • IOS
  • 安卓