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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 求助，对于未relax的结构测试ecutwfc时不收敛
141/2返回列表12下一页
查看: 3449  |  回复: 13
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

frank_zhan

金虫 (正式写手)

[求助] 求助，对于未relax的结构测试ecutwfc时不收敛

请教各位虫友，我要计算一个(13,0)的碳纳米管上修饰10个金原子，我建了一个很大的胞，减小径向的image。在轴向是周期性的。
在没有做relax之前，我先用scf来测试ecutwfc取值的收敛性，但是我选用ecutwfc=30D0, ecutrho=300.D0,计算了100不步后，发现不收敛啊。
代码下：
&control
    calculation = 'scf'
    restart_mode= 'from_scratch'
    prefix= 'SWNT-Au'
    tstress = .true.
    tprnfor = .true.
    pseudo_dir = '/home/zhan0273/espresso-4.1.2/pseudo',
    outdir= '/scratch/zhan0273/PWSCF/CNT/tmp'
/
&system
    ibrav   = 8,
celldm(1) = 39.6842,
celldm(2) = 1.1904,
celldm(3) = 0.8778,
    nat    = 217,
    ntyp   = 2,
  ecutwfc   = 30.D0,
ecutrho=300.D0,
  occupations= 'smearing', smearing= 'fd', degauss=0.02D0
/
&electrons
    mixing_mode = 'plain'
    mixing_beta = 0.7
    conv_thr =  1.0d-7
