| 查看: 4937 | 回复: 8 | ||
alick_cxj新虫 (小有名气)
|
[求助]
高斯考虑溶剂效应后优化出错
|
|
高斯刚刚开始学不久,基本上都不怎么懂,故向各位请教。 我在之前不考虑溶剂效应时是可以优化成功的,但考虑后就一直出错,不知道怎么回事,请各位帮忙看看应该怎么解决,谢谢了 错误: Polarizable Continuum Model (PCM) ================================= Model : C-PCM. Atomic radii : UFF (Universal Force Field). Polarization charges : Total charges. Charge compensation : None. Solution method : Matrix inversion. Cavity : GePol (RMin=0.200 OFac=0.890). Default sphere list used, NSphG= 52. Tesserae with average area of 0.200 Ang**2. 1st derivatives : Analytical q*S(x)*q algorithm for C-PCM (CSMder). Cavity 1st derivative terms included. Solvent : Ethanol, Eps = 24.550000 Eps(inf)= 1.847000 RSolv = 2.180000 Ang. ------------------------------------------------------------------------------ BldSpC: Error generating genealogic tree for sphere 365 at level 15 Error termination via Lnk1e in /home/xjchen/g03/l301.exe at Mon Nov 7 18:59:31 2011. Job cpu time: 5 days 17 hours 29 minutes 31.5 seconds. File lengths (MBytes): RWF= 1318 Int= 0 D2E= 0 Chk= 38 Scr= 1 输入文件: %chk=P1_cpcm.chk %mem=500MW %nproc=8 # opt freq pbe1pbe/6-31+g(d) scrf=(cpcm,solvent=ethanol,read) iop(1/8=5) P1_cpcm 0 1 坐标…… Radii=UFF |
» 猜你喜欢
研究发现一个新的临界系数,独立于临界压缩因子
已经有7人回复
夏威夷大学生物胶体与界面实验室诚招全奖博士生和博士后
已经有1人回复
物理化学论文润色/翻译怎么收费?
已经有248人回复
200年来人们一直无法破解的难题: 气液临界点状态方程
已经有7人回复
关于引力延迟下双星系统角动量不守恒的疑惑---求高手解答!
已经有10人回复
求助立方晶体In(OH)3的CIF卡片
已经有2人回复
七嗪类物质合成求助
已经有0人回复
[电子教材]天津大学物理化学(第七版)上、下册
已经有1人回复
阴离子交换膜电解二氧化碳还原
已经有0人回复
» 本主题相关价值贴推荐,对您同样有帮助:
请问Gauss中溶剂化的计算怎样指定溶剂
已经有7人回复
高斯计算单点能出错
已经有7人回复
请各位大侠看下高斯出错
已经有3人回复
高斯03中PCM溶剂模型 优化
已经有9人回复
高斯09的溶剂化能到底包含哪几项?
已经有6人回复
【求助】计算G09荧光例子时报错
已经有13人回复
【求助】高斯计算td的问题
已经有9人回复
【整理】今天整理硬盘上的资料,偶然发现的好玩的小东西就直接发上来吧
已经有105人回复
【求助】咨询:高斯优化中局域能量最小值和全局最小值的问题
已经有6人回复
【求助】高斯09优化激发态几何后算频率出错。求教改正方案。
已经有4人回复
【求助】高斯优化过渡态出错
已经有14人回复
【求助】加溶剂与不加溶剂,高斯计算的机时上的差异??
已经有8人回复
【活动】征集:高斯09能做但高斯03做不了的研究实例
已经有8人回复
alick_cxj
新虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 902
- 帖子: 52
- 在线: 16.6小时
- 虫号: 1307761
- 注册: 2011-05-27
- 性别: GG
- 专业: 理论和计算化学
|
额,还是没人回复啊 我把Gaussian给我的回复贴出来,大家看看吧,高斯公司回复还挺快的,昨天下午问,今天早上就收到邮件了。 Thank you for your inquiry. This is a geometrical problem regarding the generation of the solute cavity by the PCM models. Most of the difficulties with PCM in G03 were of this nature. The generation of the solute cavity in the G03 PCM implementation was made out of a series of interlocking spheres centered in the atoms of the molecule and determined by the chosen radii. In addition to this, the cavity generator included a number of additional spheres to smooth the surface, which are important around the edges where the atom-centered spheres interlock. The added spheres however can add too much complexity and problems associated with very small tesserae, which can make geometry optimizations difficult to converge. The most problematic cases for geometry optimizations in G03 are cage-like structures, weakly bound complexes consisting of several molecules, organometallic compounds whose ligands are close to each other, etc. The G03-style cavity tried to fill in some of the empty space with added spheres. Also, on some symmetric systems, the generation of the solute cavity may be problematic. Note that G09 features a new formalism for computing the reaction field that ensures continuity, smoothness and robustness, including continuous derivatives with respect to atomic positions and external perturbing fields. This formalism, initially proposed in 1999 by Karplus and York [York99], has been developed and generalized for whole the PCM family of solvation methods [Scalmani09]. This is achieved by expanding the apparent surface charge that builds up at the solute-solvent interface ("Polarization charges" in terms of spherical Gaussian functions located at each surface element in which the cavity surface is discretized. Discontinuities in the surface derivatives are removed by effectively smoothing the regions where the spheres intersect, which also eliminates the need for the problematic "added spheres" used in the default G03 solute cavities to smooth the cavity surfaces. For new studies, we recommend the use of PCM in the new G09 implementation in order to avoid this type of problems. If you still use G03, you may try to reduce the number of added spheres. The addition of spheres to smooth the surface cavity in PCM models is controlled by the "RMin" and "OFac" PCM options (additional PCM keywords to be included at the end of your input file). These parameters are the minimum radius for any added sphere and the maximum overlap between spheres, respectively. The default values for these parameters are "RMin=0.2" and "OFac=0.89", so, for instance, #p RMin=0.5 OFac=0.8 will increase the minimum size of the additional spheres (RMin > default), and reduce the number of these spheres (OFac < default) resulting in a smaller total number of added spheres than in the default case. By modifying the parameters this way with respect to the default values (that is by reducing the number of added spheres), the solute cavity will become a bit rougher but it will not change the overall shape of the cavity. The default OFac and RMin values ("RMin=0.2" and "OFac=0.89" are generally appropriate for most cases. The values shown above ("RMin=0.5" and "OFac=0.8" generally work well when problems with the cavity appear. Any intermediate values between the default values and the ones given here may also be acceptable. Keep in mind that the modification of the "OFac" and "RMin" parameters will make changes in the total energy that may be noticeable. Thus, the safest choice will be to use the same "OFac" and "RMin" values for the generation of the cavities in all the PCM calculations for this same study, so the cavities in each case are comparable, and so will be the energies. Sincerely, Fernando R. Clemente, Ph.D. Technical Support Gaussian, Inc. help@gaussian.com |
2楼2011-11-09 10:22:32
3楼2011-11-09 13:17:43
xiongb123456
铜虫 (小有名气)
- 应助: 17 (小学生)
- 金币: 207.8
- 散金: 134
- 红花: 5
- 帖子: 183
- 在线: 112.7小时
- 虫号: 896637
- 注册: 2009-11-07
- 专业: 电化学
【答案】应助回帖
★ ★ ★
alick_cxj(金币+1): 没人回答就给你了吧 2011-11-10 09:33:36
zhou2009:编辑内容 2011-11-14 08:47
zhou2009(金币+3): 2011-11-14 08:49:23
alick_cxj(金币+1): 没人回答就给你了吧 2011-11-10 09:33:36
zhou2009:编辑内容 2011-11-14 08:47
zhou2009(金币+3): 2011-11-14 08:49:23
|
高斯公司的回复 是让你 改参数嘛!?“RMin=0.5 OFac=0.8”不要用默认值了 你要不听听高斯的 自己写个RMin=0.5 OFac=0.8在后面 试试吧 我最近算溶剂效应用的是PCM 不是CPCM 我也是新手 原来不会的还可以问高斯公司 学习了 [ Last edited by zhou2009 on 2011-11-14 at 08:47 ] |
4楼2011-11-09 16:43:28
5楼2011-11-09 19:12:31
alick_cxj
新虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 902
- 帖子: 52
- 在线: 16.6小时
- 虫号: 1307761
- 注册: 2011-05-27
- 性别: GG
- 专业: 理论和计算化学
6楼2011-11-09 22:43:22
alick_cxj
新虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 902
- 帖子: 52
- 在线: 16.6小时
- 虫号: 1307761
- 注册: 2011-05-27
- 性别: GG
- 专业: 理论和计算化学
7楼2011-11-09 22:44:46
xiongb123456
铜虫 (小有名气)
- 应助: 17 (小学生)
- 金币: 207.8
- 散金: 134
- 红花: 5
- 帖子: 183
- 在线: 112.7小时
- 虫号: 896637
- 注册: 2009-11-07
- 专业: 电化学
8楼2011-11-14 08:14:32
alick_cxj
新虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 902
- 帖子: 52
- 在线: 16.6小时
- 虫号: 1307761
- 注册: 2011-05-27
- 性别: GG
- 专业: 理论和计算化学
9楼2011-11-14 10:21:19












回复此楼
in terms of spherical Gaussian functions located at each surface element in which the cavity surface is discretized. Discontinuities in the surface derivatives are removed by effectively smoothing the regions where the spheres intersect, which also eliminates the need for the problematic "added spheres" used in the default G03 solute cavities to smooth the cavity surfaces.