24小时热门版块排行榜    

查看: 2794  |  回复: 20
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

怒冲星河

木虫 (正式写手)

[求助] castep声子计算未知原因失败

用castep对Ni3Pt进行声子计算,不知什么原因失败了

文件如下:

Pseudo atomic calculation performed for Ni 3s2 3p6 3d8 4s2

Converged in 39 iterations to a total energy of -4270.8633 eV


Pseudo atomic calculation performed for Pt 5d9 6s1

Converged in 28 iterations to a total energy of -716.8235 eV

Calculation not parallelised.

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

***************************** General Parameters ******************************
  
output verbosity                               : normal  (1)
write checkpoint data to                       : Ni3Pt.check
type of calculation                            : single point energy
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
  
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
output        entropy unit                     : J/mol/K
  
wavefunctions paging                           : none
random number generator seed                   : randomised (212502906)
data distribution                              : optimal for this architecture
optimization strategy                          : balance speed and memory

*********************** Exchange-Correlation Parameters ***********************
  
using functional                               : Perdew Burke Ernzerhof
Divergence correction                          : off
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                   :   860.0000   eV
size of standard grid                          :     1.5000
size of   fine   gmax                          :    22.5361   1/A
largest prime factor in FFT                    :          5
finite basis set correction                    : none

**************************** Electronic Parameters ****************************
  
number of  electrons                           :  64.00   
net charge of system                           :  0.000   
net spin   of system                           :  0.000   
number of  up  spins                           :  32.00   
number of down spins                           :  32.00   
treating system as non-spin-polarized
number of bands                                :         39

********************* 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.5128E-07   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.5128E-08   eV

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

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

                           -------------------------------
                                      Unit Cell
                           -------------------------------
        Real Lattice(A)                      Reciprocal Lattice(1/A)
   3.5930392   0.0000000   0.0000000        1.7487105   0.0000000   0.0000000
   0.0000000   3.5930392   0.0000000        0.0000000   1.7487105   0.0000000
   0.0000000   0.0000000   3.5930392        0.0000000   0.0000000   1.7487105

                       Lattice parameters(A)       Cell Angles
                    a =    3.593039          alpha =   90.000000
                    b =    3.593039          beta  =   90.000000
                    c =    3.593039          gamma =   90.000000

                       Current cell volume =   46.385888       A**3

                           -------------------------------
                                     Cell Contents
                           -------------------------------
                         Total number of ions in cell =    4
                      Total number of species in cell =    2
                        Max number of any one species =    3

            xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
            x  Element    Atom        Fractional coordinates of atoms  x
            x            Number           u          v          w      x
            x----------------------------------------------------------x
            x  Ni           1         0.000000   0.500000   0.500000   x
            x  Ni           2         0.500000   0.000000   0.500000   x
            x  Ni           3         0.500000   0.500000   0.000000   x
            x  Pt           1         0.000000   0.000000   0.000000   x
            xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx


                         No user defined ionic velocities

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

                               Mass of species in AMU
                                    Ni   58.7099991
                                    Pt  195.0899963

                          Electric Quadrupole Moment (Barn)
                                    Ni    0.1620000 Isotope 61
                                    Pt    1.0000000 No Isotope Defined

                          Files used for pseudopotentials:
                                    Ni Ni_00PBE_OP.recpot
                                    Pt Pt_00PBE_OP.recpot

                           -------------------------------
                              k-Points For BZ Sampling
                           -------------------------------
                       MP grid size for SCF calculation is  6  6  6
                         Number of kpoints used =          10

             +++++++++++++++++++++++++++++++++++++++++++++++++++++++
             +  Number       Fractional coordinates        Weight  +
             +-----------------------------------------------------+
             +     1   0.416667   0.416667   0.416667   0.0370370  +
             +     2   0.416667   0.416667   0.250000   0.1111111  +
             +     3   0.416667   0.416667   0.083333   0.1111111  +
             +     4   0.416667   0.250000   0.250000   0.1111111  +
             +     5   0.416667   0.250000   0.083333   0.2222222  +
             +     6   0.416667   0.083333   0.083333   0.1111111  +
             +     7   0.250000   0.250000   0.250000   0.0370370  +
             +     8   0.250000   0.250000   0.083333   0.1111111  +
             +     9   0.250000   0.083333   0.083333   0.1111111  +
             +    10   0.083333   0.083333   0.083333   0.0370370  +
             +++++++++++++++++++++++++++++++++++++++++++++++++++++++

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

                      Number of symmetry operations   =          48
         There are no ionic constraints specified or generated for this cell
                      Maximum deviation from symmetry =  0.00000         ANG
               Point group of crystal =    32: Oh, m-3m, 4/m -3 2/m
           Set iprint > 1 for details on symmetry rotations/translations

                        Centre of mass is NOT constrained

                         Number of cell constraints= 5
                         Cell constraints are:  1 1 1 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 NODE -----------------+