/
ATOMIC_SPECIES
C  12.0107  C.pbe-rrkjus.UPF
Au 196.966  Au.pbe-nd-rrkjus.UPF
ATOMIC_POSITIONS {bohr}
C        28.6714        26.2825        0.6743
C        29.2262        24.0316        4.6994
C        29.2262        24.0316        7.3828
C        26.0568        30.0705        0.6743
C        28.6714        26.2825        3.3577
C        27.594        28.3352        4.6994
C        27.594        28.3352        7.3828
C        21.9813        32.2094        0.6743
C        26.0568        30.0705        3.3577
C        24.1489        31.3874        4.6994
C        24.1489        31.3874        7.3828
C        17.3786        32.2094        0.6743
C        21.9813        32.2094        3.3577
C        19.6799        32.4889        4.6994
C        19.6799        32.4889        7.3828
C        13.3031        30.0705        0.6743
C        17.3786        32.2094        3.3577
C        15.211        31.3874        4.6994
C        15.211        31.3874        7.3828
C        10.6885        26.2825        0.6743
C        13.3031        30.0705        3.3577
C        11.7658        28.3352        4.6994
C        11.7658        28.3352        7.3828
C        10.1337        21.7134        0.6743
C        10.6885        26.2825        3.3577
C        10.1337        24.0316        4.6994
C        10.1337        24.0316        7.3828
C        11.7658        17.4098        0.6743
C        10.1337        21.7134        3.3577
C        10.6885        19.4625        4.6994
C        10.6885        19.4625        7.3828
C        15.211        14.3576        0.6743
C        11.7658        17.4098        3.3577
C        13.3031        15.6746        4.6994
C        13.3031        15.6746        7.3828
C        19.6799        13.2561        0.6743
C        15.211        14.3576        3.3577
C        17.3786        13.5356        4.6994
C        17.3786        13.5356        7.3828
C        24.1489        14.3576        0.6743
C        19.6799        13.2561        3.3577
C        21.9813        13.5356        4.6994
C        21.9813        13.5356        7.3828
C        27.594        17.4098        0.6743
C        24.1489        14.3576        3.3577
C        26.0568        15.6746        4.6994
C        26.0568        15.6746        7.3828
C        29.2262        21.7134        0.6743
C        27.594        17.4098        3.3577
C        28.6714        19.4625        4.6994
C        28.6714        19.4625        7.3828
C        29.2262        21.7134        3.3577
C        28.6714        26.2825        8.7245
C        29.2262        24.0316        12.7497
C        29.2262        24.0316        15.4331
C        26.0568        30.0705        8.7245
C        28.6714        26.2825        11.408
C        27.594        28.3352        12.7497
C        27.594        28.3352        15.4331
C        21.9813        32.2094        8.7245
C        26.0568        30.0705        11.408
C        24.1489        31.3874        12.7497
C        24.1489        31.3874        15.4331
C        17.3786        32.2094        8.7245
C        21.9813        32.2094        11.408
C        19.6799        32.4889        12.7497
C        19.6799        32.4889        15.4331
C        13.3031        30.0705        8.7245
C        17.3786        32.2094        11.408
C        15.211        31.3874        12.7497
C        15.211        31.3874        15.4331
C        10.6885        26.2825        8.7245
C        13.3031        30.0705        11.408
C        11.7658        28.3352        12.7497
C        11.7658        28.3352        15.4331
C        10.1337        21.7134        8.7245
C        10.6885        26.2825        11.408
C        10.1337        24.0316        12.7497
C        10.1337        24.0316        15.4331
C        11.7658        17.4098        8.7245
C        10.1337        21.7134        11.408
C        10.6885        19.4625        12.7497
C        10.6885        19.4625        15.4331
C        15.211        14.3576        8.7245
C        11.7658        17.4098        11.408
C        13.3031        15.6746        12.7497
C        13.3031        15.6746        15.4331
C        19.6799        13.2561        8.7245
C        15.211        14.3576        11.408
C        17.3786        13.5356        12.7497
C        17.3786        13.5356        15.4331
C        24.1489        14.3576        8.7245
C        19.6799        13.2561        11.408
C        21.9813        13.5356        12.7497
C        21.9813        13.5356        15.4331
C        27.594        17.4098        8.7245
C        24.1489        14.3576        11.408
C        26.0568        15.6746        12.7497
C        26.0568        15.6746        15.4331
C        29.2262        21.7134        8.7245
C        27.594        17.4098        11.408
C        28.6714        19.4625        12.7497
C        28.6714        19.4625        15.4331
C        29.2262        21.7134        11.408
C        28.6714        26.2825        16.7748
C        29.2262        24.0316        20.7999
C        29.2262        24.0316        23.4833
C        26.0568        30.0705        16.7748
C        28.6714        26.2825        19.4582
C        27.594        28.3352        20.7999
C        27.594        28.3352        23.4833
C        26.0568        30.0705        19.4582
C        24.1489        31.3874        20.7999
C        24.1489        31.3874        23.4833
C        17.3786        32.2094        16.7748
C        21.9813        32.2094        19.4582
C        19.6799        32.4889        20.7999
C        19.6799        32.4889        23.4833
C        13.3031        30.0705        16.7748
C        17.3786        32.2094        19.4582
C        15.211        31.3874        20.7999
C        15.211        31.3874        23.4833
C        10.6885        26.2825        16.7748
C        13.3031        30.0705        19.4582
C        11.7658        28.3352        20.7999
C        11.7658        28.3352        23.4833
C        10.1337        21.7134        16.7748
C        10.6885        26.2825        19.4582
C        10.1337        24.0316        20.7999
C        10.1337        24.0316        23.4833
C        11.7658        17.4098        16.7748
C        10.1337        21.7134        19.4582
C        10.6885        19.4625        20.7999
C        10.6885        19.4625        23.4833
C        15.211        14.3576        16.7748
C        11.7658        17.4098        19.4582
C        13.3031        15.6746        20.7999
C        13.3031        15.6746        23.4833
C        19.6799        13.2561        16.7748
C        15.211        14.3576        19.4582
C        17.3786        13.5356        20.7999
C        17.3786        13.5356        23.4833
C        24.1489        14.3576        16.7748
C        19.6799        13.2561        19.4582
C        21.9813        13.5356        20.7999
C        21.9813        13.5356        23.4833
C        27.594        17.4098        16.7748
C        24.1489        14.3576        19.4582
C        26.0568        15.6746        20.7999
C        26.0568        15.6746        23.4833
C        29.2262        21.7134        16.7748
C        27.594        17.4098        19.4582
C        28.6714        19.4625        20.7999
C        28.6714        19.4625        23.4833
C        29.2262        21.7134        19.4582
C        28.6714        26.2825        24.825
C        29.2262        24.0316        28.8501
C        29.2262        24.0316        31.5335
C        26.0568        30.0705        24.825
C        28.6714        26.2825        27.5084
C        27.594        28.3352        28.8501
C        27.594        28.3352        31.5335
C        21.9813        32.2094        24.825
C        26.0568        30.0705        27.5084
C        24.1489        31.3874        28.8501
C        24.1489        31.3874        31.5335
C        17.3786        32.2094        24.825
C        21.9813        32.2094        27.5084
C        19.6799        32.4889        28.8501
C        19.6799        32.4889        31.5335
C        13.3031        30.0705        24.825
C        17.3786        32.2094        27.5084
C        15.211        31.3874        28.8501
C        15.211        31.3874        31.5335
C        10.6885        26.2825        24.825
C        13.3031        30.0705        27.5084
C        11.7658        28.3352        28.8501
C        11.7658        28.3352        31.5335
C        10.1337        21.7134        24.825
C        10.6885        26.2825        27.5084
C        10.1337        24.0316        28.8501
C        10.1337        24.0316        31.5335
C        11.7658        17.4098        24.825
C        10.1337        21.7134        27.5084
C        10.6885        19.4625        28.8501
C        10.6885        19.4625        31.5335
C        15.211        14.3576        24.825
C        11.7658        17.4098        27.5084
C        13.3031        15.6746        28.8501
C        13.3031        15.6746        31.5335
C        19.6799        13.2561        24.825
C        15.211        14.3576        27.5084
C        17.3786        13.5356        28.8501
C        17.3786        13.5356        31.5335
C        24.1489        14.3576        24.825
C        19.6799        13.2561        27.5084
C        21.9813        13.5356        28.8501
C        21.9813        13.5356        31.5335
C        27.594        17.4098        24.825
C        24.1489        14.3576        27.5084
C        26.0568        15.6746        28.8501
C        26.0568        15.6746        31.5335
C        29.2262        21.7134        24.825
C        27.594        17.4098        27.5084
C        28.6714        19.4625        28.8501
C        28.6714        19.4625        31.5335
C        29.2262        21.7134        27.5084
Au        22.3609        34.5351        17.2871
Au        25.1602        37.8218        17.9816
Au        21.3376        38.6366        18.3464
Au        23.0098        37.9754        15.0892
Au        23.5247        41.8599        12.567
Au        19.6975        42.3466        19.4269
Au        28.4041        41.327        20.0947
Au        23.9107        42.0014        20.0721
Au        21.1661        42.1764        16.0619
Au        25.7285        41.9286        16.0269
K_POINTS {Gamma}

