24小时热门版块排行榜    

查看: 2451  |  回复: 16

daiyulan85

木虫 (正式写手)

[求助] 这种激发态扫描图是怎么做的(Linear interpolation step不知道具体是指什么)

如题,看起来应该是横坐标是扫描的步数,但是这个Linear interpolation step具体是指什么呢?google查到了Linear interpolation 的意思(Linear interpolation is a method of curve fitting using linear polynomials),仅仅是设定某一个参数进行激发态扫描吗?求解。

Energies of the S0,S1,S2,S3 excited states along the LIIC connecting the FCpoint to the S1(2APB) for 5’-2APT-3'
回复此楼
计算化学光谱
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

070241078

金虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
daiyulan85(金币+10): ★★★很有帮助 2012-03-07 09:17:23
应该就是在初始值的基础上加上去的步数,我的势能面曲线图横坐标是以一定的步长间隔的二面角的度数
2楼2012-03-07 09:10:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

daiyulan85

木虫 (正式写手)

引用回帖:
2楼: Originally posted by 070241078 at 2012-03-07 09:10:08:
应该就是在初始值的基础上加上去的步数,我的势能面曲线图横坐标是以一定的步长间隔的二面角的度数

你的是激发态扫描吗?你说的那种情况也可以,就是不知道是不是图中的那一种。感觉他是在两个点之间之间取值,并且在这个两个点的直线上。
计算化学光谱
3楼2012-03-07 09:20:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

070241078

金虫 (小有名气)

【答案】应助回帖

★ ★
御剑江湖(金币+2): 呵呵,谢谢 2012-03-08 19:26:23
是啊,也是激发态势能曲线扫描。应该都是一样的吧,不管横坐标怎么样,表示的意义都差不多。你要是用上面这种图,可能需要注明你的步长值?我也是初学者,所以也就了解到这种程度啦!
4楼2012-03-07 09:27:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pwzhou

铁杆木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
gmy1990(金币+3): 2012-03-07 11:26:28
daiyulan85(金币+10): ★★★很有帮助 2012-03-07 17:00:36
这个是在第一个构型和最后一个构型的所有内坐标(3N-6个)之间插值扫描出来的。
比如第一个构型的某个键长为1.2,而最后一个构型的该键长为1.3,那么中间9个点的构型该键长为1.21,1.22,1.23,1.24,1.25,1.26,1.27,1.28,1.29。依次类推,对所有的3N-6个坐标都做一样的处理。

这个实现起来也很简单,只需要知道第一个结构和最后一个结构的Z-Matrix,然后对Z-Matrix中的3N-6个值进行差值即可。
几点注意事项:第一个结构和最后一个结构的Z-Matrix中的原子序号和内坐标必须完全一致;键长和键角可以之间用N_start+(N_end-N_start)/n*i 的方式来处理,但是二面角不可以这么处理,需要根据差值的大小做一点特殊处理,否则插值出来的结构就会完全乱套,几个原子会完全挤在一起。
5楼2012-03-07 10:22:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pwzhou

铁杆木虫 (正式写手)

引用回帖:
2楼: Originally posted by 070241078 at 2012-03-07 09:10:08:
应该就是在初始值的基础上加上去的步数,我的势能面曲线图横坐标是以一定的步长间隔的二面角的度数

这个只是对某一个参数的扫描,而LIIC是对所有内坐标做差值扫描出来的。
6楼2012-03-07 10:26:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

viger87

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★
感谢参与,应助指数 +1
御剑江湖(金币+4): 呵呵,谢谢 2012-03-08 19:26:36
All intermediate geometries were produced by the equation Q i =Qi (r) + n .{Qi (p) -Q i (r)}, where Qi (r) is any coordinate of the reactant, and Qi(p) is the same coordinate of the product – in this case the conical intersection. The interpolation parameter e is changed in the interval from 0 (for reactant) to 1 (for product – conical intersection S0–S1 ) with a step of 0.1.
也就是把产物的坐标减去反应物的坐标,然后n等分这个值,再一份一份的把这个等分的值加上反应物的坐标,就是要的每个状态的坐标了
有花堪折直须折,末待花落空折枝!
7楼2012-03-07 18:59:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

daiyulan85

木虫 (正式写手)

引用回帖:
5楼: Originally posted by pwzhou at 2012-03-07 10:22:05:
这个是在第一个构型和最后一个构型的所有内坐标(3N-6个)之间插值扫描出来的。
比如第一个构型的某个键长为1.2,而最后一个构型的该键长为1.3,那么中间9个点的构型该键长为1.21,1.22,1.23,1.24,1.25,1.26,1.27 ...

那请问在扫描的时候是怎么处理的呢?是每一个值都要扫一遍(比如说键长从1.2扫到1.27,每一步为0.01),还是处理完所有的值,得到所有构型再扫描?这么多构型怎么操作呢?
计算化学光谱
8楼2012-03-08 08:43:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

daiyulan85

木虫 (正式写手)

引用回帖:
7楼: Originally posted by viger87 at 2012-03-07 18:59:47:
All intermediate geometries were produced by the equation Q i =Qi (r) + n .{Qi (p) -Q i (r)}, where Qi (r) is any coordinate of the reactant, and Qi(p) is the same coordinate of the product – in t ...

不知道这个conical intersection 的坐标是怎么得到的,我看原文有一句这样的话:The LIIC between initial FC point and S1(2APB) are shown in Figure 3,这个FC point 就是S0的结构吗?
计算化学光谱
9楼2012-03-08 08:52:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

070241078

金虫 (小有名气)

引用回帖:
6楼: Originally posted by pwzhou at 2012-03-07 10:26:04:
这个只是对某一个参数的扫描,而LIIC是对所有内坐标做差值扫描出来的。

谢谢你的意见啊
10楼2012-03-08 09:19:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 daiyulan85 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 6GojgJvkDudM 2026-08-16 4/200 2026-08-17 13:49 by EYK67XJMj64E
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +33 医学老男孩 2026-08-13 78/3900 2026-08-17 13:22 by ziyangfang
[基金申请] filecode=后面第一个是大写字母 +7 wangze12014 2026-08-14 8/400 2026-08-17 12:45 by a441914426
[基金申请] 今天系统多次维护,明天很可能放榜! +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 i7NFEVbjQMM5 2026-08-16 4/200 2026-08-17 12:10 by 4GBAYCdVQoK3
[考研] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 6GojgJvkDudM 2026-08-16 4/200 2026-08-17 11:46 by 4GBAYCdVQoK3
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 6GojgJvkDudM 2026-08-16 3/150 2026-08-17 09:49 by xLVPIRuSvCUe
[论文投稿] 岩土工程学报什么时候才能终审完呐 25+3 yeager111 2026-08-10 4/200 2026-08-16 20:42 by tfang
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 6/300 2026-08-16 20:07 by 奶牛小黑
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +3 Tide man 2026-08-14 3/150 2026-08-16 17:47 by jurkat.1640
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 26/1300 2026-08-16 12:35 by 启萌科技
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +10 majunge000 2026-08-11 12/600 2026-08-16 08:18 by xli1984
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 是这周出结果还是下周出结果? +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[基金申请] filecode +15 documentary 2026-08-10 17/850 2026-08-14 10:08 by kissu88
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by 且听虎啸
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
信息提示
请填处理意见