24小时热门版块排行榜    

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

1164028932

金虫 (初入文坛)

[求助] materials-studio软件计算氧气分子能量。 已有3人参与

步骤如下:
      1、先构建一个3D Atomistic文件(即构建一个3D物件视窗);2、在Build——>Crystals——>Build Crystals,这样就构建出了一个立方晶格空间。这个你肯定会。3、Build——>Add Atoms——>选择O原子,在把xyz坐标改为0.5、0.5、0.5,这样就在立方晶格中央构建出了一个O原子,你再用添加原子的笔在那个O原子上加上去一个O原子,就构建出了可以计算氧分子能量的模型。4、单独计算O原子的话,不再添加第二个O原子就行。

计算结果:我用这种方法算了氧气分子能量。castep文件中final energy=-866eV,是看这个数据吗?还是看其他的。谢谢了(看别人计算的,还有文献上给的能量大约都是8~9ev)
   

castep文件:

*********************************** Forces ***********************************
*                                                                            *
*                        Cartesian components (eV/A)                         *
* -------------------------------------------------------------------------- *
*                   x                    y                    z              *
*                                                                            *
* O         1     -0.00085              0.00568             -0.00079         *
* O         2      0.00085             -0.00568              0.00079         *
*                                                                            *
******************************************************************************

***************** Stress Tensor *****************
*                                               *
*          Cartesian components (GPa)           *
* --------------------------------------------- *
*             x             y             z     *
*                                               *
*  x     -0.018289      0.002225     -0.000730  *
*  y      0.002225     -0.006909     -0.006006  *
*  z     -0.000730     -0.006006     -0.013130  *
*                                               *
*  Pressure:    0.0128                          *
*                                               *
*************************************************
  
+------------+-------------+-------------+-----------------+ <-- min BFGS
|    Step    |   lambda    |   F.delta   |    enthalpy     | <-- min BFGS
+------------+-------------+-------------+-----------------+ <-- min BFGS
|  previous  |    0.000000 |    0.000000 |     -866.417494 | <-- min BFGS
| trial step |    1.000000 |    0.000000 |     -866.417497 | <-- min BFGS
+------------+-------------+-------------+-----------------+ <-- min BFGS
  
BFGS: finished iteration     4 with enthalpy= -8.66417497E+002 eV
  
+-----------+-----------------+-----------------+------------+-----+ <-- BFGS
| Parameter |      value      |    tolerance    |    units   | OK? | <-- BFGS
+-----------+-----------------+-----------------+------------+-----+ <-- BFGS
|  dE/ion   |   1.577937E-006 |   5.000000E-006 |         eV | Yes | <-- BFGS
|  |F|max   |   5.800782E-003 |   1.000000E-002 |       eV/A | Yes | <-- BFGS
|  |dR|max  |   3.427978E-005 |   5.000000E-004 |          A | Yes | <-- BFGS
|   Smax    |   1.828935E-002 |   2.000000E-002 |        GPa | Yes | <-- BFGS
+-----------+-----------------+-----------------+------------+-----+ <-- BFGS
  
BFGS: Geometry optimization completed successfully.

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

                           -------------------------------
                                      Unit Cell
                           -------------------------------
        Real Lattice(A)              Reciprocal Lattice(1/A)
  10.0007554   0.0004342   0.0009489        0.6282711  -0.0000273  -0.0000596
   0.0004342  10.0010894   0.0018433       -0.0000273   0.6282501  -0.0001158
   0.0009489   0.0018433  10.0034738       -0.0000596  -0.0001158   0.6281004

                       Lattice parameters(A)       Cell Angles
                    a =   10.000755          alpha =   89.978882
                    b =   10.001090          beta  =   89.989128
                    c =   10.003474          gamma =   89.995024

                Current cell volume = 1000.531886 A**3

                           -------------------------------
                                     Cell Contents
                           -------------------------------

            xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
            x  Element    Atom        Fractional coordinates of atoms  x
            x            Number           u          v          w      x
            x----------------------------------------------------------x
            x  O            1         0.515446   0.528133   0.552192   x
            x  O            2        -0.515446  -0.528133  -0.552192   x
            xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx


