计算化学公社

标题: 计算小白从哪里学习计算优化结构和激发态的具体操作流程 [打印本页]

作者
Author: 身在盘丝洞 时间: 2023-7-20 16:00
标题: 计算小白从哪里学习计算优化结构和激发态的具体操作流程
本人需要用激发态S1的NTO分析来解释实验现象，但博文Gaussian中用TDDFT计算激发态和吸收、荧光、磷光光谱的方法[url=http://sobereva.com/314（http://bbs.keinsci.com/thread-2413-1-1.html）]http://sobereva.com/314（http://bbs.keinsci.com/thread-2413-1-1.html）[/url]中的操作流程不熟悉（之前用的菜鸟用的gaussview），除了高斯使用说明外有什么途径可以了解到具体详细的操作步骤

作者
Author: wzkchem5 时间: 2023-7-20 16:27
具体列出对于哪几步不熟悉，我们好针对性地给你推荐更进一步的学习材料

作者
Author: 身在盘丝洞 时间: 2023-7-20 17:00

wzkchem5 发表于 2023-7-20 16:27
具体列出对于哪几步不熟悉，我们好针对性地给你推荐更进一步的学习材料

第一步：在(1)得到的结构下，先在溶剂下做基态单点计算，把基态的溶剂信息写到chk里，同时得到基态在平衡溶剂下的能量（SCF Done后面的值）：
# PBE1PBE/6-311G* scrf(read,solvent=ethanol)，末尾空一行写noneq=write
  第二步：做电子激发计算，并读取上一步在chk里写入的基态溶剂信息，由此得到激发态在非平衡溶剂下的能量：
# PBE1PBE/6-311G* TD scrf(externaliteration,read,solvent=ethanol) geom=check guess=read，末尾空一行写noneq=read

在优化结构之后我把基态单点以下面的内容提交计算
%nprocshared=64
%mem=16GB
%chk=bo-cz18.chk
# pbepbe/6-31g(d,p) scrf=(solvent=acetonitrile,read)

[No Title]

0 1
H                2.66915683 -0.59021315 2.79316130
H                3.22945407 0.47130159 4.09292306

noneq=read
然后算完就导出fch文件用gauss view打开再输出一个gjf文件再加入博文的关键词却计算失败，不知道我的操作流程有什么问题吗老师

作者
Author: wzkchem5 时间: 2023-7-20 18:09

身在盘丝洞发表于 2023-7-20 10:00
第一步：在(1)得到的结构下，先在溶剂下做基态单点计算，把基态的溶剂信息写到chk里，同时得到基态在平 ...

但凡说“计算失败”，必须给出详细的报错信息，并上传完整的输出文件。失败的原因又不是唯一的，不可能你一说失败我们就知道该怎么解决

作者
Author: 身在盘丝洞 时间: 2023-7-20 18:44

wzkchem5 发表于 2023-7-20 18:09
但凡说“计算失败”，必须给出详细的报错信息，并上传完整的输出文件。失败的原因又不是唯一的，不可能你 ...

我想知道我的流程对吗老师我想先排除一下这个
opt的fch文件→用gaussview打开加入基态单点导出gjf文件参照博文改关键词在最后加noneq=write→提交计算平台得到转换后的fch文件→继续用gaussview打开改关键词，最后一行加noneq=read→提交计算得到的fch文件用multwfn做NTO分析

作者
Author: sobereva 时间: 2023-7-20 19:03

身在盘丝洞发表于 2023-7-20 18:44
我想知道我的流程对吗老师我想先排除一下这个
opt的fch文件→用gaussview打开加入基态单点导出gjf文件参 ...

正经用Gaussian做计算的人根本不用gview的calculation setup界面创建输入文件

博文里说得够清楚的了，根本就没让你用gview

北京科音基础量子化学培训班（http://www.keinsci.com/workshop/KBQC_content.html）里有非常完整的考虑非平衡溶剂效应的例子，输入输出文件直接提供了。如果之前一点Gaussian使用经验都没有，先从初级班学起、正经入门计算再说：北京科音初级量子化学培训班（http://www.keinsci.com/workshop/KEQC_content.html）

用pbepbe完全莫名其妙，平时谁拿纯泛函算TDDFT。不从头系统、正确学一遍，会犯各种低级错误。

作者
Author: wzkchem5 时间: 2023-7-20 20:04

身在盘丝洞发表于 2023-7-20 11:44
我想知道我的流程对吗老师我想先排除一下这个
opt的fch文件→用gaussview打开加入基态单点导出gjf文件参 ...

我还是那句话，看到报错信息以后可能我们马上就知道你哪里做错了，把输出文件从你们服务器上下载下来再上传到论坛又不麻烦。
如果你只告诉我们你的计算流程，很可能有的细节部分你会漏掉，比如输入文件最后空行数目不对、提交错了计算节点、环境变量设置错误，等等等等。但是这些问题只要我们一看到完整的输出文件，基本立马就能知道是怎么回事。
建议先认真学习一下http://sobereva.com/620，了解怎么正确地提问题

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