24小时热门版块排行榜    

查看: 3999  |  回复: 7

li_sen

新虫 (小有名气)

[求助] castep结构优化报错

我在用CASTEP做结构优化的时候,刚开始一分钟就报错,这个CIF文件是CM发表出来的,感觉应该没有问题,究竟是何原因,请高手指点。我也计算了其他的物质,软件没有问题。
报错信息:
Time Module Type Message Text
2018/09/07@14:26:58  Diagnostics INFO Log file created
2018/09/07@14:28:44  [0AK3O] - cm7b00702_si_002 CASTEP GeomOpt ERROR CASTEP terminated with unknown error status.
2018/09/07@14:29:09  [0AK3O] - cm7b00702_si_002 CASTEP GeomOpt ERROR Error on Job Complete. Failed to import document from C:\Date of software\Matrrials studio\Calculation date\Bi13S18I2-8-8\Bi13S18I2-8-8_Files\Documents\cm7b00702_si_002 CASTEP GeomOpt\cm7b00702_si_002.xsd. - Binary file with CASTEP results not present.
2018/09/07@14:29:10  [0AK3O] - cm7b00702_si_002 CASTEP GeomOpt ERROR Files not removed from server due to previous errors.
2018/09/07@14:29:10  [0AK3O] - cm7b00702_si_002 CASTEP GeomOpt WARN As instructed, files have been left on remote server, these may need to be manually archived.
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mhslee

木虫 (正式写手)

李明憲
2楼2018-09-07 18:49:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

li_sen

新虫 (小有名气)

引用回帖:
2楼: Originally posted by mhslee at 2018-09-07 18:49:59
要查看 .castep 文件。

Job started on host T50GDMG
at Sat Sep  8 14:57:41 2018

+-------------------------------------------------+
|                                                 |
|      CCC   AA    SSS  TTTTT  EEEEE  PPPP        |
|     C     A  A  S       T    E      P   P       |
|     C     AAAA   SS     T    EEE    PPPP        |
|     C     A  A     S    T    E      P           |
|      CCC  A  A  SSS     T    EEEEE  P           |
|                                                 |
+-------------------------------------------------+
|                                                 |
| Welcome to Materials Studio CASTEP version 8.0  |
| Ab Initio Total Energy Program                  |
|                                                 |
| Authors:                                        |
| M. Segall, M. Probert, C. Pickard, P. Hasnip,   |
| S. Clark, K. Refson, J. R. Yates, M. Payne      |
|                                                 |
| Contributors:                                   |
| P. Lindan, P. Haynes, J. White, V. Milman,      |
| N. Govind, M. Gibson, P. Tulip, V. Cocula,      |
| B. Montanari, D. Quigley, M. Glover,            |
| L. Bernasconi, A. Perlov, M. Plummer,           |
| E. McNellis, J. Meyer, J. Gale, D. Jochym       |
| J. Aarons, B. Walker, R. Gillen, D. Jones       |
| T. Green                                        |
|                                                 |
| Copyright (c) 2000 - 2014                       |
|                                                 |
| Please cite                                     |
|                                                 |
|     "First principles methods using CASTEP"     |
|                                                 |
|         Zeitschrift fuer Kristallographie       |
|           220(5-6) pp. 567-570 (2005)           |
|                                                 |
| S. J. Clark, M. D. Segall, C. J. Pickard,       |
| P. J. Hasnip, M. J. Probert, K. Refson,         |
| M. C. Payne                                     |
|                                                 |
|       in all publications arising from          |
|              your use of CASTEP                 |
|                                                 |
+-------------------------------------------------+


This version was compiled for x86_64-windows-msvc2008 on Dec 04 2014
Code version:    6546
Intel(R) Math Kernel Library Version 11.1.2
Fundamental constants values: CODATA 2010

License checkout of MS_castep successful


Pseudo atomic calculation performed for S 3s2 3p4

