计算化学公社

标题: castep表面模型结构优化问题 只迭代了一步就结束了 [打印本页]
作者
Author: lzc    时间: 2015-9-12 16:23
标题: castep表面模型结构优化问题 只迭代了一步就结束了
各位老师好，我使用castep进行结构优化对Ag（111）面进行优化时，（6层结构，最下两层固定，提交任务时选择no，即选择非对称优化）结构如图： (, 下载次数 Times of downloads: 65) 请问为什么我的只优化了一部就结束了呢？能量值上面显示只有一步， (, 下载次数 Times of downloads: 66) 而且各层之间也没有发生变化，下面是输出的castep结果。提示是说优化成功了。麻烦老师帮我看一下，我出错在哪里？应该如何修改，谢谢谢谢 Job started on host node01
at Thu Sep 10 18:26:00 2015

+-------------------------------------------------+
|                                                 |
|      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-rhel5-intel14.0 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 Ag 4d10 5s1

Converged in 32 iterations to a total energy of -1024.8999 eV

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

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

***************************** General Parameters ******************************

output verbosity                               : normal  (1)
write checkpoint data to                       : AgGGA__1_1_1___2_.check
type of calculation                            : geometry optimization
stress calculation                             : off
density difference calculation                 : off
electron localisation func (ELF) calculation   : off
Hirshfeld analysis                             : 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
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 (182600939)
data distribution                              : optimal for this architecture
optimization strategy                          : balance speed and memory

*********************** Exchange-Correlation Parameters ***********************

using functional                               : Perdew Wang (1991)
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                   :   380.0000   eV
size of standard grid                          :     1.5000
size of   fine   gmax                          :    14.9803   1/A
largest prime factor in FFT                    :          5
finite basis set correction                    : none

**************************** Electronic Parameters ****************************

number of  electrons                           :  594.0   
net charge of system                           :  0.000   
net spin   of system                           :  0.000   
number of  up  spins                           :  297.0   
number of down spins                           :  297.0   
treating system as non-spin-polarized
number of bands                                :        357

********************* 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.1000E-05   eV
eigen-energy convergence tolerance             : 0.1513E-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                      :  380.0       eV
charge density mixing g-vector                 :  1.500       1/A

********************** Geometry Optimization Parameters ***********************

