24小时热门版块排行榜    

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

landian666

新虫 (小有名气)

[求助] castep计算Fe2O3声子谱失败,原因求教~~

Fe2O3用castep算声子性质失败,对这个物质优化成功了,算能带、光学性质成功了,并且算Al2O3的声子谱时是成功的,这说明不是软件本身的问题,应该是对这个物质的设置问题吧,请指教了~~~多谢~~~

文件错误提示如下,
文件haematite_Efield.param内容
task : Phonon+Efield
continuation : default
spin_polarized : false
elec_energy_tol :   1.000000000000000e-006
cut_off_energy :      900.000000000000000
efield_ignore_molec_modes : crystal
efield_calc_ion_permittivity : true
phonon_max_cycles : 100
phonon_sum_rule : true
efield_energy_tol :   1.000000000000000e-005
excited_state_scissors :   0.000000000000000e+000
calculate_born_charges : true
born_charge_sum_rule : true
bs_max_iter : 200
efield_max_cycles : 500
calculate_raman : true
backup_interval : 900
pdos_calculate_weights : false
num_dump_cycles : 0
bs_write_eigenvalues : false


文件haematite_PhonDisp.param内容
task : Phonon
continuation : default
spin_polarized : false
elec_energy_tol :   1.000000000000000e-006
cut_off_energy :      900.000000000000000
calculate_stress : false
calculate_ELF : false
popn_calculate : false
calculate_densdiff : false
phonon_calc_lo_to_splitting : true
phonon_sum_rule : true
born_charge_sum_rule : true
backup_interval : 900
phonon_energy_tol :   1.000000000000000e-005
phonon_fine_method : interpolate
phonon_force_constant_cutoff :       11.890774075710638  ang
phonon_max_cycles : 100
bs_max_iter : 250
bs_max_cg_steps : 25
bs_eigenvalue_tol : 1.0e-9
pdos_calculate_weights : false
num_dump_cycles : 0
bs_write_eigenvalues : false


文件haematite_PhonDOS.castep后半段内容
      Symmetry and Constraints
                           -------------------------------

                      Number of symmetry operations   =          12
         There are no ionic constraints specified or generated for this cell
                      Maximum deviation from symmetry = 0.329003E-14     ANG
               Point group of crystal =    20: D3d, -3m, -3 2/m
           Set iprint > 1 for details on symmetry rotations/translations

                        Centre of mass is NOT constrained

                         Number of cell constraints= 4
                         Cell constraints are:  1 1 1 4 4 4

                         External pressure/stress (GPa)
                          0.00000   0.00000   0.00000
                                    0.00000   0.00000
                                              0.00000
Phonon perturbation    1 (of    5): O  ion   1 X at q-point (  0.000  0.000  0.000) :       1.89 s
  Parallel strategy changed for PHONON calculation
      32 K-points are distributed over    2 groups, each containing    1 nodes.
  +---------------- MEMORY AND SCRATCH DISK ESTIMATES PER NODE -----------------+
  |                                                     Memory          Disk    |
  | Model and support data                               54.8 MB         0.0 MB |
  | DFPT phonon calculation requirements                 39.3 MB       629.0 MB |
  |                                               ----------------------------- |
  | Approx. total storage required per node              94.1 MB       629.0 MB |
  |                                                                             |
  | Requirements will fluctuate during execution and may exceed these estimates |
  +-----------------------------------------------------------------------------+
  Secondd_get_u0: Warning: Band structure did not converge for    1 out of   10 kpoints
                   Consider increasing BS_MAX_ITER or BS_EIGENVALUE_TOL.

Secondd_minimise_phonon_E2: Can not converge to requested precision in 100 iterations
        Best convergence achievable is 0.541902E-02eV/A**2   
        Consider increasing PHONON_MAX_CYCLES or convergence tolerance PHONON_ENERGY_TOL


************************************************************************
Writing checkpoint file
************************************************************************

Writing model to haematite_PhonDOS.check
Error in PHONON_LINEAR_RESPONSE.  DFPT convergence failed - aborting
Current trace stack:
phonon_linear_response
phonon_dynamical_matrix
phonon_calculate_std_kpoints
phonon_calculate
castep
Error in PHONON_LINEAR_RESPONSE.  DFPT convergence failed - aborting
Current trace stack:
phonon_linear_response
phonon_dynamical_matrix
phonon_calculate_std_kpoints
phonon_calculate
castep
'-max_procs' 不是内部或外部命令,也不是可运行的程序
或批处理文件。
OMP: Warning #2: Cannot open message catalog "2052\libiomp5ui.dll":
OMP: System error #126: 找不到指定的模块。
OMP: Warning #2: Cannot open message catalog "2052\libiomp5ui.dll":
OMP: System error #126: 找不到指定的模块。
OMP: Info #3: Default messages will be used.
OMP: Info #3: Default messages will be used.


************************************************************************
Writing checkpoint file
************************************************************************
MPI Application rank 0
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

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

landian666

新虫 (小有名气)

引用回帖:
4楼: Originally posted by Toapollo at 2011-10-01 19:11:34:
提示的主要问题是不收敛!
你可以从以下几个方面试一下:
固定坐标参数;
增大循环迭代的次数,如180或以上,默认的次数对你的体系太小了;
适当修改一下smear参数。
此问题应该就可以解决了。

我再试试,我看下面问题提示好像是平行计算问题
5楼2011-10-01 23:25:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答

landian666

新虫 (小有名气)

顶一下,急需求解,希望高人指点呀~~多谢,多谢~~

我这两天又试了一下,用一个核算,还是出现类似上面的错误提示~~
2楼2011-10-01 11:54:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)

【答案】应助回帖


uuv2010(金币+1): 多谢提示 2011-10-01 14:47:14
landian666(金币+2): 多谢~~ 2011-10-01 23:24:32
使用CASTEP的线性响应?CASTEP的线性响应只能计算模守赝势,而CASTEP提供的模守赝势烂到不行……不收敛是正常的……如果提高K点,提高截断能还不收敛,就换Abinit或者PWscf吧~
或者,LZ可以考虑使用finited displacement方法计算,那个可以使用超软赝势,只不过对内存的要求比较高~
思想重于技巧,内涵重于表象
3楼2011-10-01 13:07:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Toapollo

铁杆木虫 (正式写手)

西部牛仔

【答案】应助回帖

★ ★
youzhizhe(金币+2): 谢谢交流。 2011-10-01 20:55:29
landian666(金币+3): 多谢~~ 2011-10-01 23:24:58
提示的主要问题是不收敛!
你可以从以下几个方面试一下:
固定坐标参数;
增大循环迭代的次数,如180或以上,默认的次数对你的体系太小了;
适当修改一下smear参数。
此问题应该就可以解决了。
偷得一日闲便是福,枉作千年计岂非愚!
4楼2011-10-01 19:11:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见