Converged in 20 iterations to a total energy of -273.8401 eV


Pseudo atomic calculation performed for I 5s2 5p5

Converged in 18 iterations to a total energy of -311.6166 eV


Pseudo atomic calculation performed for Bi 6s2 6p3

Converged in 18 iterations to a total energy of -147.1723 eV

Calculation parallelised over 10 processes.
Data is distributed by k-point(10-way)

************************************ Title ************************************
CASTEP calculation from Materials Studio

***************************** General Parameters ******************************
  
output verbosity                               : normal  (1)
write checkpoint data to                       : cm7b00702_si_002.check
type of calculation                            : geometry optimization
stress calculation                             : on
density difference calculation                 : off
electron localisation func (ELF) calculation   : off
Hirshfeld analysis                             : on
unlimited duration calculation
timing information                             : on
memory usage estimate                          : on
write final potential to formatted file        : off
write final density to formatted file          : off
write BibTeX reference list                    : on
checkpoint writing                             : both castep_bin and check files
  
output         length unit                     : A
output           mass unit                     : amu
output           time unit                     : ps
output         charge unit                     : e
output           spin unit                     : hbar/2
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
output        entropy unit                     : J/mol/K
  
wavefunctions paging                           : none
random number generator seed                   : randomised (145743955)
data distribution                              : optimal for this architecture
optimization strategy                          : balance speed and memory

*********************** Exchange-Correlation Parameters ***********************
  
using functional                               : Perdew Burke Ernzerhof
Divergence correction                          : off
relativistic treatment                         : Koelling-Harmon
DFT+D: Semi-empirical dispersion correction    : off

************************* Pseudopotential Parameters **************************
  
pseudopotential representation                 : reciprocal space
<beta|phi> representation                      : reciprocal space

**************************** Basis Set Parameters *****************************
  
plane wave basis set cut-off                   :   250.0000   eV
size of standard grid                          :     2.0000
size of   fine   grid                          :     3.0000
size of   fine   gmax                          :    24.3013   1/A
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                           :  202.0   
net charge of system                           :  0.000   
net spin   of system                           :  0.000   
number of  up  spins                           :  101.0   
number of down spins                           :  101.0   
treating system as non-spin-polarized
number of bands                                :        122

********************* Electronic Minimization Parameters **********************
  
Method: Treating system as metallic with density mixing treatment of electrons,
         and number of  SD  steps               :          1
         and number of  CG  steps               :          4
  
total energy / atom convergence tol.           : 0.5000E-06   eV
eigen-energy convergence tolerance             : 0.1475E-06   eV
max force / atom convergence tol.              : ignored
convergence tolerance window                   :          3   cycles
max. number of SCF cycles                      :        100
number of fixed-spin iterations                :         10
smearing scheme                                : Gaussian
smearing width                                 : 0.1000       eV
Fermi energy convergence tolerance             : 0.2721E-13   eV
periodic dipole correction                     : NONE

************************** Density Mixing Parameters **************************
  
density-mixing scheme                          : Pulay
max. length of mixing history                  :         20
charge density mixing amplitude                : 0.5000   
cut-off energy for mixing                      :  250.0       eV
charge density mixing g-vector                 :  1.500       1/A

*********************** Population Analysis Parameters ************************
  
Population analysis with cutoff                :  3.000       A

********************** Geometry Optimization Parameters ***********************
  
optimization method                            : BFGS
variable cell method                           : fixed basis quality
max. number of steps                           :        500
estimated bulk modulus                         :  500.0       GPa
estimated <frequency>                          :  1668.       cm-1
geom line minimiser                            : on
with line minimiser tolerance                  :     0.4000
total energy convergence tolerance             : 0.5000E-05   eV/atom
        max ionic |force| tolerance             : 0.1000E-01   eV/A
max ionic |displacement| tolerance             : 0.5000E-03   A
   max |stress component| tolerance             : 0.2000E-01   GPa