optimization method                            : BFGS
variable cell method                           : fixed basis quality
max. number of steps                           :        200
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.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)
   7.6127983  -4.3952512   0.0000000        0.8253450  -0.0000000  -0.0000000
   0.0000000   8.7905023  -0.0000000        0.4126725   0.7147698  -0.0000000
   0.0000000   0.0000000  26.1153584       -0.0000000   0.0000000   0.2405935

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

                       Current cell volume = 1747.648181       A**3

                           -------------------------------
                                     Cell Contents
                           -------------------------------
                         Total number of ions in cell =   54
                      Total number of species in cell =    1
                        Max number of any one species =   54

            xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
            x  Element    Atom        Fractional coordinates of atoms  x
            x            Number           u          v          w      x
            x----------------------------------------------------------x
            x  Ag           1         0.222222   0.111111   0.000000   x
            x  Ag           2         0.111111   0.222222   0.091612   x
            x  Ag           3        -0.000000  -0.000000   0.183223   x
            x  Ag           4         0.222222   0.111111   0.274835   x
            x  Ag           5         0.111111   0.222222   0.366447   x
            x  Ag           6         0.000000   0.000000   0.458058   x
            x  Ag           7         0.555556   0.111111   0.000000   x
            x  Ag           8         0.444444   0.222222   0.091612   x
            x  Ag           9         0.333333  -0.000000   0.183223   x
            x  Ag          10         0.555556   0.111111   0.274835   x
            x  Ag          11         0.444444   0.222222   0.366447   x
            x  Ag          12         0.333333   0.000000   0.458058   x
            x  Ag          13         0.888889   0.111111   0.000000   x
            x  Ag          14         0.777778   0.222222   0.091612   x
            x  Ag          15         0.666667   0.000000   0.183223   x
            x  Ag          16         0.888889   0.111111   0.274835   x
            x  Ag          17         0.777778   0.222222   0.366447   x
            x  Ag          18         0.666667   0.000000   0.458058   x
            x  Ag          19         0.222222   0.444444   0.000000   x
            x  Ag          20         0.111111   0.555556   0.091612   x
            x  Ag          21        -0.000000   0.333333   0.183223   x
            x  Ag          22         0.222222   0.444444   0.274835   x
            x  Ag          23         0.111111   0.555556   0.366447   x
            x  Ag          24         0.000000   0.333333   0.458058   x
            x  Ag          25         0.555556   0.444444   0.000000   x
            x  Ag          26         0.444444   0.555556   0.091612   x
            x  Ag          27         0.333333   0.333333   0.183223   x
            x  Ag          28         0.555556   0.444444   0.274835   x
            x  Ag          29         0.444444   0.555556   0.366447   x
            x  Ag          30         0.333333   0.333333   0.458058   x
            x  Ag          31         0.888889   0.444444   0.000000   x
            x  Ag          32         0.777778   0.555556   0.091612   x
            x  Ag          33         0.666667   0.333333   0.183223   x
            x  Ag          34         0.888889   0.444444   0.274835   x
            x  Ag          35         0.777778   0.555556   0.366447   x
            x  Ag          36         0.666667   0.333333   0.458058   x
            x  Ag          37         0.222222   0.777778   0.000000   x
            x  Ag          38         0.111111   0.888889   0.091612   x
            x  Ag          39        -0.000000   0.666667   0.183223   x
            x  Ag          40         0.222222   0.777778   0.274835   x
            x  Ag          41         0.111111   0.888889   0.366447   x
            x  Ag          42         0.000000   0.666667   0.458058   x
            x  Ag          43         0.555556   0.777778   0.000000   x
            x  Ag          44         0.444444   0.888889   0.091612   x
            x  Ag          45         0.333333   0.666667   0.183223   x
            x  Ag          46         0.555556   0.777778   0.274835   x
            x  Ag          47         0.444444   0.888889   0.366447   x
            x  Ag          48         0.333333   0.666667   0.458058   x
            x  Ag          49         0.888889   0.777778   0.000000   x
            x  Ag          50         0.777778   0.888889   0.091612   x
            x  Ag          51         0.666667   0.666667   0.183223   x
            x  Ag          52         0.888889   0.777778   0.274835   x
            x  Ag          53         0.777778   0.888889   0.366447   x
            x  Ag          54         0.666667   0.666667   0.458058   x
            xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx


                         No user defined ionic velocities

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

                               Mass of species in AMU
                                    Ag  107.8679962

                          Electric Quadrupole Moment (Barn)
                                    Ag    1.0000000 No Isotope Defined

                          Files used for pseudopotentials:
                                    Ag Ag_00PBE.usp

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

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

                      Cell is a supercell containing 9 primitive cells
                      Maximum deviation from symmetry =  0.00000         ANG

         There are no symmetry operations specified or generated for this cell
                      Number of ionic constraints     =          54
                      Point group of crystal =     1: C1, 1, 1
                      Space group of crystal =   164: P-3m1, -P 3 2"


                         Centre of mass is NOT constrained
             Set iprint > 1 for details of linear ionic constraints

                         Number of cell constraints= 6
                         Cell constraints are:  0 0 0 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                              867.6 MB       150.7 MB |
| Electronic energy minimisation requirements        1173.7 MB         0.0 MB |
|                                               ----------------------------- |
| Approx. total storage required per process         2041.3 MB       150.7 MB |
|                                                                             |
| Requirements will fluctuate during execution and may exceed these estimates |
+-----------------------------------------------------------------------------+
------------------------------------------------------------------------ <-- SCF
SCF loop      Energy           Fermi           Energy gain       Timer   <-- SCF
                               energy          per atom          (sec)   <-- SCF