|                                                     Memory          Disk    |
| Model and support data                               39.0 MB         0.0 MB |
| Electronic energy minimisation requirements          21.5 MB        15.9 MB |
|                                               ----------------------------- |
| Approx. total storage required per node              60.5 MB        15.9 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.31061150E+004  0.00000000E+000                         0.78  <-- SCF
      1  -1.32402763E+004  1.14506873E+001   3.35403379E+001       6.39  <-- SCF
      2  -1.35404367E+004  6.40275175E+000   7.50401079E+001      10.64  <-- SCF
      3  -1.35541799E+004  5.90166761E+000   3.43577928E+000      13.70  <-- SCF
      4  -1.35549192E+004  5.88483534E+000   1.84829656E-001      18.59  <-- SCF
      5  -1.35572208E+004  7.21857059E+000   5.75405641E-001      22.23  <-- SCF
      6  -1.35527108E+004  5.83193807E+000  -1.12751156E+000      25.86  <-- SCF
      7  -1.35526964E+004  5.80084663E+000  -3.59850256E-003      30.64  <-- SCF
      8  -1.35526994E+004  5.82002587E+000   7.72961361E-004      34.08  <-- SCF
      9  -1.35526997E+004  5.80707011E+000   5.73070765E-005      38.94  <-- SCF
     10  -1.35526996E+004  5.80906619E+000  -1.79040618E-005      42.38  <-- SCF
     11  -1.35526996E+004  5.80850429E+000  -3.48392411E-006      45.08  <-- SCF
     12  -1.35526996E+004  5.80832691E+000  -1.31979697E-009      47.19  <-- SCF
     13  -1.35526996E+004  5.80804465E+000   2.44199177E-008      49.28  <-- SCF
------------------------------------------------------------------------ <-- SCF

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

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


Writing model to Ni3Pt.check

***************** Symmetrised Forces *****************
*                                                    *
*            Cartesian components (eV/A)             *
* -------------------------------------------------- *
*                   x            y            z      *
*                                                    *
* Ni        1      0.00000      0.00000      0.00000 *
* Ni        2      0.00000      0.00000      0.00000 *
* Ni        3      0.00000      0.00000      0.00000 *
* Pt        1      0.00000      0.00000      0.00000 *
*                                                    *
******************************************************

Writing model to Ni3Pt.check

Writing analysis data to Ni3Pt.castep_bin

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

Initialisation time =      0.59 s
Calculation time    =     49.19 s
Finalisation time   =      0.16 s
Total time          =     49.94 s


请各位高手指教
回复此楼

» 猜你喜欢

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

怒冲星河

木虫 (正式写手)

引用回帖:
12楼: Originally posted by cswithme at 2013-09-10 10:58:55
你把phonon*.CASTEP那个文件的内容帖出来

Ni3Pt_PhonDOS.castep:

Job started on host NGZVDWXG9UJWYWE
at Fri Sep  6 21:25:53 2013

+-------------------------------------------------+
|                                                 |
|      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 6.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                            |
|                                                 |
| Copyright (c) 2000 - 2011                       |
|                                                 |
| 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 i686-windows-msvc2008 on Dec 08 2011

