求助，HSE06计算吸收光谱出错

慕容轩逸

Job started on host dell-PC
at Thu Dec 10 20:20:04 2020

+-------------------------------------------------+
|                                                 |
|      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           |
|                                                 |
+-------------------------------------------------+
|                                                 |
| Materials Studio CASTEP version 2017            |
| 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, I. J. Bush, C. J. Armstrong,          |
| E. J. Higgins, E. L. Brown, M. S. McFly         |
|                                                 |
| Copyright (c) 2000 - 2016                       |
|                                                 |
| 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 May 05 2016
Intel(R) Math Kernel Library Version 11.1.2
Fundamental constants values: CODATA 2010

License checkout of MS_castep successful

Reading continuation data from model file shoudong_Optics.check

Pseudo atomic calculation performed for O 2s2 2p4

Converged in 24 iterations to a total energy of -421.4504 eV


Pseudo atomic calculation performed for V 3d3 4s2

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


Pseudo atomic calculation performed for Ag 4s2 4p6 4d10 5s1

Converged in 44 iterations to a total energy of -3944.4304 eV

Initialising basis set for model.
Reinitialising basis set for current geometry.
Continuing from previous run with ground state wavefunction.
Continuing from previous run with ground state density.
Info: CALCULATE_HIRSHFELD modified from T to F
Info: POPN_CALCULATE modified from T to F
Info: SPECTRAL_TASK modified from BANDSTRUCTURE to Optics
Info: SPECTRAL_NBANDS modified from        238 to        291
Info: SPECTRAL_EIGENVALUE_TOL modified from    0.1000000000000000E-05 eV to    0.1000000000000000E-04 eV
Calculation parallelised over 32 processes.
Data is distributed by G-vector(32-way)

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

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

output verbosity                               : normal  (1)
continuing from                                : shoudong_Optics.check
write checkpoint data to                       : shoudong_Optics.check
type of calculation                            : Electronic Spectroscopy
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
write OTFG pseudopotential files               : 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                   :   83709515
data distribution                              : optimal for this architecture
optimization strategy                          : maximize speed(+++)

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

using functional                               : hybrid HSE06
k-screen averaging scheme                      : none - use natural jellium density
k-screening update                             : off
relativistic treatment                         : Koelling-Harmon
DFT+D: Semi-empirical dispersion correction    : off

************************* Pseudopotential Parameters **************************

pseudopotential representation                 : reciprocal space
<beta|phi> representation                      : reciprocal space
spin-orbit coupling                            : off

**************************** Basis Set Parameters *****************************

plane wave basis set cut-off                   :   765.0000   eV
size of standard grid                          :     1.5000
size of   fine   gmax                          :    21.2550   1/A
finite basis set correction                    : none

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

number of  electrons                           :  344.0   
net charge of system                           :  0.000   
treating system as non-spin-polarized
number of bands                                :        172

********************* 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
max force / atom convergence tol.              : ignored
convergence tolerance window                   :          3   cycles
max. number of SCF cycles                      :        250
periodic dipole correction                     : NONE

********************* Electronic Spectroscopy Parameters **********************

