ms和高斯计算氮化碳对氧气的吸附能求助

成都张雨绮

高斯计算模型 输入文件如下 （氧气是开壳的计算也考虑了 结果差异不大 就没用）
%chk=CN-O2-4.chk
%mem=8GB
%nprocshared=8
# opt=readoptimize freq b3lyp/6-31g(d)
empiricaldispersion=gd3

Title Card Required

0 1
C                 -8.41898357    1.06489234   -0.46417751
C                -10.26127157    2.54977634   -0.46437151
C                -10.62466757    0.24636334   -0.46508551
C                -10.99037657   -2.09119066   -0.46539651
C                 -8.78457357   -1.23866066   -0.46460551
C                 -6.55188957   -0.41144966   -0.46346851
N                -11.09226757    1.53229834   -0.46491651
N                 -8.91606157    2.31590334   -0.46405051
N                 -9.27626257    0.02143934   -0.46466551
N                -11.48316057   -0.81563166   -0.46545751
N                 -9.64037957   -2.27789866   -0.46495051
N                 -7.44147457   -1.44943066   -0.46403251
N                 -7.07741957    0.85047934   -0.46360451
N                 -5.16391657   -0.63012266   -0.46276351
C                 -1.68307557    3.65864134   -0.45774751
C                 -3.54033457    5.12286934   -0.45778651
C                 -3.87840157    2.81530834   -0.45966251
C                 -4.28099557    0.46539534   -0.46161151
C                 -2.06712657    1.34107634   -0.45966151
C                  0.12505043    2.13271634   -0.45778651
N                 -4.34215257    4.08478334   -0.45900851
N                 -2.18450957    4.92692234   -0.45711751
N                 -2.54407757    2.60402534   -0.45903251
N                 -4.73490357    1.75778634   -0.46092651
N                 -2.92526157    0.28889234   -0.46091551
N                 -0.73292957    1.14187334   -0.45901051
N                 -0.33903557    3.42200434   -0.45712051
N                  1.47550043    1.88157234   -0.45715951
C                 -3.16261557   -5.75372766   -0.46360551
C                 -0.95080657   -4.89068066   -0.46333851
C                 -2.80553957   -3.44796266   -0.46338951
C                 -4.65634057   -1.94291966   -0.46309851
C                 -4.99611957   -4.30569466   -0.46349251
C                 -5.36854757   -6.62476766   -0.46381851
N                 -1.46127357   -3.67513166   -0.46330751
N                 -1.82407857   -5.94361066   -0.46349751
N                 -3.65644257   -4.49878266   -0.46352351
N                 -3.30628357   -2.18538866   -0.46319851
N                 -5.48733357   -3.03624866   -0.46327751
N                 -5.82655357   -5.38864666   -0.46362051
N                 -4.01325657   -6.80331866   -0.46373351
C                  3.53123543   -3.14407266   -0.45959551
C                  5.72701443   -2.27645166   -0.45740851
C                  3.89772043   -0.82340166   -0.45759351
C                  2.02067543    0.62369334   -0.45778451
C                  1.70771243   -1.71231866   -0.45981851
C                  1.32205043   -4.01738266   -0.46181451
N                  5.25071643   -1.04813566   -0.45681251
N                  4.86144943   -3.33036966   -0.45878651
N                  3.04614943   -1.88876166   -0.45900651
N                  3.38044243    0.44301034   -0.45698051
N                  1.20888943   -0.46760166   -0.45920651
N                  0.87669843   -2.78506766   -0.46119351
N                  2.67863043   -4.19750866   -0.46097451
N                  0.43918143   -5.10721566   -0.46309351
C                  2.78983343   -7.92082666   -0.46452951
C                  0.94880943   -6.41747766   -0.46393951
C                  0.61825543   -8.73638866   -0.46543051
C                  0.21308243  -11.06318166   -0.46625951
C                  2.43879043  -10.24440866   -0.46565551
C                  4.62168643   -9.41083066   -0.46451451
N                  0.15163143   -7.45980166   -0.46476751
N                  2.29589543   -6.65831266   -0.46388751
N                  1.94784043   -8.97170766   -0.46526751
N                 -0.23375257   -9.77601566   -0.46593051
N                  1.56467043  -11.29692266   -0.46614051
N                  3.79661443  -10.43786066   -0.46526151
N                  4.11464543   -8.14483766   -0.46417851
C                 -9.54358057   -6.03084866   -0.46501951
C                -11.36613157   -4.53917966   -0.46511151
C                -11.75084957   -6.84144466   -0.46237051
C                -12.11713357   -9.18299566   -0.45934051
C                 -9.91440257   -8.33179666   -0.46143151
C                 -7.69288657   -7.50354666   -0.46368651
N                -12.20989357   -5.54926866   -0.46336951
N                -10.02846457   -4.78169766   -0.46614951
N                -10.40397857   -7.07114866   -0.46297351
N                -12.61199557   -7.90496766   -0.46059151
N                -10.76897657   -9.36825466   -0.45963951
N                 -8.57340657   -8.54541166   -0.46168151
N                 -8.20826157   -6.24444166   -0.46548651
N                 -6.23464957   -7.72471566   -0.46388351
N                -11.83658857   -3.16934666   -0.46555351
N                 -0.68266157  -12.10494466   -0.46624551
C                 -4.23645257  -12.85656966   -0.46330051
C                 -2.04751457  -11.94590866   -0.46571451
C                 -3.89699657  -10.53390166   -0.46524251
C                 -5.73435857   -9.03045866   -0.46400051
C                 -6.07237957  -11.38391766   -0.46241351
C                 -6.45702357  -13.68162866   -0.46038751
N                 -2.56687357  -10.74532466   -0.46620451
N                 -2.88064557  -13.03743066   -0.46432851
N                 -4.73607857  -11.59078866   -0.46370551
N                 -4.39189557   -9.27393066   -0.46540951
N                 -6.56595457  -10.11597466   -0.46258351
N                 -6.90200357  -12.44910066   -0.46076551
N                 -5.10395257  -13.91145566   -0.46163351
H                 -0.30214957  -13.07871466   -0.46612651
H                -12.86796557   -2.99969766   -0.46556451
H                  2.12765043    2.69874334   -0.45612351
C                -10.63805457  -13.11677766   -0.45735951
C                -12.47015057  -11.63563266   -0.45677951
C                -12.82607757  -13.94071266   -0.45294051
C                -13.18091357  -16.26862266   -0.44893551
C                -10.97753257  -15.43818466   -0.45315351
C                 -8.77956757  -14.55898666   -0.45733951
H                -13.86126157  -17.10901966   -0.44675551
N                -13.29992857  -12.65342766   -0.45409551
N                -11.13278657  -11.87000666   -0.45853951
N                -11.48068357  -14.17018166   -0.45448651
N                -13.67576857  -14.99855366   -0.45016851
N                -11.83586757  -16.48886866   -0.45038051
N                 -9.62314557  -15.63463866   -0.45454951
N                 -9.30904757  -13.31152766   -0.45879151
N                 -7.34046657  -14.72990066   -0.45864451
N                -12.95550957  -10.26787166   -0.45777051
C                 -7.67928457    5.85948134   -0.46037151
C                 -5.47483557    6.67539134   -0.45799151
C                 -7.27276957    8.16015534   -0.45789451
C                 -9.09344657    9.65300334   -0.45776451
C                 -9.50378857    7.33394734   -0.46020251
C                 -9.86398657    4.99443634   -0.46256351
H                 -9.46847357   10.66715334   -0.45708251
N                 -5.92474857    7.90930334   -0.45715351
N                 -6.35623157    5.64539034   -0.45963251
N                 -8.15615357    7.12088734   -0.45949751
N                 -7.74854457    9.43068834   -0.45703451
N                 -9.96956257    8.60840634   -0.45933051
N                -10.36316857    6.26756134   -0.46172051
N                 -8.52241857    4.81851234   -0.46188251
N                -10.73452157    3.83718034   -0.46392951
N                 -4.05493857    6.39369034   -0.45719051
C                  7.31394443   -6.17471466   -0.45998151
C                  6.92974143   -8.49938266   -0.46229651
C                  9.11044443   -7.66903966   -0.45961451
C                 11.30956543   -6.82679866   -0.45682751
C                  9.49905743   -5.32088166   -0.45713051
C                  7.64802743   -3.84402466   -0.45745751
H                 12.37632343   -7.00320966   -0.45600451
N                  8.22502943   -8.71724866   -0.46129851
N                  6.47293243   -7.22075766   -0.46165551
N                  8.64555743   -6.38598566   -0.45891151
N                 10.44996043   -7.88481666   -0.45856651
N                 10.83659643   -5.54831466   -0.45609551
N                  9.00038443   -4.04510666   -0.45638851
N                  6.82846143   -4.92348266   -0.45927651
N                  7.07905643   -2.50866766   -0.45666851
N                  5.98086443   -9.59861166   -0.46386951
H                -11.76734657    3.99811534   -0.46424451
H                  6.35570343  -10.57449766   -0.46420551
H                  7.73333243   -1.69355966   -0.45550851
H                 -6.96241557  -15.70423966   -0.45812051
H                -13.98819057  -10.10700066   -0.45727351
H                 -3.39265757    7.20242534   -0.45613051
O                 -1.38429737   -2.01584783    0.17510509
O                 -0.23842455   -2.12435764    0.17656370

这是我用ms计算的  用的castep板块 2×2×1晶胞




delete 一个氮原子的结果和输入分别如下

用H饱和了C的成键




两者计算的结果如下

高斯结果

ms结果

两者计算的共同点：在 氮空位的氮化碳表面的氧气键长都变的更长了  
但是我感觉ms这个板块  原始氮化碳的相对能量是不是变化太大了   至今没找到原因  请大家指教

sobereva

如置顶的新社员必读贴、论坛首页的公告栏、版头的红色大字非常明确所示，求助帖必须在帖子标题明确体现出此帖内容是求助或提问，并清楚、准确反映出帖子具体内容，避免有任何歧义，仔细看http://bbs.keinsci.com/thread-9348-1-1.html。我已把你的不恰当标题“ms和高斯计算氮化碳氧气的吸附能”改了，以后务必注意，下次将删帖+扣分处理。

当前帖子和“量化理论”毫无直接联系，仔细看清楚本版块版头的分类说明，什么样的话题才适合选择“量化理论”，别随便选择“量化理论”分类。这次给你改了，以后注意，下次乱设分类将扣分处理。

太长的输入文件别直接贴，应当当做附件上传，以后注意。