求助用orca模拟分子内键断裂

lingchen

使用约束（如固定距离或施加拉伸力）来模拟键断裂过程。这个代码应该怎么写啊，查手册也没查到

wxyhgk

1. Constraining bond distances        :   { B N1 N2 value C }
   
2. Constraining bond angles           :   { A N1 N2 N1 value C }
   
3. Constraining dihedral angles       :   { D N1 N2 N3 N4 value C }
   
4. Constraining cartesian coordinates :   { C N1 C }
   
5. 
   
6. 
   
7. 约束键长：   { B N1 N2 value C }
   
8. # 其中 N1、N2 是原子编号，value 是约束的键长值
   
9. 
   
10. 约束键角：   { A N1 N2 N3 value C }
    
11. # 其中 N1、N2、N3 是原子编号，value 是约束的键角值（单位：度）
    
12. 
    
13. 约束二面角： { D N1 N2 N3 N4 value C }
    
14. # 其中 N1、N2、N3、N4 是原子编号，value 是约束的二面角值（单位：度）
    
15. 
    
16. 约束笛卡尔坐标： { C N1 C }
    
17. # 其中 N1 是原子编号，约束该原子的 x、y、z 坐标

复制代码

lingchen

wxyhgk 发表于 2024-9-28 13:42

谢谢谢谢，我去手册里面再找一下

lingchen

wxyhgk 发表于 2024-9-28 13:42

我还有个问题，如果我想让他一次计算多种约束条件下的状态该怎么写代码啊。比如我想知道某个键长分别是1.2，1.25，1.3，1.35，1.4，1.45，1.5条件下的能量就只能写多个输入文件一个一个计算吗，不能在一个输入文件里一起计算吗

wxyhgk

lingchen 发表于 2024-9-28 20:25
我还有个问题，如果我想让他一次计算多种约束条件下的状态该怎么写代码啊。比如我想知道某个键长分别是1. ...

直接扫描不就行了吗？

1. %geom Scan
   
2. # n1,n2换成你需要扫描的原子
   
3.         B n1 n2 = 1.2, 1.5, 7 # 这将扫描从1.2到1.5的7个点
   
4. 
   
5.         end
   
6.       end

复制代码

sobereva

利用Gaussian或ORCA程序把分子结构拉直的几种方法
http://sobereva.com/594（http://bbs.keinsci.com/thread-22945-1-1.html）

了解势能面扫描常识，ORCA同理
详谈使用Gaussian做势能面扫描
http://sobereva.com/474（http://bbs.keinsci.com/thread-12660-1-1.html）