计算化学公社

标题: ORCA6.1版本怎么才能计算二极化率或更高阶？ [打印本页]

作者
Author: wenbo6558 时间: 2025-9-27 23:52
标题: ORCA6.1版本怎么才能计算二极化率或更高阶？
各位老师好，

我用HyperPol 1 能计算出一阶极化率，HyperPol 2也没得到二极化率的结果，只得到了下面结果的数据。

我的计算代码是下面这段，我怎么样才能让orca计算二阶极化率呢？

跟polar 2/3被废除有关系吗？我用文档Compound的实例试一下也没有二阶极化率的。

输入文件
!freq wb97x-d3 6-311+G(d,p) RIJCOSX AutoAux DefGrid3 PrintBasis PrintMOs PrintGap
! KeepInts

%pal nprocs 6 end

%MaxCore 4500

%scf
  Convergence          # The default convergence is between medium and strong
            Tight
  IntName "orca_lv-1" #Read the existing integrals from a previous calculation in conventional mode (use together with %scf IntName "PreviousBaseName" end)
  MaxIter=500
end

%moinp "orca_lv-1.gbw" # 复用之前的波函数

%elprop
Polar ture    # 计算偶极矩μ和极化率α (基础)
HyperPol ture # 一阶超极化率
   # er阶超极化率
PrintLevel 2
maxiter 500  # 增加CPHF迭代次数

PolarVelocity true # polarizability w.r.t. velocity perturbations
PolarDipQuad  true # dipole-quadrupole polarizability
PolarQuadQuad true # quadrupole-quadrupole polarizability

#原子偶极矩和极化率
PolarAtom true
Dipole true
DipoleAtom true
Quadrupole true
QuadrupoleAtom true
end

*xyz 0 1 #带电数量自旋多重度

结果：
Hyperpolarizability calculation:

---------------------------------------------------
      STATIC HYPERPOLARIZABILITY TENSOR
---------------------------------------------------

Method          : SCF
Type of density : Electron Density
Type of derivative : Electric Field (Direction=X)
Multiplicity    : 1
Irrep             : 0
Basis             : AO

  The raw Cartesian tensor (atomic units):

   ( x x x ):       -124.237138
   ( x x y ):       270.830135
   ( x x z ):       66.861209
   ( x y x ):       270.830138
   ( x y y ):       -462.956555
   ( x y z ):       84.733171
   ( x z x ):       66.861206
   ( x z y ):       84.733170
   ( x z z ):       -20.110770
   ( y x x ):       270.830138
   ( y x y ):       -462.956555
   ( y x z ):       84.733171
   ( y y x ):       -462.956554
   ( y y y ):       78.963468
   ( y y z ):       -241.588169
   ( y z x ):       84.733170
   ( y z y ):       -241.588172
   ( y z z ):       84.903195
   ( z x x ):       66.861206
   ( z x y ):       84.733170
   ( z x z ):       -20.110770
   ( z y x ):       84.733170
   ( z y y ):       -241.588172
   ( z y z ):       84.903195
   ( z z x ):       -20.110771
   ( z z y ):       84.903192
   ( z z z ):       -121.774743

Hyperpolarizability calculation done in 1193.0 sec

-----------
SCF HESSIAN
-----------

作者
Author: sobereva 时间: 2025-9-28 00:25
不要自己在标题里手写【求助】这种碍眼的标签，http://bbs.keinsci.com/thread-9348-1-1.html里明确说了。这次给你改了，以后注意

作者
Author: sobereva 时间: 2025-9-28 00:26
没法直接给出第二超极化率，需要自己基于解析的第一超极化率做有限差分算

作者
Author: wenbo6558 时间: 2025-9-28 11:03

sobereva 发表于 2025-9-28 00:25
不要自己在标题里手写【求助】这种碍眼的标签，http://bbs.keinsci.com/thread-9348-1-1.html里明确说了。 ...

好的，谢谢老师

作者
Author: wenbo6558 时间: 2025-9-28 13:38
本帖最后由 wenbo6558 于 2025-10-1 10:34 编辑

sobereva 发表于 2025-9-28 00:26
没法直接给出第二超极化率，需要自己基于解析的第一超极化率做有限差分算

老师，有限差用Multiwfn可以计算吗？

我把.out转换成。molden后，

24 (Hyper)polarizability analysis

就出现了这些