------------------------------------------------------------------------ <-- SCF
Initial  -1.05486933E+004  0.00000000E+000                        19.51  <-- SCF
      1  -5.24873051E+004 -1.90733536E+000   7.76640959E+002     372.75  <-- SCF
      2  -5.53050034E+004 -7.53699450E+000   5.21795983E+001     674.71  <-- SCF
      3  -5.55423071E+004 -6.65380676E+000   4.39451316E+000     930.81  <-- SCF
      4  -5.55455166E+004 -6.64236598E+000   5.94352070E-002    1176.35  <-- SCF
      5  -5.55330970E+004 -7.01058113E+000  -2.29993888E-001    1437.44  <-- SCF
      6  -5.55329972E+004 -6.99465611E+000  -1.84687932E-003    1736.83  <-- SCF
      7  -5.55330167E+004 -6.99198927E+000   3.61612640E-004    1993.65  <-- SCF
      8  -5.55330080E+004 -6.98358368E+000  -1.61795532E-004    2254.83  <-- SCF
      9  -5.55329997E+004 -6.97057082E+000  -1.54679370E-004    2482.19  <-- SCF
     10  -5.55329761E+004 -6.91745678E+000  -4.36070562E-004    2736.23  <-- SCF
     11  -5.55329810E+004 -6.91034215E+000   8.97093823E-005    2967.24  <-- SCF
     12  -5.55329753E+004 -6.88023009E+000  -1.05269193E-004    3201.62  <-- SCF
     13  -5.55329702E+004 -6.85392965E+000  -9.44344457E-005    3427.39  <-- SCF
     14  -5.55329693E+004 -6.82216103E+000  -1.53214510E-005    3648.00  <-- SCF
     15  -5.55329836E+004 -6.81061576E+000   2.63948796E-004    3875.94  <-- SCF
     16  -5.55329901E+004 -6.79344966E+000   1.19950234E-004    4078.55  <-- SCF
     17  -5.55329936E+004 -6.78816752E+000   6.45324590E-005    4276.89  <-- SCF
     18  -5.55329938E+004 -6.79033674E+000   4.68553161E-006    4465.78  <-- SCF
     19  -5.55329938E+004 -6.78925382E+000   2.92318672E-007    4645.07  <-- SCF
     20  -5.55329938E+004 -6.78898424E+000   5.55688803E-008    4822.06  <-- SCF
------------------------------------------------------------------------ <-- SCF

Final energy, E             =  -55532.87693541     eV
Final free energy (E-TS)    =  -55532.99382952     eV
(energies not corrected for finite basis set)

NB est. 0K energy (E-0.5TS)      =  -55532.93538246     eV


Writing analysis data to AgGGA__1_1_1___2_.castep_bin

Writing model to AgGGA__1_1_1___2_.check

+---------------- MEMORY AND SCRATCH DISK ESTIMATES PER PROCESS --------------+
|                                                     Memory          Disk    |
| Model and support data                              867.6 MB       150.7 MB |
| Electronic energy minimisation requirements        1173.7 MB         0.0 MB |
| Geometry minimisation requirements                 1667.7 MB         0.0 MB |
|                                               ----------------------------- |
| Approx. total storage required per process         3709.0 MB       150.7 MB |
|                                                                             |
| Requirements will fluctuate during execution and may exceed these estimates |
+-----------------------------------------------------------------------------+


