记录学习用ORCA结合Multiwfn计算RESP、RESP2脚本计算原子电荷

主要是记录一下踩坑的过程,只是整理过程,其中具体参数不做过多解释

1、复合物PDB准备

受体和配体提取:可选择的软件有很多,maestro、pymol、DS、MOE这些软件都可以

注意要把配体也就是小分子加氢保存为mol2,比如test.mol2,受体蛋白需要保留一条链且去除溶剂和金属离子,保存为PBD,test.pdb格式

2、配体拓扑文件准备并计算RESP2电荷

使用sobtop生成配体的拓扑文件,然后使用相应的量子化学计算软件,(ORCA或高斯)与Multifwn联用计算RESP2电荷。

这里使用了sobtop+ORCA+Multifwn+linux的方法,可以参考分子动力学模拟详细过程记录 - 哔哩哔哩

sobtop安装及使用

下载和安装:Sobtop

使用下载压缩包后,直接解压缩,在windows系统直接双击sobtop.exe使用

打开后出现如下界面:

a、把tes.mol2文件拖进去,或者输入文件路径,回车

b、选择2回车,输入生成目录:E:\XX\test.gro,意思是生成gro文件,结束后回到主页

c、输入1回车,生成拓扑文件,选择分子力场2或其他力场,然后回车,选择4回车

d、后面会让输入路径也就是生成目录:E:\XX\test.top,回车

e、再输入路径:E:\XX\test.itp回车

f、直接叉掉或者输入q,回车返回主页

现在已经存在了下面的文件,

ORCA安装

ORCA安装分为windows和linux,为了计算方便l,我选择的是linux版本

ORCA下载:https://orcaforum.kofo.mpg.de/,需要注册才能登陆下载

安装方法下面这篇博客说的很清楚

量子化学程序ORCA的安装方法 - 思想家公社的门口:量子化学·分子模拟·二次元 (sobereva.com)

需要注意的点:一定要先看看自己电脑有没有openmpi,以及把路径加入环境变量,要退出后,重新打开终端设置才能生效

Multiwfn安装

Multiwfn下载地址:Multiwfn,尽量选择高的版本,而且还区分了普通版和no_GUI版本,区别在于Multiwfn在Linux下安装的中文说明 - 思想家公社的门口:量子化学·分子模拟·二次元 (sobereva.com)

里面说的很清楚如果没有安装motif图形库提供的libXm.so.4文件,普通版无法启动,且motif库安装需要root权限,Multiwfn在官网上特意提供了noGUI版,由于它在启动时不需要motif库,因此也就不需要弄到root权限来安装它,故普通用户都可以很容易地安装。noGUI版的缺点在于没法使用任何Multiwfn与图形有关的功能

计算RESP2电荷:

看下面这篇就够了,主要是要把RESP2_ORCA.sh文件里的路径写对

http://sobereva.com/637

注意可能在计算过程中会涉及权限问题:需要赋权,chmod +x  <路径>

把RESP2_ORCA.sh 和test.mol2放在一文件夹内运行

注意:在使用没有图形界面的版本时,RESP2_ORCA.sh文件中Multiwfn替换成Multiwfn_noGUI才能正常运行

后面得到带.chg后缀的文件就是得到的电荷数据

### Abaqus 中装配体的拓扑优化方法 #### 创建模型并定义材料属性 在进行装配体的拓扑优化之前,需先建立精确的几何模型,并赋予合适的材料特性。对于复杂的装配件而言,建议采用精细化建模策略来捕捉所有重要的细节特征。 #### 设置边界条件与载荷工况 合理施加约束以及外部作用力是确保计算收敛性的前提之一。针对不同类型的连接方式(如螺栓联接、焊接等),应考虑其特有的力学行为模式[^1]。 #### 定义接触相互作用 当处理由多个部件组成的系统时,正确描述各组件间的相对运动至关重要。这通常涉及到指定恰当形式的摩擦系数及其变化规律;同时也要注意区分法向和切向方向上的响应差异[^4]。 #### 构造初始网格布局 为了提高求解效率,在不影响最终结果精度的基础上尽可能简化离散化方案不失为一种明智之举。特别是面对那些具有不规则外形轮廓的对象来说,“自底向上”的思路往往能够带来意想不到的效果——即优先填充内部区域然后再向外扩展直至覆盖整个实体表面为止[^3]。 #### 启动作业流程并配置参数选项 进入Job Manager界面之后新建一个名为Optimization的任务实例,随后按照提示逐步完成各项设定工作直到准备就绪可以提交给后台引擎去执行实际运算任务为止。期间务必仔细检查每一个输入框内的数值是否准确无误以免造成不必要的麻烦。 ```python from abaqus import * from abaqusConstants import * # 假设已经完成了前面几步的操作... mdb.models['Model-1'].optimizationTasks.create( name='Topology Optimization', region=mdb.models['Model-1'].parts['Part-1'].sets['All Elements'], densityLowerBound=0.0, densityUpperBound=1.0, objectiveFunctionType=MINIMIZE_COMPLIANCE, numberOfIterations=200) session.writeToReplay('submit_job.py') ```
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值