最后的结果很奇怪：

     Parallel version (MPI)

     Number of processors in use:      16
     R & G space division:  proc/pool =   16

     For Norm-Conserving or Ultrasoft (Vanderbilt) Pseudopotentials or PAW

     Current dimensions of program pwscf are:
     Max number of different atomic species (ntypx) = 10
     Max number of k-points (npk) =  40000
     Max angular momentum in pseudopotentials (lmaxx) =  3
     Waiting for input...
     file C.pbe-rrkjus.UPF: wavefunction(s)  2S 2P renormalized

     gamma-point specific algorithms are used


     Subspace diagonalization in iterative solution of the eigenvalue problem:
     a parallel distributed memory algorithm will be used,
     eigenstates matrixes will be distributed block like on
     ortho sub-group =    4*   4 procs


     Planes per process (thick) : nr3 =200 npp =  13 ncplane =60750
     Planes per process (smooth): nr3s=125 npps=   8 ncplanes=25920


     bravais-lattice index     =            8
     lattice parameter (a_0)   =      39.6842  a.u.
     unit-cell volume          =   65304.2406 (a.u.)^3
     number of atoms/cell      =          217
     number of atomic types    =            2
     number of electrons       =       938.00
     number of Kohn-Sham states=          563
     kinetic-energy cutoff     =      30.0000  Ry
     charge density cutoff     =     300.0000  Ry
     convergence threshold     =      1.0E-07
     mixing beta               =       0.7000
     number of iterations used =            8  plain     mixing
     Exchange-correlation      =  SLA  PW   PBE  PBE (1434)

     celldm(1)=  39.684200  celldm(2)=   1.190400  celldm(3)=   0.877800
     celldm(4)=   0.000000  celldm(5)=   0.000000  celldm(6)=   0.000000

     crystal axes: (cart. coord. in units of a_0)
               a(1) = (  1.000000  0.000000  0.000000 )  
               a(2) = (  0.000000  1.190400  0.000000 )  
               a(3) = (  0.000000  0.000000  0.877800 )  

     reciprocal axes: (cart. coord. in units 2 pi/a_0)
               b(1) = (  1.000000  0.000000  0.000000 )  
               b(2) = (  0.000000  0.840054  0.000000 )  
               b(3) = (  0.000000  0.000000  1.139212 )  


     PseudoPot. # 1 for C  read from file C.pbe-rrkjus.UPF
     Pseudo is Ultrasoft, Zval =  4.0
     Generated by new atomic code, or converted to UPF format
     Using radial grid of  627 points,  4 beta functions with:
                l(1) =   0
                l(2) =   0
                l(3) =   1
                l(4) =   1
     Q(r) pseudized with 0 coefficients


     PseudoPot. # 2 for Au read from file Au.pbe-nd-rrkjus.UPF
     Pseudo is Ultrasoft + core correction, Zval = 11.0
     Generated by new atomic code, or converted to UPF format
     Using radial grid of 1279 points,  3 beta functions with:
                l(1) =   2
                l(2) =   2
                l(3) =   1
     Q(r) pseudized with 0 coefficients


     atomic species   valence    mass     pseudopotential
        C              4.00    12.01070     C ( 1.00)
        Au            11.00   196.96600     Au( 1.00)

     ...............
  Self-consistent Calculation

     iteration #  1     ecut=    30.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  1.00E-02,  avg # of iterations =  2.0

     negative rho (up, down):  0.993E-03 0.000E+00

     total cpu time spent up to now is   1113.87 secs

     total energy              =   -2977.22341292 Ry
     Harris-Foulkes estimate   =   -3236.40939116 Ry
     estimated scf accuracy    <     294.26077390 Ry

     iteration #  2     ecut=    30.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  1.00E-02,  avg # of iterations = 12.0

     negative rho (up, down):  0.117E+00 0.000E+00

     total cpu time spent up to now is   3979.61 secs

     total energy              =    -497.62887573 Ry
     Harris-Foulkes estimate   =   -6586.82676250 Ry
     estimated scf accuracy    <   87542.15719839 Ry

     iteration #  3     ecut=    30.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  1.00E-02,  avg # of iterations = 11.0

     negative rho (up, down):  0.697E-01 0.000E+00

     total cpu time spent up to now is   6846.40 secs

     total energy              =   -3156.18188028 Ry
     Harris-Foulkes estimate   =   -3287.62595780 Ry
     estimated scf accuracy    <     875.21197632 Ry

     iteration #  4     ecut=    30.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  1.00E-02,  avg # of iterations =  6.0

     negative rho (up, down):  0.253E-01 0.000E+00

     total cpu time spent up to now is   7879.45 secs

     total energy              =   -3203.59796421 Ry
     Harris-Foulkes estimate   =   -3241.65719464 Ry
     estimated scf accuracy    <     224.69123842 Ry

     iteration #  5     ecut=    30.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  1.00E-02,  avg # of iterations =  4.0

     negative rho (up, down):  0.200E-01 0.000E+00

     total cpu time spent up to now is   8442.13 secs

     total energy              =   -3215.43269596 Ry
     Harris-Foulkes estimate   =   -3229.12440698 Ry
     estimated scf accuracy    <     103.54941645 Ry

     iteration #  6     ecut=    30.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  1.00E-02,  avg # of iterations =  2.0

     negative rho (up, down):  0.157E-01 0.000E+00

     total cpu time spent up to now is   8955.15 secs

     total energy              =   -3217.71950617 Ry
     Harris-Foulkes estimate   =   -3227.87637503 Ry
     estimated scf accuracy    <     104.30799301 Ry

     iteration #  7     ecut=    30.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  1.00E-02,  avg # of iterations =  2.0

     negative rho (up, down):  0.108E-01 0.000E+00

     total cpu time spent up to now is   9481.27 secs

     total energy              =   -3218.72993801 Ry
     Harris-Foulkes estimate   =   -3224.43398830 Ry
     estimated scf accuracy    <      40.06120485 Ry

     iteration #  8     ecut=    30.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  4.27E-03,  avg # of iterations =  2.0

     negative rho (up, down):  0.121E-01 0.000E+00

     total cpu time spent up to now is   9990.08 secs

     total energy              =   -3219.49626617 Ry
     Harris-Foulkes estimate   =   -3224.58195503 Ry
     estimated scf accuracy    <      41.17142700 Ry

     iteration #  9     ecut=    30.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  4.27E-03,  avg # of iterations =  2.0

     negative rho (up, down):  0.485E-01 0.000E+00

     total cpu time spent up to now is  10515.47 secs

     total energy              =   -3219.48018063 Ry
     Harris-Foulkes estimate   =   -3223.40913785 Ry
     estimated scf accuracy    <      31.70397119 Ry

     iteration # 10     ecut=    30.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  3.38E-03,  avg # of iterations =  1.0

     negative rho (up, down):  0.143E+00 0.000E+00

     total cpu time spent up to now is  11035.10 secs

     total energy              =   -3219.58368562 Ry
     Harris-Foulkes estimate   =   -3224.09297601 Ry
     estimated scf accuracy    <      58.18087981 Ry

     iteration # 11     ecut=    30.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  3.38E-03,  avg # of iterations =  1.0

     negative rho (up, down):  0.715E-01 0.000E+00

     total cpu time spent up to now is  11549.73 secs

     total energy              =   -3220.88759528 Ry
     Harris-Foulkes estimate   =   -3222.28928438 Ry
     estimated scf accuracy    <      10.40542082 Ry

     iteration # 12     ecut=    30.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  1.11E-03,  avg # of iterations =  5.0

     negative rho (up, down):  0.112E-01 0.000E+00

     total cpu time spent up to now is  12101.65 secs

     total energy              =   -3220.64662380 Ry
     Harris-Foulkes estimate   =   -3222.30862685 Ry
     estimated scf accuracy    <      13.50569556 Ry

     iteration # 13     ecut=    30.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  1.11E-03,  avg # of iterations =  6.0

     negative rho (up, down):  0.132E-01 0.000E+00

     total cpu time spent up to now is  12649.82 secs

     total energy              =   -3221.07090471 Ry
     Harris-Foulkes estimate   =   -3222.30118620 Ry
     estimated scf accuracy    <       9.17029961 Ry

     iteration # 14     ecut=    30.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  9.78E-04,  avg # of iterations =  3.0

     negative rho (up, down):  0.981E-02 0.000E+00

     total cpu time spent up to now is  13178.03 secs

     total energy              =   -3221.34879009 Ry
     Harris-Foulkes estimate   =   -3222.07646819 Ry
     estimated scf accuracy    <       5.92444643 Ry

     iteration # 15     ecut=    30.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  6.32E-04,  avg # of iterations =  1.0

     negative rho (up, down):  0.625E-02 0.000E+00

     total cpu time spent up to now is  13703.17 secs

     total energy              =   -3221.56124590 Ry
     Harris-Foulkes estimate   =   -3221.83904556 Ry
     estimated scf accuracy    <       3.21641224 Ry