electronic spectroscopy with theory level      : DFT      
spectroscopy calculation                       : optics
max. number of iterations                      :         60
max. steps per iteration                       :          5
number of bands / k-point                      :        291
band convergence tolerance                     : 0.1000E-04   eV
write orbitals file                            : off
parallelepiped integration                     : on
x-size of parallelepiped                       :          8
y-size of parallelepiped                       :          8
z-size of parallelepiped                       :          8
k-screen averaging scheme                      : none - use natural jellium density
k-screening update                             : off

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


                           -------------------------------
                                      Unit Cell
                           -------------------------------
        Real Lattice(A)                      Reciprocal Lattice(1/A)
   9.4209605   0.0000000  -4.9986036        0.6669368   0.0000000   0.0000000
   0.0000000   5.0242871   0.0000000        0.0000000   1.2505626   0.0000000
   0.0000000   0.0000000  10.4078449        0.3203115   0.0000000   0.6036971

                       Lattice parameters(A)       Cell Angles
                    a =   10.664921          alpha =   90.000000
                    b =    5.024287          beta  =  117.949672
                    c =   10.407845          gamma =   90.000000

                       Current cell volume =  492.640876       A**3

                           -------------------------------
                                     Cell Contents
                           -------------------------------
                         Total number of ions in cell =   32
                      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  O            1         0.408095  -0.158924   0.098568   x
            x  O            2         0.381968   0.247824   0.284806   x
            x  O            3         0.908095   0.341076   0.098568   x
            x  O            4         0.881968   0.747824   0.284806   x
            x  O            5        -0.408095  -0.158924   0.401432   x
            x  O            6        -0.381968   0.247824   0.215194   x
            x  O            7         0.091905   0.341076   0.401432   x
            x  O            8         0.118032   0.747824   0.215194   x
            x  O            9        -0.408095   0.158924  -0.098568   x
            x  O           10        -0.381968  -0.247824  -0.284806   x
            x  O           11         0.091905   0.658924  -0.098568   x
            x  O           12         0.118032   0.252176  -0.284806   x
            x  O           13         0.408095   0.158924   0.598568   x
            x  O           14         0.381968  -0.247824   0.784806   x
            x  O           15         0.908095   0.658924   0.598568   x
            x  O           16         0.881968   0.252176   0.784806   x
            x  V            1         0.500000   0.047967   0.250000   x
            x  V            2         1.000000   0.547967   0.250000   x
            x  V            3        -0.500000  -0.047967  -0.250000   x
            x  V            4         0.000000   0.452033  -0.250000   x
            x  Ag           1         0.262338  -0.044666  -0.593981   x
            x  Ag           2         0.762338   0.455334  -0.593981   x
            x  Ag           3        -0.262338  -0.044666   1.093981   x
            x  Ag           4         0.237662   0.455334   1.093981   x
            x  Ag           5        -0.262338   0.044666   0.593981   x
            x  Ag           6         0.237662   0.544666   0.593981   x
            x  Ag           7         0.262338   0.044666  -0.093981   x
            x  Ag           8         0.762338   0.544666  -0.093981   x
            x  Ag           9         0.000000   0.000000   0.500000   x
            x  Ag          10         0.500000   0.500000   0.500000   x
            x  Ag          11         0.000000   0.000000   0.000000   x
            x  Ag          12         0.500000   0.500000   0.000000   x
            xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx


                         No user defined ionic velocities

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

                               Mass of species in AMU
                                    O   15.9989996
                                    V   50.9410019
                                    Ag  107.8679962

                          Electric Quadrupole Moment (Barn)
                                    O   -0.0255800 Isotope 17
                                    V   -0.0520000 Isotope 51
                                    Ag    1.0000000 No Isotope Defined

                          Files used for pseudopotentials:
                                    O O_00PBE_OP.recpot
                                    V V_00.recpot
                                    Ag Ag_00PBE_OP.recpot

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

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

                      Cell is a supercell containing 2 primitive cells
                      Maximum deviation from symmetry = 0.231101E-14     ANG

                      Number of symmetry operations   =           8
         There are no ionic constraints specified or generated for this cell
                      Point group of crystal =     5: C2h, 2/m, 2/m
                      Space group of crystal =    15: C2/c, -C 2yc

             Set iprint > 1 for details on symmetry rotations/translations

                         Centre of mass is NOT constrained

                         Number of cell constraints= 2
                         Cell constraints are:  1 2 3 0 5 0

                         External pressure/stress (GPa)
                          0.00000   0.00000   0.00000
                                    0.00000   0.00000
                                              0.00000
