MS模拟物质Raman谱，分析时出现错误：.castep file without Raman intensities

Tanmin

各位大家好：

最近在用MS模拟物质Raman谱时，软件显示计算成功，但分析时，出现下图错误：


请各位多多指导

下面是计算参数，经过convergence testing
                task : GeometryOptimization
                comment : CASTEP calculation from Materials Studio
                xc_functional : PW91
                spin_polarized : false
                opt_strategy : Default
                page_wvfns :        0
                cut_off_energy :      700.000000000000000
                grid_scale :        2.000000000000000
                fine_grid_scale :        3.000000000000000
                finite_basis_corr :        2
                finite_basis_npoints :        3
                elec_energy_tol :   2.000000000000000e-006
                max_scf_cycles :      200
                fix_occupancy : true
                nextra_bands : 0
                num_dump_cycles : 0
                geom_energy_tol :   5.000000000000000e-006
                geom_force_tol :        0.010000000000000
                geom_stress_tol :        0.020000000000000
                geom_disp_tol :   5.000000000000000e-004
                geom_max_iter :      100
                geom_method : BFGS
                fixed_npw : false
                geom_modulus_est :      500.000000000000000  GPa
                calculate_ELF : false
                calculate_stress : true
                popn_calculate : false
                calculate_hirshfeld : false
                calculate_densdiff : false
                pdos_calculate_weights : false

                task : Phonon+Efield
                continuation : default
                comment : CASTEP calculation from Materials Studio
                xc_functional : PW91
                spin_polarized : false
                opt_strategy : Default
                page_wvfns :        0
                cut_off_energy :      700.000000000000000
                grid_scale :        2.000000000000000
                fine_grid_scale :        3.000000000000000
                finite_basis_corr :        2
                finite_basis_npoints :        3
                elec_energy_tol :   1.000000000000000e-010
                elec_eigenvalue_tol :   1.000000000000000e-009
                max_scf_cycles :      200
                fix_occupancy : true
                nextra_bands : 0
                num_dump_cycles : 0
                efield_ignore_molec_modes : crystal
                efield_calc_ion_permittivity : true
                phonon_max_cycles : 100
                phonon_sum_rule_method : reciprocal
                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
                raman_range_low :   0.000000000000000e+000 cm-1
                raman_range_high :   1.000000000000000e+004 cm-1
                backup_interval : 900
                pdos_calculate_weights : false
                spectral_write_eigenvalues : false

下面是-efiled.castep 文件中的内容：
                Warning in parameters_validate: current value of
                   ELEC_EIGENVALUE_TOL = 0.100000E-08eV
                   is too large to achieve desired level of convergence of response properties.
                   This may cause convergence failures and/or inaccuracy of results
                   of PHONON+EFIELD calculations - recommend you use a smaller value, e.g.
                   ELEC_EIGENVALUE_TOL = 0.583333E-10eV
               
                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: ELEC_ENERGY_TOL modified from    0.2000000000000000E-05 eV to    0.1000000000000000E-09 eV
                 Info: ELEC_EIGENVALUE_TOL modified from    0.1000000000000000E-05 eV to    0.1000000000000000E-08 eV
                 Info: BS_MAX_ITER modified from         60 to        200
                 Info: PHONON_MAX_CYCLES modified from         50 to        100
                 Info: CALCULATE_RAMAN modified from F to T
                 Info: RAMAN_RANGE_HIGH modified from     1.000000000000000     eV to     1.239841929200420     eV
                 Info: EFIELD_MAX_CYCLES modified from         50 to        500
                 Info: SPECTRAL_WRITE_EIGENVALUES modified from T to F
                 
                 WARNING: Ground state wavefunctions and density data in check file marked as unconverged.
                          because some parameters changed on continuation.
                 
                 Calculation parallelised over 12 processes.
                 Data is distributed by k-point(12-way)
               
                Vibrational Frequencies                          +
                 +                           -----------------------                          +
                 +                                                                            +
                 + Performing frequency calculation at  1 wavevector  (q-pt )                 +
                 + ========================================================================== +
                 +                                                                            +
                 +============================================================================+
                 +                       Raman Intensity Progress Report                      +
                 +----------------------------------------------------------------------------+
                 + RAMAN: Not a continuation; generating Raman intensities from beginning     +
                 +----------------------------------------------------------------------------+
                 + Mode     1: not Raman active; calculation not required                     +
                 + Mode     2: not Raman active; calculation not required                     +
                 + Mode     3: not Raman active; calculation not required                     +
                 + Mode     4: performing Raman intensity calculation                         +
                 + Mode number:   4 completed;                   time taken:  1130.47 seconds +
                 + Checkpointing.....                                                         +
               
                Writing analysis data to Al4C3_mp-632442_Efield.castep_bin
               
                Writing model to Al4C3_mp-632442_Efield.check
                 + Mode     5: performing Raman intensity calculation                         +
                 + Mode number:   5 completed;                   time taken:  1112.86 seconds +
                 + Checkpointing.....                                                         +
……………………………………………………为便于虫友查看，此次省略部分计算中间内容…………………………………………………………

+----------------------------------------------------------------------------+
                 +   E N D    O F    R A M A N    I N T E N S I T Y    C A L C U L A T I O N  +
                 +         Raman intensity calculation completed in   20340.52 seconds        +
                 +============================================================================+
                 + -------------------------------------------------------------------------- +
                 +  q-pt=    1 (  0.000000  0.000000  0.000000)     1.0000000000              +
                 + -------------------------------------------------------------------------- +
                 +  Acoustic sum rule correction <  12.385066 cm-1 applied                    +
                 +     N      Frequency           ir intensity         raman activity         +
                 +                (cm-1)         ((D/A)**2/amu)              (A**4/amu)       +
                 +                                                                            +
                 +     1      -0.038186              0.0000000               0.0000000        +
                 +     2      -0.031196              0.0000000               0.0000000        +
                 +     3      -0.022058              0.0000000               0.0000000        +
                 +     4     119.962794              6.7878338            4278.3122460        +
                 +     5     156.579251              2.9536536            1714.4486269        +
                 +     6     237.309866              0.8729072            3720.6925752        +
                 +     7     239.260406              5.4126247            2858.4290772        +
                 +     8     269.800560              9.1693299           16985.9192068        +
                 +     9     307.592119              4.9199842            5630.1892928        +
                 +    10     324.508650             20.2511256            9473.4043425        +
                 +    11     390.981051             28.3528238           30248.0328004        +
                 +    12     430.409968              7.0370798            3611.2560606        +
                 +    13     458.187024              3.9590171           48080.5970025        +
                 +    14     483.372487             17.9475001           15893.8888328        +
                 +    15     496.447888              7.3642806           12008.9516531        +
                 +    16     522.189379             48.2565305           33424.7536118        +
                 +    17     551.058141             36.1193259            6076.5938836        +
                 +    18     583.399120              4.2745716           12692.2891445        +
                 +    19     601.070436              8.0126592            4388.6698821        +
                 +    20     661.413769             28.3433146           15918.7213986        +
                 +    21     746.522448             21.3010125           35454.5852783        +
                 ==============================================================================
                上图虽然显示了intensity，但经过与能够得出Raman图谱的数据相比，其少了很多内容。
经过自己对上面错误的分析，认为可能是计算精度不够，故对下面的参数做了如下修改。但任务运行后，显示任务运算成功，但分析Raman谱图时，任报相同的错误。
                修改scf: medium to fine
                修改k-point: 5,4,2 to 11,7,7
                修改elec_eigenvalue_tol ：0.58333333333332e-011
                修改max_scf_cycles : 500