|
不好意思,这是从一个已经计算成功的使用vasp软件计算过得CeO2例子,想说明模型并没有错,
Pseudo atomic calculation performed for O 2s2 2p4
Converged in 23 iterations to a total energy of -429.5705 eV
Pseudo atomic calculation performed for Cu 3d10 4s1
Converged in 52 iterations to a total energy of -1472.4645 eV
Pseudo atomic calculation performed for Ce 4f1 5s2 5p6 5d1 6s2
Converged in 25 iterations to a total energy of -1055.3538 eV
Calculation parallelised over 12 processes.
Data is distributed by G-vector(6-way) and k-point(2-way)
************************************ Title ************************************
CASTEP calculation from Materials Studio
***************************** General Parameters ******************************
output verbosity : normal (1)
write checkpoint data to : 111.check
type of calculation : geometry optimization
stress calculation : off
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 (80028079)
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 : 400.0000 eV
size of standard grid : 1.5000
size of fine gmax : 15.3695 1/A
largest prime factor in FFT : 5
finite basis set correction : none
**************************** Electronic Parameters ****************************
number of electrons : 299.0
net charge of system : 0.000
net spin of system : 1.000
number of up spins : 150.0
number of down spins : 149.0
treating system as spin-polarized
number of bands : 196
********************* 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.2000E-05 eV
eigen-energy convergence tolerance : 0.3776E-06 eV
max force / atom convergence tol. : ignored
convergence tolerance window : 3 cycles
max. number of SCF cycles : 2500
number of fixed-spin iterations : 6
smearing scheme : Gaussian
smearing width : 0.2000 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
spin density mixing amplitude : 2.000
cut-off energy for mixing : 400.0 eV
charge density mixing g-vector : 1.500 1/A
spin 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
with spin fully able to relax for all steps
total energy convergence tolerance : 0.2000E-04 eV/atom
max ionic |force| tolerance : 0.5000E-01 eV/A
max ionic |displacement| tolerance : 0.2000E-02 A
max |stress component| tolerance : 0.1000 GPa
convergence tolerance window : 2 steps
backup results every : 5 steps
*******************************************************************************
-------------------------------
Unit Cell
-------------------------------
Real Lattice(A) Reciprocal Lattice(1/A)
6.6270862 -3.8261500 0.0000000 0.9481068 0.0000000 0.0000000
0.0000000 7.6523000 0.0000000 0.4740534 0.8210846 0.0000000
0.0000000 0.0000000 22.8100999 -0.0000000 0.0000000 0.2754563
Lattice parameters(A) Cell Angles
a = 7.652300 alpha = 90.000000
b = 7.652300 beta = 90.000000
c = 22.810100 gamma = 120.000000
Current cell volume = 1156.756085 A**3
-------------------------------
Cell Contents
-------------------------------
Total number of ions in cell = 37
Total number of species in cell = 3
Max number of any one species = 24
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
x Element Atom Fractional coordinates of atoms x
x Number u v w x
x----------------------------------------------------------x
x O 1 0.000419 0.000956 -0.010232 x
x O 2 0.167712 0.335530 0.133622 x
x O 3 0.315672 0.157702 0.279705 x
x O 4 0.190629 0.381343 0.360466 x
x O 5 -0.000686 -0.001272 0.206105 x
x O 6 0.330919 0.165538 0.063395 x
x O 7 0.500546 0.001202 -0.009550 x
x O 8 0.666778 0.333675 0.131393 x
x O 9 0.841929 0.157709 0.279693 x
x O 10 0.665673 0.331397 0.344189 x
x O 11 0.504730 0.009571 0.206619 x
x O 12 0.834504 0.165549 0.063396 x
x O 13 -0.000148 0.500547 -0.009528 x
x O 14 0.167795 0.832251 0.133450 x
x O 15 0.333226 0.666508 0.271275 x
x O 16 0.190462 0.809997 0.360171 x
x O 17 -0.009627 0.495365 0.206869 x
x O 18 0.333279 0.666682 0.063111 x
x O 19 0.500574 0.500547 -0.009526 x
x O 20 0.664357 0.832256 0.133448 x
x O 21 0.841710 0.683491 0.280735 x
x O 22 0.619567 0.810055 0.360157 x
x O 23 0.504901 0.495356 0.206869 x
x O 24 0.834684 0.669477 0.063332 x
x Cu 1 0.334493 0.668936 0.391632 x
x Ce 1 -0.000442 -0.000721 0.313126 x
x Ce 2 0.332024 0.165351 0.167754 x
x Ce 3 0.166907 0.333943 0.026651 x
x Ce 4 0.502493 0.005074 0.308839 x
x Ce 5 0.833227 0.165368 0.167752 x
x Ce 6 0.666906 0.333944 0.026829 x
x Ce 7 -0.004808 0.496315 0.308797 x
x Ce 8 0.333096 0.666274 0.168643 x
x Ce 9 0.166989 0.833846 0.026533 x
x Ce 10 0.501063 0.496320 0.308788 x
x Ce 11 0.833734 0.667560 0.167585 x
x Ce 12 0.666732 0.833849 0.026534 x
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
No user defined ionic velocities
-------------------------------
Details of Species
-------------------------------
Mass of species in AMU
O 15.9989996
Cu 63.5460014
Ce 140.1199951
Electric Quadrupole Moment (Barn)
O -0.0255800 Isotope 17
Cu -0.2200000 Isotope 63
Ce 1.0000000 No Isotope Defined
Files used for pseudopotentials:
O O_00PBE.usp
Cu Cu_00PBE.usp
Ce Ce_00PBE.usp
-------------------------------
k-Points For BZ Sampling
-------------------------------
MP grid size for SCF calculation is 2 2 1
with an offset of 0.000 0.000 0.000
Number of kpoints used = 2
-------------------------------
Symmetry and Constraints
-------------------------------
Maximum deviation from symmetry = 0.00000 ANG
There are no symmetry operations specified or generated for this cell
Number of ionic constraints = 72
Point group of crystal = 1: C1, 1, 1
Space group of crystal = 1: P1, P 1
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 126.2 MB 113.4 MB |
| Electronic energy minimisation requirements 60.1 MB 0.0 MB |
| ----------------------------- |
| Approx. total storage required per process 186.3 MB 113.4 MB |
| |
| Requirements will fluctuate during execution and may exceed these estimates |
*Warning* max. SCF cycles performed but system has not reached the groundstate.
Current total energy, E = -24764.49124797 eV
Current free energy (E-TS) = -24764.50236165 eV
(energies not corrected for finite basis set)
NB est. 0K energy (E-0.5TS) = -24764.49680481 eV
****************************************************************************
Warning: electronic minimisation did not converge when finding ground state.
****************************************************************************
Writing checkpoint file...
Writing analysis data to 111.castep_bin
Writing model to 111.check
Overall parallel efficiency rating: Very good (84%)
Data was distributed by:-
G-vector (6-way); efficiency rating: Very good (85%)
k-point (2-way); efficiency rating: Ideal (100%)
Error check_elec_ground_state : electronic_minimisation of initial cell failed.
Current trace stack:
check_elec_ground_state
castep
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 1
您看这个也是说的空带问题是吧,空带我曾经加到30%了,还显示这样的错误,这个也不能增加太多吧 |
|