0 Return to main menu
1 Parse output file of (hyper)polarizability task of Gaussian and calculate various related quantities
2 Study (hyper)polarizability by sum-over-states (SOS) method and perform two/three-level analysis
3 (hyper)polarizability density analysis
4 Calculate atomic polarizabilities in molecule
5 Visualize (hyper)polarizability via unit sphere and vector representations

不知道我这个对不对。没有出现有限差的选项

作者
Author: sobereva 时间: 2025-9-28 14:16

wenbo6558 发表于 2025-9-28 13:38
老师，有限差用Multiwfn可以计算吗？

我把.out转换成。molden后，

跟Multiwfn没任何关系
自己按照第二超极化率和第一超极化率之间的关系，基于不同电场下的第一超极化率用计算器手算第二超极化率

作者
Author: wenbo6558 时间: 2025-9-28 14:38

sobereva 发表于 2025-9-28 14:16
跟Multiwfn没任何关系
自己按照第二超极化率和第一超极化率之间的关系，基于不同电场下的第一超极化率用 ...

好的，谢谢老师

作者
Author: wenbo6558 时间: 2025-9-30 20:38
本帖最后由 wenbo6558 于 2025-10-1 10:34 编辑

sobereva 发表于 2025-9-28 14:16
跟Multiwfn没任何关系
自己按照第二超极化率和第一超极化率之间的关系，基于不同电场下的第一超极化率用 ...

老师好，我用官网提供的numericalPolarizability.cmp计算各方向的场的第一极化率

但是调用过程中发现，只能调用一个CPU单元，无法并行计算，这个可以怎么添加？直接加inp里面调用不了。

另外，如何才能把各个场的计算一起进行，不用一个一个算？

谢谢老师指教

%Compound "numericalPolarizability.cmp"
with
molecule = "orca_lv-1.xyz";
method    = "wb97x-d3";
basis    = "6-311+G(d,p)";
restOfInput = "VeryTightSCF";
E_Field    = 0.0 0.0 -0.0005;

end

作者
Author: sobereva 时间: 2025-10-1 05:32

wenbo6558 发表于 2025-9-30 20:38
老师好，我用官网提供的numericalPolarizability.cmp计算各方向的场的第一极化率

发帖时别用那么大字号，看着很难受。如果嫌默认的文字小，你自己的浏览器里增加放大比例就完了

总共只需要+x -x +y -y +z -z算6次，没必要非得用什么组合脚本
怎么并行下文说了
量子化学程序ORCA的安装方法
http://sobereva.com/451（http://bbs.keinsci.com/thread-11697-1-1.html）

非要同时算所有电场的，同时提交六个任务就完了

顺带一提，6-311+G(d,p)算超极化率是极其糟糕的，该用什么看
谈谈量子化学中基组的选择
http://sobereva.com/336（http://bbs.keinsci.com/thread-3545-1-1.html）

作者
Author: wenbo6558 时间: 2025-10-1 10:45

sobereva 发表于 2025-10-1 05:32
发帖时别用那么大字号，看着很难受。如果嫌默认的文字小，你自己的浏览器里增加放大比例就完了

总共只 ...

好的，谢谢老师。字号问题已更改

作者
Author: wenbo6558 时间: 2025-10-1 22:22
本帖最后由 wenbo6558 于 2025-10-1 22:27 编辑

sobereva 发表于 2025-10-1 05:32
发帖时别用那么大字号，看着很难受。如果嫌默认的文字小，你自己的浏览器里增加放大比例就完了

总共只 ...

老师好，

并行的基本文件我都有，结构优化的时候也能并行运算，没有问题。

但是用inp文件调用numericalPolarizability.cmp后，在restOfInput或blocksInput中加入 %pal nprocs 6 end或者！pal6，程序都还是单CPU计算，没有并行计算。

这种情况，我应该怎么写才能让他启动并行计算？

下面是我的inp文件
%Compound "numericalPolarizability.cmp"
with
molecule = "orca_lv-1.xyz";
method    = "wb97x-d3";
basis    = "6-311+G(d,p)";
restOfInput = "VeryTightSCF";
blocksInput = "";
E_Field    = 0.0 0.0 -0.0005;

end

谢谢老师指点！

作者
Author: xiaobanma99 时间: 2025-10-9 15:50
本帖最后由 xiaobanma99 于 2025-10-10 07:58 编辑

学习了

欢迎光临计算化学公社 (http://bbs.keinsci.com/)