| 查看: 1758 | 回复: 7 | ||
[求助]
求大神们帮忙看看这个程序哪儿需要改啊 已有1人参与
|
|
初学molpro，学着计算NO基态的势能，想要单点能扫描，这下面是计算程序，但是出来的结果不对啊，能不能请大神们帮忙看看，万分感谢！ ***,NO potential geometry={h n; o,n,r(i); } basis=vtz distances=[0.4,0.41,0.42,0.43,0.44,0.45,0.46,0.47,0.48,0.49,0.5,0.51,0.52,0.53,0.54,0.55,0.56,0.57,0.58,0.59,0.6,0.61,0.62,0.63,0.64,0.65,0.66,0.67,0.68,0.69,0.7,0.71,0.72,0.73,0.74,0.75,0.76,0.77,0.78,0.79,0.8,0.81,0.82,0.83,0.84,0.85,0.86,0.87,0.88,0.89,0.9,0.91,0.92,0.93,0.94,0.95,0.96,0.97,0.98,0.99,1.0,1.01,1.02,1.03,1.04,1.05,1.06,1.07,1.08,1.09,1.1,1.11,1.12,1.13,1.14,1.15,1.16,1.17,1.18,1.19,1.2,1.21,1.22,1.23,1.24,1.25,1.26,1.27,1.28,1.29,1.3,1.31,1.32,1.33,1.34,1.35,1.36,1.37,1.38,1.39,1.4,1.41,1.42,1.43,1.44,1.45,1.46,1.47,1.48,1.49,1.5,1.51,1.52,1.53,1.54,1.55,1.56,1.57,1.58,1.59,1.6,1.61,1.62,1.63,1.64,1.65,1.66,1.67,1.68,1.69,1.7,1.71,1.72,1.73,1.74,1.75,1.76,1.77,1.78,1.79,1.8,1.81,1.82,1.83,1.84,1.85,1.86,1.87,1.88,1.89,1.9,1.91,1.92,1.93,1.94,1.95,1.96,1.97,1.98,1.99,2.0,2.01,2.02,2.03,2.04,2.05,2.06,2.07,2.08,2.09,2.1,2.2,2.3,2.4,2.5,2.6,2.7,2.8,2.9,3.0,3.1,3.2,3.3,3.4,3.5,3.6] i=0 do ir=1,#distances i=i+1 r(i)=distances(ir) {hf; wf,15,1,1 occ,4,2,2 } escf(i)=energy enddo {table,r,escf title,results for NO } --- |
回复此楼» 猜你喜欢
华南师范大学广州市光谱分析与功能探针重点实验室诚招青年英才和博士后岗位研究人员
已经有10人回复
-
200年来人们一直无法破解的难题： 气液临界点状态方程
已经有7人回复
-
物理化学论文润色/翻译怎么收费?
已经有181人回复
-
关于引力延迟下双星系统角动量不守恒的疑惑---求高手解答！
已经有10人回复
-
求标准卡片Cu0.4In0.4Zn1.2S2 PDF#32-0340
已经有2人回复
-
量子点电致发光器件（QLED）：蓝光量子点（QDs）和纳米氧化锌电子传输材料
已经有0人回复
-
需要一个圆偏振光诱导化学反应，请问使用圆偏振光作为光源，应该需要哪些设备搭建
已经有0人回复
-
求ENDNOTE破解软件安装包
已经有4人回复
-
河海大学 （211.双一流高校）钙钛矿光电实验室丁勇团队欢迎有志青年！
已经有3人回复
-
黑龙江省原子能研究院化学课题组招收2026级化学专业或化学相关专业硕士生
已经有0人回复
-
26博士申请
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 图形求助，大神帮忙~~~
已经有19人回复
- 求大神帮忙翻译
已经有2人回复
- 哪位大神能帮我把这个图画到word里哈？不是用截图的形式哦~然后用word回复给我！
已经有5人回复
- 急急急急急急，TEM测试的图太模糊了，求大神帮忙分析原因！
已经有19人回复
- 请大神帮忙看看这是什么植物？？？？？
已经有6人回复
- 新手求助~~求大神帮忙！！
已经有3人回复
- 积分求助 求大神帮忙！
已经有7人回复
- 哪位大神帮帮忙
已经有4人回复
- 很简单的Fortran赋值出错啦，求大神们帮忙找bug...
已经有11人回复
- 请各位大神帮忙看看！！！急急
已经有7人回复
- 求大神帮忙解答！
已经有8人回复
- 请问有没有哪位环境类的大神做分子模拟之类的
已经有5人回复
- 热重-差热分析的图，哪个大神帮忙解释一下呗……
已经有6人回复
- 修改图标问题，请各位大神帮忙！
已经有5人回复
- 问题解决，谢谢关注！
已经有20人回复
- 第一次投英文的，作者指导有点没看明白，求大神帮忙，谢谢
已经有17人回复
- 求两个值中最大的值但编写的代码不知道哪里错了，求大神帮忙。
已经有10人回复
- 求大神帮忙，怎么用tecplot出这种图
已经有6人回复
- 求大神帮我改一下下面的程序！matlab迭代法的 不知道错在哪里！
已经有7人回复
luqing6879
木虫 (著名写手)
- 应助: 8 (幼儿园)
- 金币: 3566
- 散金: 249
- 红花: 3
- 帖子: 1908
- 在线: 808.5小时
- 虫号: 614833
- 注册: 2008-09-28
- 专业: 理论和计算化学
2楼2014-05-07 01:01:49
3楼2014-05-07 08:17:26
luqing6879
木虫 (著名写手)
- 应助: 8 (幼儿园)
- 金币: 3566
- 散金: 249
- 红花: 3
- 帖子: 1908
- 在线: 808.5小时
- 虫号: 614833
- 注册: 2008-09-28
- 专业: 理论和计算化学
4楼2014-05-07 08:44:42
5楼2014-05-07 08:48:30
【答案】应助回帖
|
你这里有两个问题。 1、Molpro中，内坐标的键长单位默认是Bohr，你小于2Bohr的算了那么多都没用的，只要看2到5Bohr的就差不多了。 2、你说的没有一个阱的结构，是因为你的态算错了，你算的这个态不是成键的态，所以键长越长，能量越低。 解决办法有两个， 一是不使用对称性，算出来的就是最低的态。在geom上一行加上nosym，然后occ后面的参数改成8就行了。需要注意的是，这种方法不一定能得到正确的结果，需要旋转轨道尝试很多次。 二是把态给选对。N-O的最大的abel群是c2v群，有4个不可约表示，需要搞清楚每个不可约表示的原子轨道占据数，以及波函数属于哪个不可约表示。 我用wf,15,1,1; occ,5,2,1,0或者wf,15,2,1; occ,5,2,1,0算出来都比你原来的要低，你的occ,4,2,2肯定是不对的。 [ 发自手机版 http://muchong.com/3g ] |
6楼2014-05-09 21:49:48
【答案】应助回帖
★ ★ ★ ★ ★
kdswx: 金币+5, ★★★★★最佳答案 2014-05-09 22:34:40
kdswx: 金币+5, ★★★★★最佳答案 2014-05-09 22:34:40
|
1.5 -128.6871215 1.6 -128.9296432 1.7 -129.0904337 1.8 -129.1926184 1.9 -129.2528221 2.0 -129.2830743 2.1 -129.2921377 2.2 -129.2864358 2.3 -129.2707005 2.4 -129.2484223 2.5 -129.2221507 2.6 -129.1937002 2.7 -129.1643240 2.8 -129.1348963 2.9 -129.1061222 3.0 -129.0788913 3.1 -129.0547177 3.2 -129.0348454 3.3 -129.0191149 3.4 -129.0067277 3.5 -128.9968953 3.6 -128.9889930 3.8 -128.9772466 4.0 -128.9690866 4.2 -128.9632158 4.4 -128.9588936 4.6 -128.9556628 4.8 -128.9532197 5.0 -128.9513529 [ 发自手机版 http://muchong.com/3g ] |
7楼2014-05-09 21:55:59
8楼2014-05-09 22:34:06