*********************************** Forces ***********************************
*                                                                            *
*                        Cartesian components (eV/A)                         *
* -------------------------------------------------------------------------- *
*                   x                    y                    z              *
*                                                                            *
* Ag        1     -0.00023(cons'd)      0.00040(cons'd)      0.03095(cons'd) *
* Ag        2     -0.00060(cons'd)      0.00103(cons'd)     -0.01494(cons'd) *
* Ag        3     -0.00005              0.00008             -0.01818         *
* Ag        4      0.00005             -0.00009              0.02067         *
* Ag        5      0.00060             -0.00103             -0.00233         *
* Ag        6      0.00023             -0.00040             -0.01616         *
* Ag        7     -0.00023(cons'd)      0.00041(cons'd)      0.03095(cons'd) *
* Ag        8     -0.00060(cons'd)      0.00104(cons'd)     -0.01494(cons'd) *
* Ag        9     -0.00005              0.00008             -0.01818         *
* Ag       10      0.00005             -0.00008              0.02067         *
* Ag       11      0.00060             -0.00103             -0.00233         *
* Ag       12      0.00023             -0.00040             -0.01616         *
* Ag       13     -0.00023(cons'd)      0.00041(cons'd)      0.03095(cons'd) *
* Ag       14     -0.00060(cons'd)      0.00104(cons'd)     -0.01494(cons'd) *
* Ag       15     -0.00005              0.00008             -0.01818         *
* Ag       16      0.00005             -0.00008              0.02067         *
* Ag       17      0.00060             -0.00103             -0.00233         *
* Ag       18      0.00023             -0.00040             -0.01616         *
* Ag       19     -0.00023(cons'd)      0.00040(cons'd)      0.03095(cons'd) *
* Ag       20     -0.00060(cons'd)      0.00103(cons'd)     -0.01494(cons'd) *
* Ag       21     -0.00005              0.00008             -0.01818         *
* Ag       22      0.00005             -0.00008              0.02067         *
* Ag       23      0.00060             -0.00103             -0.00233         *
* Ag       24      0.00023             -0.00040             -0.01616         *
* Ag       25     -0.00023(cons'd)      0.00040(cons'd)      0.03095(cons'd) *
* Ag       26     -0.00060(cons'd)      0.00103(cons'd)     -0.01494(cons'd) *
* Ag       27     -0.00005              0.00008             -0.01818         *
* Ag       28      0.00005             -0.00008              0.02067         *
* Ag       29      0.00060             -0.00103             -0.00233         *
* Ag       30      0.00023             -0.00040             -0.01616         *
* Ag       31     -0.00023(cons'd)      0.00041(cons'd)      0.03095(cons'd) *
* Ag       32     -0.00060(cons'd)      0.00104(cons'd)     -0.01494(cons'd) *
* Ag       33     -0.00005              0.00008             -0.01818         *
* Ag       34      0.00005             -0.00008              0.02067         *
* Ag       35      0.00060             -0.00103             -0.00233         *
* Ag       36      0.00023             -0.00040             -0.01616         *
* Ag       37     -0.00023(cons'd)      0.00041(cons'd)      0.03095(cons'd) *
* Ag       38     -0.00060(cons'd)      0.00103(cons'd)     -0.01494(cons'd) *
* Ag       39     -0.00005              0.00008             -0.01818         *
* Ag       40      0.00005             -0.00008              0.02067         *
* Ag       41      0.00060             -0.00103             -0.00233         *
* Ag       42      0.00023             -0.00040             -0.01616         *
* Ag       43     -0.00023(cons'd)      0.00040(cons'd)      0.03095(cons'd) *
* Ag       44     -0.00060(cons'd)      0.00103(cons'd)     -0.01494(cons'd) *
* Ag       45     -0.00004              0.00008             -0.01818         *
* Ag       46      0.00005             -0.00008              0.02067         *
* Ag       47      0.00060             -0.00103             -0.00233         *
* Ag       48      0.00023             -0.00040             -0.01616         *
* Ag       49     -0.00023(cons'd)      0.00040(cons'd)      0.03095(cons'd) *
* Ag       50     -0.00060(cons'd)      0.00103(cons'd)     -0.01494(cons'd) *
* Ag       51     -0.00004              0.00008             -0.01818         *
* Ag       52      0.00005             -0.00008              0.02067         *
* Ag       53      0.00060             -0.00104             -0.00233         *
* Ag       54      0.00023             -0.00040             -0.01616         *
*                                                                            *
******************************************************************************
BFGS: finished iteration     0 with enthalpy= -5.55329938E+004 eV

+-----------+-----------------+-----------------+------------+-----+ <-- BFGS
| Parameter |      value      |    tolerance    |    units   | OK? | <-- BFGS
+-----------+-----------------+-----------------+------------+-----+ <-- BFGS
|  dE/ion   |   0.000000E+000 |   1.000000E-005 |         eV | Yes | <-- BFGS
|  |F|max   |   2.067222E-002 |   3.000000E-002 |       eV/A | Yes | <-- BFGS
|  |dR|max  |   0.000000E+000 |   1.000000E-003 |          A | Yes | <-- BFGS
+-----------+-----------------+-----------------+------------+-----+ <-- BFGS

BFGS: Geometry optimization completed successfully.

================================================================================
BFGS: Final Configuration:
================================================================================

                           -------------------------------
                                     Cell Contents
                           -------------------------------

            xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
            x  Element    Atom        Fractional coordinates of atoms  x
            x            Number           u          v          w      x
            x----------------------------------------------------------x
            x  Ag           1         0.222222   0.111111   0.000000   x
            x  Ag           2         0.111111   0.222222   0.091612   x
            x  Ag           3        -0.000000  -0.000000   0.183223   x
            x  Ag           4         0.222222   0.111111   0.274835   x
            x  Ag           5         0.111111   0.222222   0.366447   x
            x  Ag           6         0.000000   0.000000   0.458058   x
            x  Ag           7         0.555556   0.111111   0.000000   x
            x  Ag           8         0.444444   0.222222   0.091612   x
            x  Ag           9         0.333333  -0.000000   0.183223   x
            x  Ag          10         0.555556   0.111111   0.274835   x
            x  Ag          11         0.444444   0.222222   0.366447   x
            x  Ag          12         0.333333   0.000000   0.458058   x
            x  Ag          13         0.888889   0.111111   0.000000   x
            x  Ag          14         0.777778   0.222222   0.091612   x
            x  Ag          15         0.666667   0.000000   0.183223   x
            x  Ag          16         0.888889   0.111111   0.274835   x
            x  Ag          17         0.777778   0.222222   0.366447   x
            x  Ag          18         0.666667   0.000000   0.458058   x
            x  Ag          19         0.222222   0.444444   0.000000   x
            x  Ag          20         0.111111   0.555556   0.091612   x
            x  Ag          21        -0.000000   0.333333   0.183223   x
            x  Ag          22         0.222222   0.444444   0.274835   x
            x  Ag          23         0.111111   0.555556   0.366447   x
            x  Ag          24         0.000000   0.333333   0.458058   x
            x  Ag          25         0.555556   0.444444   0.000000   x
            x  Ag          26         0.444444   0.555556   0.091612   x
            x  Ag          27         0.333333   0.333333   0.183223   x
            x  Ag          28         0.555556   0.444444   0.274835   x
            x  Ag          29         0.444444   0.555556   0.366447   x
            x  Ag          30         0.333333   0.333333   0.458058   x
            x  Ag          31         0.888889   0.444444   0.000000   x
            x  Ag          32         0.777778   0.555556   0.091612   x
            x  Ag          33         0.666667   0.333333   0.183223   x
            x  Ag          34         0.888889   0.444444   0.274835   x
            x  Ag          35         0.777778   0.555556   0.366447   x
            x  Ag          36         0.666667   0.333333   0.458058   x
            x  Ag          37         0.222222   0.777778   0.000000   x
            x  Ag          38         0.111111   0.888889   0.091612   x
            x  Ag          39        -0.000000   0.666667   0.183223   x
            x  Ag          40         0.222222   0.777778   0.274835   x
            x  Ag          41         0.111111   0.888889   0.366447   x
            x  Ag          42         0.000000   0.666667   0.458058   x
            x  Ag          43         0.555556   0.777778   0.000000   x
            x  Ag          44         0.444444   0.888889   0.091612   x
            x  Ag          45         0.333333   0.666667   0.183223   x
            x  Ag          46         0.555556   0.777778   0.274835   x
            x  Ag          47         0.444444   0.888889   0.366447   x
            x  Ag          48         0.333333   0.666667   0.458058   x
            x  Ag          49         0.888889   0.777778   0.000000   x
            x  Ag          50         0.777778   0.888889   0.091612   x
            x  Ag          51         0.666667   0.666667   0.183223   x
            x  Ag          52         0.888889   0.777778   0.274835   x
            x  Ag          53         0.777778   0.888889   0.366447   x
            x  Ag          54         0.666667   0.666667   0.458058   x
            xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx


BFGS: Final Enthalpy     = -5.55329938E+004 eV
BFGS: Final <frequency>  unchanged from initial value

*********************************** Forces ***********************************
*                                                                            *
*                        Cartesian components (eV/A)                         *
* -------------------------------------------------------------------------- *
*                   x                    y                    z              *
*                                                                            *
* Ag        1     -0.00023(cons'd)      0.00040(cons'd)      0.03095(cons'd) *
* Ag        2     -0.00060(cons'd)      0.00103(cons'd)     -0.01494(cons'd) *
* Ag        3     -0.00005              0.00008             -0.01818         *
* Ag        4      0.00005             -0.00009              0.02067         *
* Ag        5      0.00060             -0.00103             -0.00233         *
* Ag        6      0.00023             -0.00040             -0.01616         *
* Ag        7     -0.00023(cons'd)      0.00041(cons'd)      0.03095(cons'd) *
* Ag        8     -0.00060(cons'd)      0.00104(cons'd)     -0.01494(cons'd) *
* Ag        9     -0.00005              0.00008             -0.01818         *
* Ag       10      0.00005             -0.00008              0.02067         *
* Ag       11      0.00060             -0.00103             -0.00233         *
* Ag       12      0.00023             -0.00040             -0.01616         *
* Ag       13     -0.00023(cons'd)      0.00041(cons'd)      0.03095(cons'd) *
* Ag       14     -0.00060(cons'd)      0.00104(cons'd)     -0.01494(cons'd) *
* Ag       15     -0.00005              0.00008             -0.01818         *
* Ag       16      0.00005             -0.00008              0.02067         *
* Ag       17      0.00060             -0.00103             -0.00233         *
* Ag       18      0.00023             -0.00040             -0.01616         *
* Ag       19     -0.00023(cons'd)      0.00040(cons'd)      0.03095(cons'd) *
* Ag       20     -0.00060(cons'd)      0.00103(cons'd)     -0.01494(cons'd) *
* Ag       21     -0.00005              0.00008             -0.01818         *
* Ag       22      0.00005             -0.00008              0.02067         *
* Ag       23      0.00060             -0.00103             -0.00233         *
* Ag       24      0.00023             -0.00040             -0.01616         *
* Ag       25     -0.00023(cons'd)      0.00040(cons'd)      0.03095(cons'd) *
* Ag       26     -0.00060(cons'd)      0.00103(cons'd)     -0.01494(cons'd) *
* Ag       27     -0.00005              0.00008             -0.01818         *
* Ag       28      0.00005             -0.00008              0.02067         *
* Ag       29      0.00060             -0.00103             -0.00233         *
* Ag       30      0.00023             -0.00040             -0.01616         *
* Ag       31     -0.00023(cons'd)      0.00041(cons'd)      0.03095(cons'd) *
* Ag       32     -0.00060(cons'd)      0.00104(cons'd)     -0.01494(cons'd) *
* Ag       33     -0.00005              0.00008             -0.01818         *
* Ag       34      0.00005             -0.00008              0.02067         *
* Ag       35      0.00060             -0.00103             -0.00233         *
* Ag       36      0.00023             -0.00040             -0.01616         *
* Ag       37     -0.00023(cons'd)      0.00041(cons'd)      0.03095(cons'd) *
* Ag       38     -0.00060(cons'd)      0.00103(cons'd)     -0.01494(cons'd) *
* Ag       39     -0.00005              0.00008             -0.01818         *
* Ag       40      0.00005             -0.00008              0.02067         *
* Ag       41      0.00060             -0.00103             -0.00233         *
* Ag       42      0.00023             -0.00040             -0.01616         *
* Ag       43     -0.00023(cons'd)      0.00040(cons'd)      0.03095(cons'd) *
* Ag       44     -0.00060(cons'd)      0.00103(cons'd)     -0.01494(cons'd) *
* Ag       45     -0.00004              0.00008             -0.01818         *
* Ag       46      0.00005             -0.00008              0.02067         *
* Ag       47      0.00060             -0.00103             -0.00233         *
* Ag       48      0.00023             -0.00040             -0.01616         *
* Ag       49     -0.00023(cons'd)      0.00040(cons'd)      0.03095(cons'd) *
* Ag       50     -0.00060(cons'd)      0.00103(cons'd)     -0.01494(cons'd) *
* Ag       51     -0.00004              0.00008             -0.01818         *
* Ag       52      0.00005             -0.00008              0.02067         *
* Ag       53      0.00060             -0.00104             -0.00233         *
* Ag       54      0.00023             -0.00040             -0.01616         *
*                                                                            *
******************************************************************************

Writing analysis data to AgGGA__1_1_1___2_.castep_bin

Writing model to AgGGA__1_1_1___2_.check

A BibTeX formatted list of references used in this run has been written to
AgGGA__1_1_1___2_.bib

Initialisation time =     14.35 s
Calculation time    =   5422.77 s
Finalisation time   =      9.03 s
Total time          =   5446.15 s
Peak Memory Use     = 2546412 kB

Overall parallel efficiency rating: Excellent (97%)                             

Data was distributed by:-
k-point (12-way); efficiency rating: Excellent (97%)                           




作者
Author: 卡开发发    时间: 2015-9-12 16:52
一般没啥问题，可能第一步算下来正好能够落在收敛标准内，castep这点还是很让人放心的。
作者
Author: lzc    时间: 2015-9-12 17:14
卡开发发 发表于 2015-9-12 16:52
一般没啥问题，可能第一步算下来正好能够落在收敛标准内，castep这点还是很让人放心的。

谢谢老师，但是这个我是从完美晶体中切出一个面，建立的这样一个带有真空层的slab模型，按理来说这样的模型进行优化应该会有层间距的变化的吧？
作者
Author: 卡开发发    时间: 2015-9-12 17:28
lzc 发表于 2015-9-12 17:14
谢谢老师，但是这个我是从完美晶体中切出一个面，建立的这样一个带有真空层的slab模型，按理来说这样的模 ...

也不一定，有时候也许变化不大，所以正好就这样了，不过这样的情况不多。

非要挑点毛病的话，这种非四方的结构还是KPOINTS取奇数比较好，还有真空层小了点。
作者
Author: jiewei    时间: 2015-9-12 18:00
cut-off energy for mixing                      :  380.0       eV
如果你实在不方向 把 cutoff energy 取高点， 收敛精度取高点再优化就是了。

另外，你重晶体中切的结构，在低收敛标准下 一步结束是正常的。 特别是对于无机分子。 因为好多的赝势就是根据实验晶体来的
作者
Author: lzc    时间: 2015-9-12 18:29
卡开发发 发表于 2015-9-12 17:28
也不一定，有时候也许变化不大，所以正好就这样了，不过这样的情况不多。

非要挑点毛病的话，这种非四 ...

谢谢老师，我再调整一下试试，谢谢谢谢
作者
Author: lzc    时间: 2015-9-12 18:29
jiewei 发表于 2015-9-12 18:00
cut-off energy for mixing                      :  380.0       eV
如果你实在不方向 把 cutoff energy  ...

谢谢谢谢，谢谢老师，我都调整一下试试，谢谢谢谢
作者
Author: 卡开发发    时间: 2015-9-12 19:04
lzc 发表于 2015-9-12 18:29
谢谢老师，我再调整一下试试，谢谢谢谢

截断能和K要测试，默认的fine一般不够。

赝势一般情形分子和固体可以通用，赝势不通用的情形在pseudopotential的选项卡可以找到专门用于分子体系的硬赝势。一般赝势会根据特定泛函的全电子计算来fitting。
作者
Author: lzc    时间: 2015-9-12 20:14
卡开发发 发表于 2015-9-12 19:04
截断能和K要测试，默认的fine一般不够。

赝势一般情形分子和固体可以通用，赝势不通用的情形在pseudop ...

谢谢老师。这样哦，老师，那么k点跟截断能我就是设置不同的k点及截断能进行计算然后进行比较是吗？另外真空层的厚度我是不是也要进行类似的测试，真空层的厚度选取是不是在真空层里面没有电荷密度就可以了？谢谢老师
作者
Author: 卡开发发    时间: 2015-9-12 21:35
lzc 发表于 2015-9-12 20:14
谢谢老师。这样哦，老师，那么k点跟截断能我就是设置不同的k点及截断能进行计算然后进行比较是吗？另外真 ...

对，直至k和截断能都收敛。真空层原则上要测试，不太可能没有密度，金属表面比较容易，上下表面排列一样，你可以当作两个表面是两块金属，拉开到一定距离就没有相互作用了，因此能量、层间距可当作一项标准，另外，work funtion也应该趋于定值，可能小幅震荡，差不多就行了。
作者
Author: lzc    时间: 2015-9-13 08:15
卡开发发 发表于 2015-9-12 21:35
对，直至k和截断能都收敛。真空层原则上要测试，不太可能没有密度，金属表面比较容易，上下表面排列一样 ...

非常感谢老师的详细简单，我先进行测试，有问题再请教您，非常感谢
作者
Author: lzc    时间: 2015-9-15 10:52
jiewei 发表于 2015-9-12 18:00
cut-off energy for mixing                      :  380.0       eV
如果你实在不方向 把 cutoff energy  ...

老师您好，您在这里说指的收敛进度提高是指在electronic\SCF tolerance 选项里面进行修改吗？我原本是选择的fine，请问我应该选到什么级别一般较理想？非常感谢老师。



欢迎光临 计算化学公社 (http://ccc.keinsci.com/) Powered by Discuz! X3.3