..............
（省略中间迭代结果）
.............
iteration # 96     ecut=    30.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     WARNING:     1 eigenvalues not converged
     c_bands:  1 eigenvalues not converged
     ethr =  5.53E-06,  avg # of iterations = 20.0

     negative rho (up, down):  0.208E-01 0.000E+00

     total cpu time spent up to now is  90461.17 secs

     total energy              =   -3217.77281804 Ry
     Harris-Foulkes estimate   =   -3226.66977913 Ry
     estimated scf accuracy    <      81.38309144 Ry

     iteration # 97     ecut=    30.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     WARNING:     1 eigenvalues not converged
     c_bands:  1 eigenvalues not converged
     ethr =  5.53E-06,  avg # of iterations = 20.0

     negative rho (up, down):  0.299E-01 0.000E+00

     total cpu time spent up to now is  92155.74 secs

     total energy              =   -3211.79720560 Ry
     Harris-Foulkes estimate   =   -3230.89684691 Ry
     estimated scf accuracy    <     267.53665710 Ry

     iteration # 98     ecut=    30.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  5.53E-06,  avg # of iterations = 18.0

     negative rho (up, down):  0.251E-01 0.000E+00

     total cpu time spent up to now is  93694.72 secs

     total energy              =   -3217.68647998 Ry
     Harris-Foulkes estimate   =   -3226.58151504 Ry
     estimated scf accuracy    <      78.81725327 Ry

     iteration # 99     ecut=    30.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  5.53E-06,  avg # of iterations = 18.0

     negative rho (up, down):  0.299E-01 0.000E+00

     total cpu time spent up to now is  95407.43 secs

     total energy              =   -3214.38628245 Ry
     Harris-Foulkes estimate   =   -3232.15352793 Ry
     estimated scf accuracy    <     187.19089394 Ry

     iteration #100     ecut=    30.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  5.53E-06,  avg # of iterations = 19.0

     negative rho (up, down):  0.155E-01 0.000E+00

     total cpu time spent up to now is  97155.45 secs

     total energy              =   -3219.75593714 Ry
     Harris-Foulkes estimate   =   -3223.84708261 Ry
     estimated scf accuracy    <      30.20838506 Ry

     End of self-consistent calculation

     convergence NOT achieved after 100 iterations: stopping
