计算化学公社

 找回密码 Forget password
 注册 Register
Views: 7317|回复 Reply: 2
打印 Print 上一主题 Last thread 下一主题 Next thread

[其它程序] PYXAID NAMD 模拟求助

[复制链接 Copy URL]

68

帖子

0

威望

639

eV
积分
707

Level 4 (黑子)

跳转到指定楼层 Go to specific reply
楼主
各位老师好。请问在Pyxaid中进行NAMD模拟的时候,需要设置basis state 、active space和initial condition(不清楚应该怎么翻译合适,见谅)。其中对于所谓的initial condition怎么理解呢?主要我看到了两种不同的表达方式
   1、I 、J 分别定义研究的基态和激发态,此为官网的教程
# Excitations from 102 (HOMO) to [103,... Nmax]
Nmin = 90
HOMO = 102
LUMO = HOMO+1
Nmax = 122

# Initial condition excitation: I->J (HOMO -> LUMO+5)
I = 102
J = 108
# This is how I,J are mapped to index of the corresponding basis state
# see module lazy.py for more details. This is not the most general way
# of such mapping though.
ex_indx = 1 + (J-LUMO)*(HOMO+1 - Nmin) + (I-Nmin) = 78

# Each entry of the list below is an initial condition. It is also a list
# but containing only 2 elements - first is the time step at which we start
# the dynamics, the second is the index of the excited state configuration
params["iconds"] = [ [0,ex_indx], [25, ex_indx], [49, ex_indx] ]



此处第一列是否代表开始将电子从HOMO激发到LUMO+5轨道的时间,还是说开始进行NAMD模拟的时间点呢?


2、第二种表达方式
# Initial conditions
nmicrost = len(params["states"]) # number of (micro)states
ic = []
i = 0
while i < 10:
    j = 0
    while j < nmicrost:
        ic.append([50*i,j])
        j = j + 1
    i = i + 1
params["iconds"] = ic
在这种情况下,似乎有更多的initial condition,而且i、 j的取值感觉不清楚(j=0,是否代表LUMO,然后类推),所以也更不清楚这两种表达方式的区别和联系。

最后希望熟悉Pyxaid程序的老师、同学不吝赐教,十分感谢!!!

23

帖子

0

威望

289

eV
积分
312

Level 3 能力者

2#
发表于 Post on 2020-8-14 09:20:25 | 只看该作者 Only view this author
您好,教程可以发给我一份吗,没有找到官网的教程;QQ邮箱:1059656898@qq.com非常感谢

68

帖子

0

威望

639

eV
积分
707

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2020-8-22 12:40:30 | 只看该作者 Only view this author
Regina 发表于 2020-8-14 09:20
您好,教程可以发给我一份吗,没有找到官网的教程;QQ邮箱:非常感谢

google 出来的第一个不就是吗

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2025-8-18 02:37 , Processed in 0.148297 second(s), 25 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list