License checkout of MS_castep successful

+-------------------------------------------------+
|                                                 |
|  D D D D    F F F F F  P P P P     T T T T T    |
|  D       D  F          P       P       T        |
|  D       D  F F F F    P P P P         T        |
|  D       D  F          P               T        |
|  D D D D    F          P               T        |
|                                                 |
+-------------------------------------------------+
|                                                 |
| Welcome to Castep Linear Response (DFPT)        |
|                                                 |
| Copyright (c) 2006 - 2011                       |
|                                                 |
| Please cite the following publications in all   |
| work arising from your use of CASTEP LR.        |
|                                                 |
| K. Refson, S. J. Clark and P. R. Tulip          |
| Variational density functional perturbation     |
| theory for dielectrics and lattice dynamics     |
| Phys. Rev. B 73(15), 155114 (2006)              |
+-------------------------------------------------+

Reading continuation data from model file Ni3Pt_PhonDOS.check

Pseudo atomic calculation performed for Ni 3s2 3p6 3d8 4s2

Converged in 39 iterations to a total energy of -4270.8633 eV


Pseudo atomic calculation performed for Pt 5d9 6s1

Converged in 28 iterations to a total energy of -716.8235 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: BACKUP_INTERVAL modified from          0 to        900
Info: BS_WRITE_EIGENVALUES modified from T to F
Info: PHONON_ENERGY_TOL modified from    0.100000E-04 eV/A**2 to     100000.     eV/A**2
Info: PHONON_FORCE_CONSTANT_CUTOFF modified from     0.00000     A to     5.07624     A
Info: PHONON_FINE_METHOD modified from NONE to supercell
Info: PHONON_METHOD modified from LINEARRESPONSE to finitedisplacement
Info: PHONON_CALC_LO_TO_SPLITTING modified from T to F
Info: CALCULATE_BORN_CHARGES modified from T to F
Calculation not parallelised.

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

***************************** General Parameters ******************************
  
output verbosity                               : normal  (1)
continuing from                                : Ni3Pt_PhonDOS.check
write checkpoint data to                       : Ni3Pt_PhonDOS.check
type of calculation                            : phonon calculation
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
  
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
output        entropy unit                     : J/mol/K
  
wavefunctions paging                           : none
data distribution                              : optimal for this architecture
optimization strategy                          : balance speed and memory

*********************** Exchange-Correlation Parameters ***********************
  
using functional                               : Perdew Burke Ernzerhof
Divergence correction                          : off
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                   :   860.0000   eV
finite basis set correction                    : none

**************************** Electronic Parameters ****************************
  
number of  electrons                           :  64.00   
net charge of system                           :  0.000   
net spin   of system                           :  0.000   
number of  up  spins                           :  32.00   
number of down spins                           :  32.00   
treating system as non-spin-polarized
number of bands                                :         39

********************* 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.5128E-07   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.5128E-08   eV

************************** Density Mixing Parameters **************************
  
density-mixing scheme                          : Pulay
max. length of mixing history                  :         20
charge density mixing amplitude                : 0.5000   
cut-off energy for mixing                      :  860.0       eV

****************************** Phonon Parameters ******************************
  