.
请教高手，我这个ecut应该怎么设置？

[ Last edited by frank_zhan on 2012-7-1 at 17:06 ]
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐，对您同样有帮助:
frank_zhan,zy2zhan@gmail.com
1楼 2012-07-01 12:45:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

frank_zhan

金虫 (正式写手)
谁能帮个忙？
回复此楼
frank_zhan,zy2zhan@gmail.com
2楼2012-07-01 17:07:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

szdtzjz

铁虫 (初入文坛)

【答案】应助回帖

★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
frank_zhan: 金币+5, ★有帮助, 谢谢回复，我先试试。 2012-07-02 09:48:49
WDD880227: 金币+1, 多谢指教 2012-07-02 10:19:41
减小mixing_beta的值，取0.1~0.3试试看，并且可以把electron_maxstep改小一点，比如20步，先看看收敛的情况怎样。如果好的话，就把restart_mode改成“restart”让它接着算；如果收敛太慢，就提高一下mixing_beta，再restart；如果estimated scf accuracy震荡还是很严重，就再减小一下electron_maxstep，再restart。
如果还是不行，可以用nbnd，增加几条空能带，也会好一点。
这个收敛的问题，在PWscf的user guide里面的Frequently Asked Questions一章有讲，你可以看看。
赞一下 回复此楼
3楼2012-07-01 19:11:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