BFGS: Final Enthalpy     = -8.66417497E+002 eV
BFGS: Final <frequency>  =       1.14214 cm-1
BFGS: Final bulk modulus unchanged from initial value
materials-studio软件计算氧气分子能量。
结构图
回复此楼

» 猜你喜欢

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

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

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

1164028932

金虫 (初入文坛)

引用回帖:
2楼: Originally posted by mywai520 at 2014-04-14 22:50:11
楼主看的没错,就是final energy

那看别人计算的,还有文献上给的能量大约都是8~9ev,怎么差距这么大呀?方法有问题吗?
3楼2014-04-15 10:37:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答

mywai520

铁杆木虫 (著名写手)


【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
1164028932: 金币+5, ★★★很有帮助 2014-04-15 10:38:15
楼主看的没错,就是final energy
2楼2014-04-14 22:50:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qlsy

银虫 (小有名气)

困惑中

[ 发自手机版 http://muchong.com/3g ]
4楼2014-04-23 15:26:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小小的西瓜籽

金虫 (著名写手)

【答案】应助回帖

我是直接 画两个氧原子出来计算的,然后几何优化就得到结果是-150.3956335Ha,DMol3 中算的。
5楼2014-04-24 11:36:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +4 8rmuugja8q 2026-02-22 8/400 2026-02-23 12:22 by alian_214
[基金申请] 体制内长辈说体制内绝大部分一辈子在底层,如同你们一样大部分普通教师忙且收入低 +10 瞬息宇宙 2026-02-20 13/650 2026-02-23 11:23 by holypower
[考研] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 khieu8v8m0 2026-02-22 8/400 2026-02-23 09:35 by w4l55oybr1
[论文投稿] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +4 khieu8v8m0 2026-02-22 8/400 2026-02-23 09:29 by w4l55oybr1
[考研] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +5 usprnugpzw 2026-02-21 11/550 2026-02-23 09:24 by w4l55oybr1
[教师之家] 为什么中国大学工科教授们水了那么多所谓的顶会顶刊,但还是做不出宇树机器人? +5 欢乐颂叶蓁 2026-02-21 8/400 2026-02-23 09:19 by 欢乐颂叶蓁
[论文投稿] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 w89i99eaeh 2026-02-22 5/250 2026-02-23 08:04 by w4l55oybr1
[博后之家] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +4 khieu8v8m0 2026-02-22 6/300 2026-02-23 07:59 by w4l55oybr1
[博后之家] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +6 3dfhjxgsh7 2026-02-22 9/450 2026-02-23 07:49 by w4l55oybr1
[考博] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +4 khieu8v8m0 2026-02-22 4/200 2026-02-23 06:46 by jsjzfl
[公派出国] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 khieu8v8m0 2026-02-22 5/250 2026-02-23 06:29 by w4l55oybr1
[硕博家园] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +4 khieu8v8m0 2026-02-22 8/400 2026-02-23 06:24 by w4l55oybr1
[考博] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +5 3dfhjxgsh7 2026-02-22 6/300 2026-02-23 02:04 by 5jlh3qtdvx
[教师之家] 版面费该交吗 +7 苹果在哪里 2026-02-22 8/400 2026-02-22 22:37 by otani
[论文投稿] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +4 usprnugpzw 2026-02-21 6/300 2026-02-22 19:48 by w89i99eaeh
[考研] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 3dfhjxgsh7 2026-02-22 4/200 2026-02-22 16:52 by khieu8v8m0
[找工作] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 usprnugpzw 2026-02-22 3/150 2026-02-22 16:37 by khieu8v8m0
[公派出国] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 usprnugpzw 2026-02-21 4/200 2026-02-22 16:27 by khieu8v8m0
[基金申请] “人文社科而论,许多学术研究还没有达到民国时期的水平” +4 苏东坡二世 2026-02-18 5/250 2026-02-22 16:07 by liangep1573
[基金申请] 什么是人一生最重要的? +4 瞬息宇宙 2026-02-21 4/200 2026-02-22 11:44 by huagongfeihu
信息提示
请填处理意见