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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » MS » 使用castep从哪儿看原子的能量
61/1返回列表
查看: 2390  |  回复: 5
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

dianfu

铜虫 (小有名气)

[求助] 使用castep从哪儿看原子的能量 已有1人参与

我做掺杂，要算结合能，需要知道原来原子与掺杂原子的能量，请问从哪儿可以看到？谢谢了

[ 发自手机版 http://muchong.com/3g ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
正在学习量子计算
1楼 2015-05-18 16:04:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pcczhouwei

铁杆木虫 (著名写手)

Wolf

【答案】应助回帖

★ ★
感谢参与，应助指数 +1
franch: 金币+2, 谢谢回帖交流。。。。 2015-05-20 11:17:49
计算完的体系，在CASTEP中倒着往上看，看到Geometry optimization success 后，上面的能量值就是所求出的能量值。
赞一下 回复此楼
知與誰同，與狼共舞。
2楼2015-05-18 16:32:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dianfu

铜虫 (小有名气)
引用回帖:
2楼: Originally posted by pcczhouwei at 2015-05-18 16:32:46
计算完的体系，在CASTEP中倒着往上看，看到Geometry optimization success 后，上面的能量值就是所求出的能量值。

Pseudo atomic calculation performed for O 2s2 2p4

Converged in 22 iterations to a total energy of -429.1638 eV


Pseudo atomic calculation performed for Ti 3s2 3p6 3d2 4s2

Converged in 32 iterations to a total energy of -1595.8791 eV


Pseudo atomic calculation performed for La 5s2 5p6 5d1 6s2

Converged in 22 iterations to a total energy of -857.7835 eV


Pseudo atomic calculation performed for Sm 4f6 5s2 5p6 6s2

Converged in 43 iterations to a total energy of -2198.6961 eV
是类似这样的数据吗？这样的话是不是意味着体系中相同原子的能量都是相同的？

[ 发自小木虫客户端 ]
赞一下 回复此楼
正在学习量子计算
3楼2015-05-18 22:06:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dianfu

铜虫 (小有名气)
引用回帖:
2楼: Originally posted by pcczhouwei at 2015-05-18 16:32:46
计算完的体系，在CASTEP中倒着往上看，看到Geometry optimization success 后，上面的能量值就是所求出的能量值。

还有一点想请教，那就是几个结构优化后的castep文件和能量计算后的castep文件中均有各原子的能量和系统能量，但是两个文件中的值不太一样，我应该选择哪一个呢？

[ 发自小木虫客户端 ]
赞一下 回复此楼
正在学习量子计算
4楼2015-05-18 22:36:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pcczhouwei

铁杆木虫 (著名写手)

Wolf

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
franch: 金币+3, 谢谢回帖交流。。。。 2015-05-20 11:18:54
dianfu: 金币+10, ★★★★★最佳答案 2015-05-20 16:26:56
引用回帖:
4楼: Originally posted by dianfu at 2015-05-18 22:36:46
还有一点想请教，那就是几个结构优化后的castep文件和能量计算后的castep文件中均有各原子的能量和系统能量，但是两个文件中的值不太一样，我应该选择哪一个呢？
...

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

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


*********************************** Forces ***********************************
*                                                                            *
*                        Cartesian components (eV/A)                         *
* -------------------------------------------------------------------------- *
*                   x                    y                    z              *
*                                                                            *
* H         1      0.04136              0.00324              0.01716         *
* H         2      0.00705              0.01903             -0.01014         *
* O         1      0.00466             -0.02132              0.00840         *
* O         2      0.00292              0.02071              0.01710         *
* O         3     -0.00531             -0.00583             -0.00846         *
* O         4      0.00095             -0.00131             -0.01021         *
* O         5     -0.00829              0.01346              0.00232         *
* O         6     -0.01514             -0.02334              0.00166         *
* O         7     -0.00668             -0.00023             -0.00191         *
* O         8      0.00167              0.00562             -0.00922         *
* O         9     -0.01073             -0.00667             -0.00483         *
* O        10      0.00630             -0.00267             -0.00738         *
* O        11     -0.00834              0.01010              0.00331         *
* O        12     -0.02338             -0.00577             -0.00092         *
* O        13      0.00619             -0.01272              0.00416         *
* O        14      0.00581              0.00571             -0.00048         *
* O        15     -0.00903             -0.00052             -0.00607         *
* O        16     -0.00424              0.00311              0.02818         *
* O        17     -0.00932              0.00094             -0.00312         *
* O        18     -0.02026             -0.00632              0.00027         *
* O        19     -0.00256              0.00048             -0.00037         *
* O        20      0.00422              0.00565             -0.00288         *
* O        21      0.01033              0.00224             -0.00273         *
* O        22      0.00041              0.00857             -0.00177         *
* O        23      0.00207             -0.00600             -0.00923         *
* O        24     -0.00350             -0.00980             -0.00061         *
* O        25      0.01560             -0.00660             -0.00235         *
* O        26      0.01558             -0.00428             -0.01675         *
* O        27     -0.00090             -0.00603              0.00304         *
* O        28      0.01362              0.00437              0.00915         *
* O        29      0.00625              0.01304              0.00034         *
* O        30      0.01011             -0.01484              0.01123         *
* O        31     -0.00952              0.01153              0.00754         *
* O        32     -0.00823             -0.01031             -0.00893         *
* O        33     -0.00928             -0.00459              0.00166         *
* O        34      0.00293             -0.00721             -0.01687         *
* O        35      0.00510             -0.00156             -0.01231         *
* O        36     -0.00074             -0.00280              0.00708         *
* O        37      0.00786             -0.00830              0.00688         *
* O        38      0.01321             -0.00795             -0.00778         *
* O        39      0.00372             -0.00254             -0.00270         *
* O        40      0.00360             -0.01011             -0.00962         *
* O        41      0.01682              0.00310             -0.00225         *
* O        42      0.01877              0.00020              0.00951         *
* O        43      0.00057              0.00336              0.00403         *
* O        44      0.00408             -0.00998              0.00327         *
* O        45      0.00297              0.00137              0.00803         *
* O        46     -0.00574             -0.00578             -0.00831         *
* O        47      0.00237              0.01414             -0.00319         *
* O        48      0.00038              0.00458             -0.00005         *
* Fe        1      0.01346              0.00453             -0.00795         *
* Fe        2      0.00983             -0.00750              0.01596         *
* Fe        3      0.01550             -0.01329             -0.00333         *
* Fe        4      0.00184              0.01390             -0.00645         *
* Fe        5      0.00279              0.00643             -0.00287         *
* Fe        6     -0.00487              0.00441              0.00494         *
* Fe        7      0.00475             -0.01923             -0.00396         *
* Fe        8     -0.00241             -0.00710             -0.00139         *
* Fe        9     -0.00620              0.01011              0.00789         *
* Fe       10     -0.00678              0.01036             -0.01054         *
* Fe       11     -0.00809              0.01176              0.01638         *
* Fe       12      0.00777             -0.00306              0.00958         *
* Fe       13     -0.01879             -0.00149             -0.00699         *
* Fe       14     -0.00890              0.00853             -0.01689         *
* Fe       15     -0.01486              0.00166              0.00669         *
* Fe       16      0.00080             -0.00442              0.01152         *
* La        1     -0.00010             -0.01021              0.00185         *
* La        2     -0.00673              0.01168              0.01778         *
* La        3     -0.01831              0.03447             -0.02515         *
* La        4      0.00607              0.00059             -0.00288         *
* La        5      0.00709             -0.00082              0.00329         *
* La        6     -0.00532              0.00117             -0.01249         *
* La        7      0.00507             -0.00118              0.00604         *
* La        8     -0.01192             -0.00091             -0.01178         *
* La        9     -0.01352             -0.00035             -0.00097         *
* La       10      0.00531              0.01436             -0.00365         *
* La       11      0.00189              0.00021             -0.00672         *
* La       12     -0.01841              0.00791              0.00715         *
* La       13     -0.00969             -0.01170              0.00391         *
* La       14      0.01087             -0.01024              0.02715         *
* La       15     -0.01248             -0.01038             -0.01199         *
* La       16     -0.00199              0.01070              0.01299         *
*                                                                            *
******************************************************************************
  
+------------+-------------+-------------+-----------------+ <-- min BFGS
|    Step    |   lambda    |   F.delta   |    enthalpy     | <-- min BFGS
+------------+-------------+-------------+-----------------+ <-- min BFGS
|  previous  |    0.000000 |    0.000012 |   -48694.478532 | <-- min BFGS
| trial step |    1.000000 |    0.000009 |   -48694.478806 | <-- min BFGS
+------------+-------------+-------------+-----------------+ <-- min BFGS
  
BFGS: finished iteration   125 with enthalpy= -4.86944788E+004 eV
  
+-----------+-----------------+-----------------+------------+-----+ <-- BFGS
| Parameter |      value      |    tolerance    |    units   | OK? | <-- BFGS
+-----------+-----------------+-----------------+------------+-----+ <-- BFGS
|  dE/ion   |   3.343676E-006 |   2.000000E-005 |         eV | Yes | <-- BFGS
|  |F|max   |   4.642571E-002 |   5.000000E-002 |       eV/A | Yes | <-- BFGS
|  |dR|max  |   8.772360E-004 |   2.000000E-003 |          A | Yes | <-- BFGS
+-----------+-----------------+-----------------+------------+-----+ <-- BFGS
  
BFGS : Geometry optimization completed successfully.

================================================================================
BFGS : Final Configuration:


我们一般取的就是这个 E  =  -48694.01692607     eV   就是我先说的几何优化成功上面的这个能量
赞一下(1人) 回复此楼
知與誰同，與狼共舞。
5楼2015-05-19 23:35:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dianfu

铜虫 (小有名气)
引用回帖:
5楼: Originally posted by pcczhouwei at 2015-05-19 23:35:01
Final energy, E             =  -48694.01692607     eV
Final free energy (E-TS)    =  -48694.47880611     eV
(energies not corrected for finite basis set)

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

好的，我懂了，谢谢你的耐心帮助。

[ 发自小木虫客户端 ]
赞一下 回复此楼
正在学习量子计算
6楼2015-05-20 16:26:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 dianfu 的主题更新
61/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