NAMD-xtb QM/MM interface

CasDonkin · 发表于 Post on 2024-4-16 10:15:20

付老师您好，我最近在使用您的NAMD/xtb interface在QM(GFN2-xtb)/MM(Charmm36m)级别下计算酶催化过程，根据前人报道，其中一步由一个亲核进攻和一个质子转移协同发生，请问colvars是否能正确地在QM/MM中处理2个CV，描绘二维反应自由能面吗？在使用NAMD作为QM/MM的interface时，1维的PMF比较常见，在此请教付老师，提前向您表示感谢

fhh2626 · 发表于 Post on 2024-4-16 14:24:14

CasDonkin 发表于 2024-4-16 10:15
付老师您好，我最近在使用您的NAMD/xtb interface在QM(GFN2-xtb)/MM(Charmm36m)级别下计算酶催化过程，根据 ...

可以的。Colvars是描述几何/集合变量的方法，和使用QM还是MM本身不直接相关

CasDonkin · 发表于 Post on 2024-7-25 18:17:04

付老师您好，我最近打算尝试使用您的Interface在QM(GFN2-xtb)/MM(Charmm36m)级别下计算活性位点含过渡金属的酶催化过程，类似于下文：https://doi.org/10.1021/acs.jpcb.1c10200。

我在阅读了您的脚本后，发现您的脚本中并未通过设置--uhf 选项加入单电子数的设置。根据我的理解，我是否可以在脚本中run_qmmm()函数部分的subprocess.call([XTBDIR, xtbxyz, '--grad', '--chrg', str(QMCHARGE[qmpart]), '--gfn', str(GFNVER)], stdout = outputRedirect, stderr = outputRedirect) 中加入'---uhf ' 选项即可运行计算，而无需其他改动？（如对六重态Mn2+设置为uhf = 5，对于单重态和三重态的Ni2+分别设置uhf为0或2）

提前向您表示感谢

fhh2626 · 发表于 Post on 2024-7-28 18:06:25

CasDonkin 发表于 2024-7-25 18:17
付老师您好，我最近打算尝试使用您的Interface在QM(GFN2-xtb)/MM(Charmm36m)级别下计算活性位点含过渡金属 ...

可以的

po390 · 发表于 Post on 2024-9-21 17:15:12

付老师,您好,您的namd-xtb脚本非常好用帮我解决了不少问题,非常感谢!
不过使用中我遇到一个问题,发现无论xtb都是单核在跑(虽然对我的单核的速度也够用了),namd如何可以多核去调用xtb?
NAMD2.14和3.0都试过,都是单核去调用xtb,无法多核调用(设置多核环境变量都无效果).
其实,之前在用namd提供的mopac interface同样存在类似的问题.
这个是否是NAMD本身的限制?

fhh2626 · 发表于 Post on 2024-9-21 21:12:15

po390 发表于 2024-9-21 17:15
付老师,您好,您的namd-xtb脚本非常好用帮我解决了不少问题,非常感谢!
不过使用中我遇到一个问题,发现无论x ...

你是不是看错了，默认xtb使用所有的核心跑的(top可以看到xtb线程占用超过100%的CPU)，要设置
export OMP_NUM_THREADS=<ncores>,1
环境变量才能更改核心数。

po390 · 发表于 Post on 2024-9-22 10:52:06

fhh2626 发表于 2024-9-21 21:12
你是不是看错了，默认xtb使用所有的核心跑的(top可以看到xtb线程占用超过100%的CPU)，要设置
export OMP ...

谢谢老师的提醒!
之前已经设置过,但还是无法多核运行xtb,
不过老师的话提醒了我,我发现是我个人使用习惯的问题,我在运行namd时加上+setcpuaffinity 导致运行时xtb无法多核,去掉之后就可以了,
不过为啥+setcpuaffinity的影响会这里?

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[程序/脚本开发] NAMD-xtb QM/MM interface