frank_zhan

金虫 (正式写手)
引用回帖:
3楼: Originally posted by szdtzjz at 2012-07-01 19:11:54
减小mixing_beta的值，取0.1~0.3试试看，并且可以把electron_maxstep改小一点，比如20步，先看看收敛的情况怎样。如果好的话，就把restart_mode改成“restart”让它接着算；如果收敛太慢，就提高一下mixing_beta，再 ...

谢谢你的建议，你认为ecutwfc需不需要再大一些？
赞一下 回复此楼
frank_zhan,zy2zhan@gmail.com
4楼2012-07-02 09:49:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

szdtzjz

铁虫 (初入文坛)

【答案】应助回帖

★ ★ ★ ★ ★ ★
frank_zhan: 金币+5, ★有帮助, 谢谢交流。 2012-07-02 14:13:13
franch: 金币+1, 谢谢回帖交流 2012-07-02 15:25:32
引用回帖:
4楼: Originally posted by frank_zhan at 2012-07-02 09:49:41
谢谢你的建议，你认为ecutwfc需不需要再大一些？...

不需要吧……老师讲ecutwfc一般20~30就足够了。。如果你对这个值不放心可以这样：
先取一个稍低一点的值，先让计算收敛，然后再提高这个值再算，比较计算的总能量。如果随着ecutwfc的提高，总能量的波动较小（mRy/每个原子），就说明ecutwfc已经足够高了
赞一下 回复此楼
5楼2012-07-02 13:44:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