convergence tolerance window                   :          2   steps
backup results every                           :          5   steps

*******************************************************************************
  

                           -------------------------------
                                      Unit Cell
                           -------------------------------
        Real Lattice(A)                      Reciprocal Lattice(1/A)
  13.5203886  -7.8060000   0.0000000        0.4647193   0.0000000   0.0000000
   0.0000000  15.6120000   0.0000000        0.2323596   0.4024587   0.0000000
   0.0000000   0.0000000   4.0168000        0.0000000   0.0000000   1.5642266

                       Lattice parameters(A)       Cell Angles
                    a =   15.612000          alpha =   90.000000
                    b =   15.612000          beta  =   90.000000
                    c =    4.016800          gamma =  120.000000

                       Current cell volume =  847.867377       A**3

                           -------------------------------
                                     Cell Contents
                           -------------------------------
                         Total number of ions in cell =   36
                      Total number of species in cell =    3
                        Max number of any one species =   18

            xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
            x  Element    Atom        Fractional coordinates of atoms  x
            x            Number           u          v          w      x
            x----------------------------------------------------------x
            x  S            1         0.201700   0.271800   0.278900   x
            x  S            2         0.148300   0.499800   1.277300   x
            x  S            3         0.316400   0.480500   0.777000   x
            x  S            4         0.937800   0.403400  -0.221300   x
            x  S            5         0.721000   0.120600   0.777600   x
            x  S            6         0.946400   0.212600   0.278000   x
            x  S            7        -0.271800  -0.070100   0.278900   x
            x  S            8        -0.499800  -0.351500   1.277300   x
            x  S            9        -0.480500  -0.164100   0.777000   x
            x  S           10        -0.403400   0.534400  -0.221300   x
            x  S           11        -0.120600   0.600400   0.777600   x
            x  S           12        -0.212600   0.733800   0.278000   x
            x  S           13         0.070100  -0.201700   0.278900   x
            x  S           14         0.351500  -0.148300   1.277300   x
            x  S           15         0.164100  -0.316400   0.777000   x
            x  S           16        -0.534400  -0.937800  -0.221300   x
            x  S           17        -0.600400  -0.721000   0.777600   x
            x  S           18        -0.733800  -0.946400   0.278000   x
            x  I            1         0.666667   0.333333   0.778500   x
            x  I            2         0.000000   0.000000   0.278600   x
            x  Bi           1         0.822130   0.277330   0.277700   x
            x  Bi           2         0.844540   0.056400   0.777700   x
            x  Bi           3         0.092480   0.370340   0.777900   x
            x  Bi           4         0.388050   0.424410   0.277900   x
            x  Bi           5        -0.277330   0.544800   0.277700   x
            x  Bi           6        -0.056400   0.788140   0.777700   x
            x  Bi           7        -0.370340  -0.277860   0.777900   x
            x  Bi           8        -0.424410  -0.036360   0.277900   x
            x  Bi           9        -0.544800  -0.822130   0.277700   x
            x  Bi          10        -0.788140  -0.844540   0.777700   x
            x  Bi          11         0.277860  -0.092480   0.777900   x
            x  Bi          12         0.036360  -0.388050   0.277900   x
            x  Bi          13         0.333333   0.666667   0.921000   x
            x  Bi          14         0.333333   0.666667   0.135000   x
            x  Bi          15         0.333333   0.666667   0.421000   x
            x  Bi          16         0.333333   0.666667   0.636000   x
            xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx


                         No user defined ionic velocities

                           -------------------------------
                                   Details of Species
                           -------------------------------

                               Mass of species in AMU
                                    S   32.0600014
                                    I  126.9039993
                                    Bi  208.9799957

                          Electric Quadrupole Moment (Barn)
                                    S   -0.0678000 Isotope 33
                                    I   -0.6960000 Isotope127
                                    Bi   -0.5160000 Isotope209

                          Files used for pseudopotentials:
                                    S S_00PBE.usp
                                    I I_00PBE.usp
                                    Bi Bi_00PBE.usp

                           -------------------------------
                              k-Points For BZ Sampling
                           -------------------------------
                       MP grid size for SCF calculation is  6  6  6
                            with an offset of   0.083  0.083  0.000
                       Number of kpoints used =            42

                           -------------------------------
                               Symmetry and Constraints
                           -------------------------------

                      Maximum deviation from symmetry = 0.181994E-13     ANG

                      Number of symmetry operations   =           3
         There are no ionic constraints specified or generated for this cell
                      Point group of crystal =    16: C3, 3, 3
                      Space group of crystal =   143: P3, P 3

             Set iprint > 1 for details on symmetry rotations/translations

                         Centre of mass is NOT constrained

                         Number of cell constraints= 4
                         Cell constraints are:  1 1 3 0 0 0

                         External pressure/stress (GPa)
                          0.00000   0.00000   0.00000
                                    0.00000   0.00000
                                              0.00000
  