phonon calculation method                      : finite displacement
finite displacement amplitude                  : 0.5292E-02   A
fine phonon calculation method                 : supercell
non-periodic force constant method             : Weighted cumulant sum
LO/TO splitting term                           : not included
acoustic phonon sum rule                       : correct D(q) in recip-space
Born effective charges                         : not calculated
Raman intensities                              : not calculated
phonon k-points                                : use reduced set
phonon DOS                                     : not calculated
backups results every                          :        900   seconds
band convergence tolerance                     : 0.1000E-05   eV

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

                           -------------------------------
                                      Unit Cell
                           -------------------------------
        Real Lattice(A)                      Reciprocal Lattice(1/A)
   3.5930392   0.0000000   0.0000000        1.7487105   0.0000000   0.0000000
   0.0000000   3.5930392   0.0000000        0.0000000   1.7487105   0.0000000
   0.0000000   0.0000000   3.5930392        0.0000000   0.0000000   1.7487105

                       Lattice parameters(A)       Cell Angles
                    a =    3.593039          alpha =   90.000000
                    b =    3.593039          beta  =   90.000000
                    c =    3.593039          gamma =   90.000000

                       Current cell volume =   46.385888       A**3

                           -------------------------------
                                     Cell Contents
                           -------------------------------
                         Total number of ions in cell =    4
                      Total number of species in cell =    2
                        Max number of any one species =    3

            xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
            x  Element    Atom        Fractional coordinates of atoms  x
            x            Number           u          v          w      x
            x----------------------------------------------------------x
            x  Ni           1         0.000000   0.500000   0.500000   x
            x  Ni           2         0.500000   0.000000   0.500000   x
            x  Ni           3         0.500000   0.500000   0.000000   x
            x  Pt           1         0.000000   0.000000   0.000000   x
            xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx


                         No user defined ionic velocities

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

                               Mass of species in AMU
                                    Ni   58.7099991
                                    Pt  195.0899963

                          Electric Quadrupole Moment (Barn)
                                    Ni    0.1620000 Isotope 61
                                    Pt    1.0000000 No Isotope Defined

                          Files used for pseudopotentials:
                                    Ni Ni_00PBE_OP.recpot
                                    Pt Pt_00PBE_OP.recpot

                           -------------------------------
                              k-Points For BZ Sampling
                           -------------------------------
                       MP grid size for SCF calculation is  6  6  6
                         Number of kpoints used =          10

             +++++++++++++++++++++++++++++++++++++++++++++++++++++++
             +  Number       Fractional coordinates        Weight  +
             +-----------------------------------------------------+
             +     1   0.416667   0.416667   0.416667   0.0370370  +
             +     2   0.416667   0.416667   0.250000   0.1111111  +
             +     3   0.416667   0.416667   0.083333   0.1111111  +
             +     4   0.416667   0.250000   0.250000   0.1111111  +
             +     5   0.416667   0.250000   0.083333   0.2222222  +
             +     6   0.416667   0.083333   0.083333   0.1111111  +
             +     7   0.250000   0.250000   0.250000   0.0370370  +
             +     8   0.250000   0.250000   0.083333   0.1111111  +
             +     9   0.250000   0.083333   0.083333   0.1111111  +
             +    10   0.083333   0.083333   0.083333   0.0370370  +
             +++++++++++++++++++++++++++++++++++++++++++++++++++++++

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

                      Number of symmetry operations   =          48
         There are no ionic constraints specified or generated for this cell
                      Maximum deviation from symmetry =  0.00000         ANG
               Point group of crystal =    32: Oh, m-3m, 4/m -3 2/m
           Set iprint > 1 for details on symmetry rotations/translations

                        Centre of mass is NOT constrained

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

                         External pressure/stress (GPa)
                          0.00000   0.00000   0.00000
                                    0.00000   0.00000
                                              0.00000
                                                                        
              The requested supercell has been generated                 
                                                                        

                           -------------------------------
                                      Unit Cell
                           -------------------------------
        Real Lattice(A)                      Reciprocal Lattice(1/A)
   7.1860784   7.1860784   0.0000000        0.4371776   0.4371776  -0.4371776
   7.1860784   0.0000000   7.1860784        0.4371776  -0.4371776   0.4371776
   0.0000000   7.1860784   7.1860784       -0.4371776   0.4371776   0.4371776

                       Lattice parameters(A)       Cell Angles
                    a =   10.162650          alpha =   60.000000
                    b =   10.162650          beta  =   60.000000
                    c =   10.162650          gamma =   60.000000

                       Current cell volume =  742.174206       A**3

                           -------------------------------
                                     Cell Contents
                           -------------------------------
                         Total number of ions in cell =   64
                      Total number of species in cell =    2
                        Max number of any one species =   48

            xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
            x  Element    Atom        Fractional coordinates of atoms  x
            x            Number           u          v          w      x
            x----------------------------------------------------------x
            x  Ni           1         0.000000   0.000000   0.250000   x
            x  Ni           2         0.000000   0.250000   0.000000   x
            x  Ni           3         0.250000   0.000000   0.000000   x
            x  Ni           4         0.500000   0.000000   0.250000   x
            x  Ni           5         0.500000   0.250000   0.000000   x
            x  Ni           6         0.750000   0.000000   0.000000   x
            x  Ni           7         0.000000   0.500000   0.250000   x
            x  Ni           8         0.000000   0.750000   0.000000   x
            x  Ni           9         0.250000   0.500000   0.000000   x
            x  Ni          10         0.000000   0.000000   0.750000   x
            x  Ni          11         0.000000   0.250000   0.500000   x
            x  Ni          12         0.250000   0.000000   0.500000   x
            x  Ni          13         0.250000   0.250000   0.500000   x
            x  Ni          14         0.250000   0.500000   0.250000   x
            x  Ni          15         0.500000   0.250000   0.250000   x
            x  Ni          16         0.500000   0.500000   0.250000   x
            x  Ni          17         0.500000   0.750000   0.000000   x
            x  Ni          18         0.750000   0.500000   0.000000   x
            x  Ni          19         0.500000   0.000000   0.750000   x
            x  Ni          20         0.500000   0.250000   0.500000   x
            x  Ni          21         0.750000   0.000000   0.500000   x
            x  Ni          22         0.750000   0.250000   0.500000   x
            x  Ni          23         0.750000   0.500000   0.250000   x
            x  Ni          24         0.000000   0.250000   0.250000   x
            x  Ni          25         0.000000   0.500000   0.750000   x
            x  Ni          26         0.000000   0.750000   0.500000   x
            x  Ni          27         0.250000   0.500000   0.500000   x
            x  Ni          28         0.250000   0.750000   0.500000   x
            x  Ni          29         0.250000   0.000000   0.250000   x
            x  Ni          30         0.500000   0.750000   0.250000   x
            x  Ni          31         0.250000   0.250000   0.000000   x
            x  Ni          32         0.250000   0.500000   0.750000   x
            x  Ni          33         0.500000   0.250000   0.750000   x
            x  Ni          34         0.500000   0.500000   0.750000   x
            x  Ni          35         0.500000   0.750000   0.500000   x
            x  Ni          36         0.750000   0.500000   0.500000   x
            x  Ni          37         0.750000   0.750000   0.500000   x
            x  Ni          38         0.750000   0.000000   0.250000   x
            x  Ni          39         0.000000   0.750000   0.250000   x
            x  Ni          40         0.750000   0.250000   0.000000   x
            x  Ni          41         0.750000   0.500000   0.750000   x
            x  Ni          42         0.000000   0.250000   0.750000   x
            x  Ni          43         0.250000   0.750000   0.000000   x
            x  Ni          44         0.250000   0.000000   0.750000   x
            x  Ni          45         0.500000   0.750000   0.750000   x
            x  Ni          46         0.750000   0.750000   0.000000   x
            x  Ni          47         0.750000   0.000000   0.750000   x
            x  Ni          48         0.000000   0.750000   0.750000   x
            x  Pt           1         0.000000   0.000000   0.000000   x
            x  Pt           2         0.500000   0.000000   0.000000   x
            x  Pt           3         0.000000   0.500000   0.000000   x
            x  Pt           4         0.000000   0.000000   0.500000   x
            x  Pt           5         0.250000   0.250000   0.250000   x
            x  Pt           6         0.500000   0.500000   0.000000   x
            x  Pt           7         0.500000   0.000000   0.500000   x
            x  Pt           8         0.750000   0.250000   0.250000   x
            x  Pt           9         0.000000   0.500000   0.500000   x
            x  Pt          10         0.250000   0.750000   0.250000   x
            x  Pt          11         0.250000   0.250000   0.750000   x
            x  Pt          12         0.500000   0.500000   0.500000   x
            x  Pt          13         0.750000   0.750000   0.250000   x
            x  Pt          14         0.750000   0.250000   0.750000   x
            x  Pt          15         0.250000   0.750000   0.750000   x
            x  Pt          16         0.750000   0.750000   0.750000   x
            xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx


                         No user defined ionic velocities

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

                               Mass of species in AMU
                                    Ni   58.7099991
                                    Pt  195.0899963

                          Electric Quadrupole Moment (Barn)
                                    Ni    0.1620000 Isotope 61
                                    Pt    1.0000000 No Isotope Defined

                          Files used for pseudopotentials:
                                    Ni Ni_00PBE_OP.recpot
                                    Pt Pt_00PBE_OP.recpot

                           -------------------------------
                              k-Points For BZ Sampling
                           -------------------------------
                       MP grid size for SCF calculation is  2  2  2
                         Number of kpoints used =           2

             +++++++++++++++++++++++++++++++++++++++++++++++++++++++
             +  Number       Fractional coordinates        Weight  +
             +-----------------------------------------------------+
             +     1  -0.250000  -0.250000  -0.250000   0.2500000  +
             +     2   0.250000  -0.250000   0.250000   0.7500000  +
             +++++++++++++++++++++++++++++++++++++++++++++++++++++++

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

                      Number of symmetry operations   =         768
         There are no ionic constraints specified or generated for this cell
                      Maximum deviation from symmetry =  0.00000         ANG
               Point group of crystal =    32: Oh, m-3m, 4/m -3 2/m
           Set iprint > 1 for details on symmetry rotations/translations

                        Centre of mass is NOT constrained

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

                         External pressure/stress (GPa)
                          0.00000   0.00000   0.00000
                                    0.00000   0.00000
                                              0.00000