frank_zhan

金虫 (正式写手)
还有是10个金币呢，不知道哪位高人再来指导一下？
赞一下 回复此楼
frank_zhan,zy2zhan@gmail.com
6楼2012-07-05 09:51:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

goldenfisher

金虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
liliangfang: 金币+1, 谢谢交流 2012-07-05 14:55:26
frank_zhan: 金币+5, ★有帮助, 谢谢交流。 2012-07-07 14:43:42
frank_zhan: 金币+5, ★★★很有帮助 2012-07-31 19:46:00
个人认为两个问题：
一是体系有金在，而且你还用的是包含d电子的PP，所以ecutwfc明显太小
而是你的smearing用的fd,这个在低温下本来就是收敛很慢的方法。可以改用mp或者mv。
赞一下 回复此楼
7楼2012-07-05 10:58:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

frank_zhan

金虫 (正式写手)
引用回帖:
7楼: Originally posted by goldenfisher at 2012-07-05 10:58:12
个人认为两个问题：
一是体系有金在，而且你还用的是包含d电子的PP，所以ecutwfc明显太小
而是你的smearing用的fd,这个在低温下本来就是收敛很慢的方法。可以改用mp或者mv。

多谢交流。那我再请教一下，怎样确定不同pp对ecutwfc值的要求呢？还有怎么才能知道那种pp是含有d电子的呢？
赞一下 回复此楼
frank_zhan,zy2zhan@gmail.com
8楼2012-07-07 14:43:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

goldenfisher

金虫 (著名写手)
引用回帖:
8楼: Originally posted by frank_zhan at 2012-07-07 14:43:21
多谢交流。那我再请教一下，怎样确定不同pp对ecutwfc值的要求呢？还有怎么才能知道那种pp是含有d电子的呢？...

要求你自己最好做收敛测试，ecut不够大的时候，收敛是很困难的，而且还有能量的收敛性也依赖于这个ecut
pp电子知道的方法很多很多，你直接打开pp文件看就知道了
赞一下 回复此楼
9楼2012-07-31 21:57:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huazhorg

铁杆木虫 (正式写手)
这个问题上次已经说过吧，你要增加electron_maxstep，默认就是100.
赞一下 回复此楼
10楼2012-08-01 16:28:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 frank_zhan 的主题更新
141/2返回列表12下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