+---------------- MEMORY AND SCRATCH DISK ESTIMATES PER PROCESS --------------+
|                                                     Memory          Disk    |
| Model and support data                              203.3 MB       425.0 MB |
| Electronic energy minimisation requirements         193.3 MB         0.0 MB |
|                                               ----------------------------- |
| Approx. total storage required per process          396.6 MB       425.0 MB |
|                                                                             |
| Requirements will fluctuate during execution and may exceed these estimates |
+-----------------------------------------------------------------------------+
Calculating finite basis set correction with  3 cut-off energies.
Calculating total energy with cut-off of  240.000eV.
------------------------------------------------------------------------ <-- SCF
SCF loop      Energy           Fermi           Energy gain       Timer   <-- SCF
                               energy          per atom          (sec)   <-- SCF
------------------------------------------------------------------------ <-- SCF
Initial   4.15208602E+003  0.00000000E+000                        10.50  <-- SCF
Error, norm_sq is negative in wave_Snormalise_slice_b
Current trace stack:
wave_Snormalise_slice_b
wave_Sorthonormalise_slice
hamiltonian_diagonalise_ks
electronic_minimisation
calculate_finite_basis_corr
check_elec_ground_state
castep
Error, norm_sq is negative in wave_Snormalise_slice_b
Current trace stack:
wave_Snormalise_slice_b
wave_Sorthonormalise_slice
hamiltonian_diagonalise_ks
electronic_minimisation
calculate_finite_basis_corr
check_elec_ground_state
castep
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 1

你好,请您帮忙看一下
3楼2018-09-08 14:06:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mhslee

木虫 (正式写手)

加大平面波截止動能看看。

发自小木虫Android客户端

» 本帖已获得的红花(最新10朵)

李明憲
4楼2018-09-09 19:53:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

li_sen

新虫 (小有名气)

送红花一朵
引用回帖:
4楼: Originally posted by mhslee at 2018-09-09 19:53:15
加大平面波截止動能看看。

我增加到了350 eV, 还是不行
5楼2018-09-10 16:12:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mhslee

木虫 (正式写手)

李明憲
6楼2018-09-12 18:31:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

li_sen

新虫 (小有名气)

送红花一朵
引用回帖:
4楼: Originally posted by mhslee at 2018-09-09 19:53:15
加大平面波截止動能看看。

谢谢您的帮助,更换赝势后,可以运行了!以后还望多多指教

发自小木虫IOS客户端
7楼2018-09-13 13:04:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mhslee

木虫 (正式写手)

引用回帖:
7楼: Originally posted by li_sen at 2018-09-13 13:04:05
谢谢您的帮助,更换赝势后,可以运行了!以后还望多多指教
...

可否提供我結構模型,我有空時對原贗勢作進一步測試?謝謝。

发自小木虫Android客户端
李明憲
8楼2018-09-16 09:44:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 li_sen 的主题更新
信息提示
请填处理意见