************************************************************************
   The memory required for this calculation exceeds machine capacity.

           Unable to continue. Run will be ABORTED.

    Suggestions:
      Check the run setup is reasonable; in particular cell dimensions,
      number of atoms, plane wave cutoff all affect memory consumption.

      Investigate the parameters "opt_strategy"/"opt_strategy_bias" and
      "page_wvfns" which may be used to reduce memory requirements.

      If the calculation is being run on a distributed memory computer
      consider increasing the number of processors.

      It may be necessary to run this calculation on another computer
      with more RAM or upgrade the physical RAM on the current system.
************************************************************************

Information for developers:
ALLOCATE request failed for array "temp_coeffs" in routine "wave_beta_phi_wv_ks"
Current trace stack:
wave_beta_phi_wv_ks
wave_calc_Soverlap_wv_ks
wave_Sorthonormalise_wv_ks
wave_Sorthonormalise_wv
wave_initialise_wv
phonon_force_constant_matrix
phonon_calculate
castep
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 0
13楼2013-09-10 11:30:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 21 个回答

怒冲星河

木虫 (正式写手)

没人吗???????????
3楼2013-09-07 19:44:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

identation

金虫 (著名写手)

★ ★ ★ ★ ★
franch: 金币+5, 谢谢回帖交流,, 2013-09-09 22:07:49
这个不太清楚,没用CASTEP算过声子谱
出错的可能原因包括:1.体系太大,内存不够;2.你计算的是金属或者合金体系,没有gap,但是你勾选了Born有效电荷选项,这也会出错;3.对称性的问题
你到log文件中看看能不能搜到error的具体信息,再分析。。。
另外,你试试把截断能改小了试试看

我用VASP计算过声子谱,没用过CASTEP,就不应助了
4楼2013-09-08 09:26:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

怒冲星河

木虫 (正式写手)

引用回帖:
4楼: Originally posted by identation at 2013-09-08 09:26:06
这个不太清楚,没用CASTEP算过声子谱
出错的可能原因包括:1.体系太大,内存不够;2.你计算的是金属或者合金体系,没有gap,但是你勾选了Born有效电荷选项,这也会出错;3.对称性的问题
你到log文件中看看能不能搜 ...

还是要谢谢你!你说的那个Born有效电荷是什么,我怎么没找到这个选项?
5楼2013-09-08 10:27:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见