spectral_calculate_dft: Cannot compute coreless spectra with nonlocal exchange
Current trace stack:
spectral_calculate_dft
castep
spectral_calculate_dft: Cannot compute coreless spectra with nonlocal exchange
Current trace stack:
spectral_calculate_dft
castep
spectral_calculate_dft: Cannot compute coreless spectra with nonlocal exchange
Current trace stack:
spectral_calculate_dft
castep
spectral_calculate_dft: Cannot compute coreless spectra with nonlocal exchange
Current trace stack:
spectral_calculate_dft
castep
spectral_calculate_dft: Cannot compute coreless spectra with nonlocal exchange
Current trace stack:
spectral_calculate_dft
castep
spectral_calculate_dft: Cannot compute coreless spectra with nonlocal exchange
Current trace stack:
spectral_calculate_dft
castep
spectral_calculate_dft: Cannot compute coreless spectra with nonlocal exchange
Current trace stack:
spectral_calculate_dft
castep
spectral_calculate_dft: Cannot compute coreless spectra with nonlocal exchange
Current trace stack:
spectral_calculate_dft
castep
spectral_calculate_dft: Cannot compute coreless spectra with nonlocal exchange
Current trace stack:
spectral_calculate_dft
castep
spectral_calculate_dft: Cannot compute coreless spectra with nonlocal exchange
Current trace stack:
spectral_calculate_dft
castep
spectral_calculate_dft: Cannot compute coreless spectra with nonlocal exchange
Current trace stack:
spectral_calculate_dft
castep
spectral_calculate_dft: Cannot compute coreless spectra with nonlocal exchange
Current trace stack:
spectral_calculate_dft
castep
spectral_calculate_dft: Cannot compute coreless spectra with nonlocal exchange
Current trace stack:
spectral_calculate_dft
castep
spectral_calculate_dft: Cannot compute coreless spectra with nonlocal exchange
Current trace stack:
spectral_calculate_dft
castep
spectral_calculate_dft: Cannot compute coreless spectra with nonlocal exchange
Current trace stack:
spectral_calculate_dft
castep
spectral_calculate_dft: Cannot compute coreless spectra with nonlocal exchange
Current trace stack:
spectral_calculate_dft
castep
spectral_calculate_dft: Cannot compute coreless spectra with nonlocal exchange
Current trace stack:
spectral_calculate_dft
castep
spectral_calculate_dft: Cannot compute coreless spectra with nonlocal exchange
Current trace stack:
spectral_calculate_dft
castep
spectral_calculate_dft: Cannot compute coreless spectra with nonlocal exchange
Current trace stack:
spectral_calculate_dft
castep
spectral_calculate_dft: Cannot compute coreless spectra with nonlocal exchange
Current trace stack:
spectral_calculate_dft
castep
spectral_calculate_dft: Cannot compute coreless spectra with nonlocal exchange
Current trace stack:
spectral_calculate_dft
castep
spectral_calculate_dft: Cannot compute coreless spectra with nonlocal exchange
Current trace stack:
spectral_calculate_dft
castep
spectral_calculate_dft: Cannot compute coreless spectra with nonlocal exchange
Current trace stack:
spectral_calculate_dft
castep
spectral_calculate_dft: Cannot compute coreless spectra with nonlocal exchange
Current trace stack:
spectral_calculate_dft
castep
spectral_calculate_dft: Cannot compute coreless spectra with nonlocal exchange
Current trace stack:
spectral_calculate_dft
castep
spectral_calculate_dft: Cannot compute coreless spectra with nonlocal exchange
Current trace stack:
spectral_calculate_dft
castep
spectral_calculate_dft: Cannot compute coreless spectra with nonlocal exchange
Current trace stack:
spectral_calculate_dft
castep
spectral_calculate_dft: Cannot compute coreless spectra with nonlocal exchange
Current trace stack:
spectral_calculate_dft
castep
spectral_calculate_dft: Cannot compute coreless spectra with nonlocal exchange
Current trace stack:
spectral_calculate_dft
castep
spectral_calculate_dft: Cannot compute coreless spectra with nonlocal exchange
Current trace stack:
spectral_calculate_dft
castep
spectral_calculate_dft: Cannot compute coreless spectra with nonlocal exchange
Current trace stack:
spectral_calculate_dft
castep
spectral_calculate_dft: Cannot compute coreless spectra with nonlocal exchange
Current trace stack:
spectral_calculate_dft
castep
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 11
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 5
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 29
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 17
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 25
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 21
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 31




各位老师，我想问一下计算一直报错是我参数设置的问题，还是函数选取的问题？