| 查看: 607 | 回复: 2 | |||
huilaoshu999金虫 (著名写手)
|
[交流]
【求助】请教关于溶质与溶剂分子间的静电作用能的问题? 已有2人参与
|
|
------------------------------------------------------------------------------ Polarizable Continuum Model (PCM) ================================= Model : PCM. Atomic radii : UA0 (Simple United Atom Topological Model). Polarization charges : Total charges. Charge compensation : None. Solution method : Matrix inversion. Cavity : GePol (RMin=0.200 OFac=0.890). Default sphere list used, NSphG= 3. Tesserae with average area of 0.200 Ang**2. Solvent : Acetonitrile, Eps= 36.640000. ------------------------------------------------------------------------------ Using symmetry in molecular cavity generation. One-electron integrals computed using PRISM. NBasis= 63 RedAO= T NBF= 46 17 NBsUse= 63 1.00D-06 NBFU= 46 17 Harris functional with IExCor= 205 diagonalized for initial guess. ExpMin= 4.38D-02 ExpMax= 4.17D+03 ExpMxC= 6.27D+02 IAcc=2 IRadAn= 4 AccDes= 0.00D+00 HarFok: IExCor= 205 AccDes= 0.00D+00 IRadAn= 4 IDoV=1 ScaDFX= 1.000000 1.000000 1.000000 1.000000 Initial guess orbital symmetries: Occupied (A1) (A1) (A1) (A1) (A1) (A1) (E) (E) (A1) (E) (E) Virtual (E) (E) (A1) (A1) (E) (E) (A1) (E) (E) (A1) (E) (E) (A1) (E) (E) (A1) (A1) (A1) (A1) (E) (E) (E) (E) (A1) (A1) (E) (E) (A1) (E) (E) (A1) (A1) (E) (E) (A1) (E) (E) (A1) (E) (E) (E) (E) (A1) (E) (E) (E) (E) (A1) (A1) (A1) (A1) (A1) The electronic state of the initial guess is 1-A1. Requested convergence on RMS density matrix=1.00D-04 within 128 cycles. Requested convergence on MAX density matrix=1.00D-02. Requested convergence on energy=1.00D-02. No special actions if energy rises. Keep R1 integrals in memory in canonical form, NReq= 2767344. Error on total polarization charges = 0.00935 SCF Done: E(RHF) = -131.940813624 A.U. after 5 cycles Convg = 0.9612D-04 -V/T = 2.0018 S**2 = 0.0000 -------------------------------------------------------------------- Variational PCM results ======================= Total free energy in solution: with all non electrostatic terms (a.u.) = -131.932196 -------------------------------------------------------------------- (Polarized solute)-Solvent (kcal/mol) = -7.40 -------------------------------------------------------------------- Cavitation energy (kcal/mol) = 9.36 Dispersion energy (kcal/mol) = -4.20 Repulsion energy (kcal/mol) = 0.25 Total non electrostatic (kcal/mol) = 5.41 -------------------------------------------------------------------- Partition over spheres: Sphere on Atom Surface Charge GEl GCav GDR 1 C1 2.32 0.004 -0.02 0.67 -0.04 2 C2 64.67 -0.263 -2.53 6.01 -2.56 3 N6 23.88 0.242 -4.71 2.68 -1.35 Added spheres: 7.51 0.008 -0.14 0.00 0.00 -------------------------------------------------------------------- ********************************************************************** Population analysis using the SCF density. ********************************************************************** Orbital symmetries: Occupied (A1) (A1) (A1) (A1) (A1) (A1) (E) (E) (A1) (E) (E) Virtual (A1) (A1) (E) (E) (E) (E) (A1) (E) (E) (A1) (A1) (E) (E) (A1) (E) (E) (A1) (A1) (A1) (E) (E) (E) (E) (A1) (A1) (A1) (E) (E) (E) (E) (A1) (A1) (A1) (E) (E) (E) (E) (A1) (E) (E) (E) (E) (A1) (E) (E) (E) (E) (A1) (A1) (A1) (A1) (A1) The electronic state is 1-A1. Alpha occ. eigenvalues -- -15.58977 -11.30022 -11.27522 -1.24896 -1.03740 Alpha occ. eigenvalues -- -0.69516 -0.62497 -0.62497 -0.57522 -0.48287 Alpha occ. eigenvalues -- -0.48287 Alpha virt. eigenvalues -- 0.06679 0.09087 0.11881 0.11881 0.14621 Alpha virt. eigenvalues -- 0.14621 0.15710 0.18266 0.18266 0.20642 Alpha virt. eigenvalues -- 0.29575 0.29850 0.29850 0.30087 0.36654 Alpha virt. eigenvalues -- 0.36654 0.40985 0.47406 0.69593 0.85120 Alpha virt. eigenvalues -- 0.85120 0.91827 0.91827 0.99891 1.06998 Alpha virt. eigenvalues -- 1.13497 1.16367 1.16367 1.24274 1.24274 Alpha virt. eigenvalues -- 1.36343 1.43002 1.65732 1.68865 1.68865 Alpha virt. eigenvalues -- 1.76572 1.76572 1.97844 2.20130 2.20130 Alpha virt. eigenvalues -- 2.26970 2.26970 2.51998 2.58696 2.58696 Alpha virt. eigenvalues -- 3.11613 3.11613 3.27249 3.67329 4.52141 Alpha virt. eigenvalues -- 4.59286 5.05432 Condensed to atoms (all electrons): 1 2 3 4 5 6 1 C 4.787052 0.138084 -0.050464 -0.050464 -0.050464 0.754286 2 C 0.138084 5.396979 0.403581 0.403581 0.403581 -0.052259 3 H -0.050464 0.403581 0.422474 -0.020543 -0.020543 -0.002631 4 H -0.050464 0.403581 -0.020543 0.422474 -0.020543 -0.002631 5 H -0.050464 0.403581 -0.020543 -0.020543 0.422474 -0.002631 6 N 0.754286 -0.052259 -0.002631 -0.002631 -0.002631 6.888662 Mulliken atomic charges: 1 1 C 0.471969 2 C -0.693548 3 H 0.268125 4 H 0.268125 5 H 0.268125 6 N -0.582796 Sum of Mulliken charges= 0.00000 Atomic charges with hydrogens summed into heavy atoms: 1 1 C 0.471969 2 C 0.110828 3 H 0.000000 4 H 0.000000 5 H 0.000000 6 N -0.582796 Sum of Mulliken charges= 0.00000 Electronic spatial extent (au): Charge= 0.0000 electrons Dipole moment (field-independent basis, Debye): X= 0.0000 Y= 0.0000 Z= -5.1515 Tot= 5.1515 Quadrupole moment (field-independent basis, Debye-Ang): XX= -17.7371 YY= -17.7371 ZZ= -21.9228 XY= 0.0000 XZ= 0.0000 YZ= 0.0000 Traceless Quadrupole moment (field-independent basis, Debye-Ang): XX= 1.3952 YY= 1.3952 ZZ= -2.7905 XY= 0.0000 XZ= 0.0000 YZ= 0.0000 Octapole moment (field-independent basis, Debye-Ang**2): XXX= 0.0000 YYY= 0.8857 ZZZ= -17.6082 XYY= 0.0000 XXY= -0.8857 XXZ= -2.8911 XZZ= 0.0000 YZZ= 0.0000 YYZ= -2.8911 XYZ= 0.0000 Hexadecapole moment (field-independent basis, Debye-Ang**3): XXXX= -23.9495 YYYY= -23.9495 ZZZZ= -184.9841 XXXY= 0.0000 XXXZ= 0.0000 YYYX= 0.0000 YYYZ= -1.4541 ZZZX= 0.0000 ZZZY= 0.0000 XXYY= -7.9832 XXZZ= -31.1307 YYZZ= -31.1307 XXYZ= 1.4541 YYXZ= 0.0000 ZZXY= 0.0000 N-N= 5.886190824948D+01 E-N=-4.255290960902D+02 KE= 1.317065684090D+02 Symmetry A' KE= 1.269810328033D+02 Symmetry A" KE= 4.725535605678D+00 这个gaussian03计算结果文件中,静电作用力从哪里看?我只看到了空穴能、散射-排斥作用能和非静电作用能总和啊,溶剂化能无法计算啊 ! |
» 猜你喜欢
书籍求助:汽车市场营销理论与实务(电子版)——章小平
已经有0人回复
东方理-中科大联合博士生招聘
已经有0人回复
物理化学论文润色/翻译怎么收费?
已经有234人回复
半透性容器进行高湿或者低湿研究
已经有0人回复
钴酸锂半电池小倍率容量上不去
已经有1人回复
吉林大学材料物理本科生求问调剂信息
已经有23人回复
2026第二届光电子与半导体器件前沿技术研讨会——光电子赋能·半导体创芯!
已经有1人回复
光电子赋能·半导体创芯!
已经有1人回复
光电子赋能·半导体创芯!
已经有0人回复
QE计算电声耦合的时候报错Error in routine lambda (100)wrong or too many modes
已经有2人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐,对您同样有帮助:
考虑溶剂效应下,如何计算两分子间的作用能
已经有3人回复
请教一些静电势的问题
已经有3人回复
请教关于溶剂和引发剂选择的问题
已经有5人回复
请教一个静电场,恒定电场的问题
已经有17人回复
【求助】静电自组装中溶解聚电解质PDDA和PSS的溶剂!
已经有11人回复
【求助】请教:空白溶剂相关问题?
已经有10人回复
【请教】胶体纳米粒子的在氨基表面的静电吸附问题
已经有22人回复
【求助】有机溶剂之间溶解机理及溶质存在形式
已经有3人回复
【求助】溶剂与溶质分子之间的作用
已经有7人回复
hzfish
金虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 1515.8
- 散金: 10
- 红花: 3
- 帖子: 604
- 在线: 358.1小时
- 虫号: 753222
- 注册: 2009-04-20
- 性别: GG
- 专业: 理论和计算化学
2楼2010-04-22 22:33:07
zhangmt
至尊木虫 (著名写手)
我叫MT
- QC强帖: 5
- 应助: 99 (初中生)
- 金币: 6961.8
- 散金: 10406
- 红花: 49
- 帖子: 1761
- 在线: 763.2小时
- 虫号: 880392
- 注册: 2009-10-22
- 性别: GG
- 专业: 理论和计算化学
★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
gavinliu7390(金币+2):谢谢交流! 2010-04-23 22:11
小木虫(金币+0.5):给个红包,谢谢回帖交流
gavinliu7390(金币+2):谢谢交流! 2010-04-23 22:11
|
空穴能、散射、排斥作用能三项加起来不就是那个非静电作用项嘛。 (Polarized solute)-Solvent那个就是静电项啊。 两个加起来就是△G(sol)啊。 也就是 Total free energy in solution with all non electrostatic terms (溶液中的能量) 这一项 减去 [ Last edited by zhangmt on 2010-4-22 at 23:28 ] |

3楼2010-04-22 23:26:52












回复此楼