应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » MS » castep声子计算未知原因失败
北京石油化工学院2026年研究生招生接收调剂公告
211/3返回列表123下一页
查看: 3006  |  回复: 20
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

怒冲星河

木虫 (正式写手)

[求助] 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


请各位高手指教
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-09-07 11:36:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

怒冲星河

木虫 (正式写手)
求指导啊
回复此楼
2楼2013-09-07 16:20:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

怒冲星河

木虫 (正式写手)
没人吗？？？？？？？？？？？
赞一下 回复此楼
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的回帖

identation

金虫 (著名写手)
★ ★
franch: 金币+2, 谢谢回帖交流，， 2013-09-09 22:07:59
引用回帖:
5楼: Originally posted by 怒冲星河 at 2013-09-08 10:27:58
还是要谢谢你！你说的那个Born有效电荷是什么，我怎么没找到这个选项？...

就是计算绝缘体或者半导体的有效电荷的选项，和电极化有关
赞一下 回复此楼
6楼2013-09-08 10:32:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

怒冲星河

木虫 (正式写手)
引用回帖:
6楼: Originally posted by identation at 2013-09-08 10:32:17
就是计算绝缘体或者半导体的有效电荷的选项，和电极化有关...

没找到
我是按金属计算的，没选那个LO-TO spilitting
赞一下 回复此楼
7楼2013-09-08 10:40:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

怒冲星河

木虫 (正式写手)
求指教啊
回复此楼
8楼2013-09-08 19:18:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mengfc

金虫 (正式写手)

【答案】应助回帖

★ ★
感谢参与，应助指数 +1
franch: 金币+2, 谢谢回帖交流，， 2013-09-09 22:08:10
内存不足。几个原子的体系算声子谱都得二十多个G内存。
赞一下 回复此楼
9楼2013-09-08 23:47:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

怒冲星河

木虫 (正式写手)
引用回帖:
9楼: Originally posted by mengfc at 2013-09-08 23:47:21
内存不足。几个原子的体系算声子谱都得二十多个G内存。

但是有人几个G内存就算出来了
赞一下 回复此楼
10楼2013-09-09 09:17:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 怒冲星河 的主题更新
211/3返回列表123下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 362求调剂一志愿中国石油大学 +4 我要考大 2026-04-06 6/300 2026-04-06 14:11 by 无际的草原
[考研] 0854求调剂 +4 assdll 2026-04-05 4/200 2026-04-06 12:29 by 中飞院空管学院�
[考研] 0703化学调剂325分 +12 15771691647 2026-04-04 13/650 2026-04-06 12:00 by lijunpoly
[考研] 085600，320分求调剂 +16 大馋小子 2026-04-04 17/850 2026-04-06 07:58 by MOF_Catal
[考研] 0703求调剂383分 +9 W55j 2026-04-03 9/450 2026-04-06 06:50 by houyaoxu
[考研] 329求调剂 +17 miaodesi 2026-04-02 20/1000 2026-04-05 18:33 by 蓝云思雨
[考研] 08专硕275调剂 +5 AaAa7420 2026-04-05 5/250 2026-04-05 18:01 by jkddd
[考研] 材料化工306分找合适调剂 +14 沧海轻舟e 2026-04-04 14/700 2026-04-05 09:53 by 朱云虎202
[考研] 283分求调剂 +7 小聂爱学习 2026-04-03 7/350 2026-04-04 21:51 by hemengdong
[考研] 272求调剂 +4 松柏常青5 2026-04-03 4/200 2026-04-04 17:03 by babysonlkd
[考研] 291求调剂 +4 迷蒙木木 2026-04-01 5/250 2026-04-04 15:59 by sihailian3
[考研] 311求调剂一志愿合肥工业大学 +15 秋二十二 2026-03-30 15/750 2026-04-03 10:19 by linyelide
[考研] 复试调剂 +3 bvzz 2026-04-01 3/150 2026-04-03 09:47 by 蓝云思雨
[考研] 285求调剂 +8 AZMK 2026-04-02 11/550 2026-04-02 20:16 by yulian1987
[考研] 279求调剂 +5 傅文秋 2026-04-02 5/250 2026-04-02 18:10 by 笔落锦州
[考研] 一志愿北京科技大学材料学硕328分求调剂 +6 1段时间 2026-03-31 7/350 2026-04-02 13:57 by 3041
[考研] 279求调剂 +7 莫xiao 2026-04-01 7/350 2026-04-01 22:05 by 客尔美德
[考研] 292求调剂 +17 木虫er12138 2026-04-01 17/850 2026-04-01 21:37 by 七度不信任
[考研] 环境工程调剂 +9 hyzzzzzzz. 2026-04-01 9/450 2026-04-01 14:20 by salamander`
[考研] 262求调剂 +9 励志一定发文章 2026-03-31 10/500 2026-04-01 12:22 by sunshine0013
信息提示
关闭
请填处理意见
关闭