计算化学公社

 找回密码 Forget password
 注册 Register
楼主 Author: sobereva
打印 Print 上一主题 Last thread 下一主题 Next thread

[综合交流] 简单量子化学问题答疑专帖(旧版,已关帖)

   关闭 [复制链接 Copy URL]

1665

帖子

5

威望

4788

eV
积分
6553

Level 6 (一方通行)

喵星人

11416#
发表于 Post on 2020-7-27 00:48:08 | 只看该作者 Only view this author
ljk 发表于 2020-7-27 00:45
装虚拟机再装linux吗?

直接搜WSL2,按照微软官方的指点装,完了选个发行版,一般是ubuntu,然后一切正常
用的是win的文件系统,/home在C盘某个文件夹里,无所谓不用管,正常的C盘D盘在/mnt/c/和/mnt/d/这样的目录里面,会用linux的话很快就能上手了,比一般的虚拟机效率高不少。

147

帖子

4

威望

1577

eV
积分
1804

Level 5 (御坂)

科研泥瓦匠

11417#
发表于 Post on 2020-7-27 01:11:49 | 只看该作者 Only view this author
老师,您好,附件是一个分子的两个构象(1与2),1与2构象互变(8-10键旋转)的过程是否有过渡态?
1,http://bbs.keinsci.com/thread-2388-3-1.html阅读该贴,通过gentor旋转8-10键180度,产生的结构逐一进行优化,其中能量最高结构作为过渡态初猜。
2,这种没有旧键断裂新键生成的过渡态文件怎么输入,一般的过渡态关键词opt=(calcfc,ts,noeigen)等是否适用?是否需要添加其他关键词吗?

TIM图片20200727005948.png (174.78 KB, 下载次数 Times of downloads: 33)

TIM图片20200727005948.png

p101.xyz

2.62 KB, 下载次数 Times of downloads: 0

r203.xyz

2.62 KB, 下载次数 Times of downloads: 0

终日寻春不见春
芒鞋踏破岭头云
归来偶把梅花嗅
春在枝头已十分

1665

帖子

5

威望

4788

eV
积分
6553

Level 6 (一方通行)

喵星人

11418#
发表于 Post on 2020-7-27 02:59:22 | 只看该作者 Only view this author
欢乐多 发表于 2020-7-27 01:11
老师,您好,附件是一个分子的两个构象(1与2),1与2构象互变(8-10键旋转)的过程是否有过渡态?
1,htt ...

就用扫描得到的最高能量的结构作初猜结构,正常优化过渡态即可

6万

帖子

99

威望

6万

eV
积分
125148

管理员

公社社长

11419#
 楼主 Author| 发表于 Post on 2020-7-27 04:36:32 | 只看该作者 Only view this author
ljk 发表于 2020-7-27 00:03
我的电脑16GB 6核12线程
INPUT:
%mem=200MW

“Windows版本gaussian09只能设置200mw内存”这是严重的误解。要分清楚32bit和64bit,认真看此文
Gaussian的安装方法及运行时的相关问题
http://sobereva.com/439http://bbs.keinsci.com/thread-10814-1-1.html
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

6万

帖子

99

威望

6万

eV
积分
125148

管理员

公社社长

11420#
 楼主 Author| 发表于 Post on 2020-7-27 04:39:26 | 只看该作者 Only view this author
xiaobchine 发表于 2020-7-26 23:06
这是输入文件:
%nprocshared=16
%chk=I-anion.chk

一个原子根本没有可以被优化的变量,显然会报错,这是最基本常识

此问题你分成两步做更是完全莫名其妙

阴离子算能量必须带弥散函数,因此你当前的计算根本没有实际意义,仔细看
谈谈弥散函数和“月份”基组
http://sobereva.com/119
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

6万

帖子

99

威望

6万

eV
积分
125148

管理员

公社社长

11421#
 楼主 Author| 发表于 Post on 2020-7-27 04:40:45 | 只看该作者 Only view this author
dingniu2 发表于 2020-7-26 23:51
请问老师,平时做高斯计算,可以一次写很多个gjf文件,将文件名写到一个bat文件里,例如:
g16 1.gjf
g16 ...

当然得cd进入相应目录下再运行xtb

参考此文脚本里的folder_runall.sh
使用Gaussian时的几个实用脚本和命令
http://sobereva.com/258http://bbs.keinsci.com/thread-190-1-1.html
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

6万

帖子

99

威望

6万

eV
积分
125148

管理员

公社社长

11422#
 楼主 Author| 发表于 Post on 2020-7-27 04:42:34 | 只看该作者 Only view this author
欢乐多 发表于 2020-7-27 01:11
老师,您好,附件是一个分子的两个构象(1与2),1与2构象互变(8-10键旋转)的过程是否有过渡态?
1,htt ...

拿个便宜的级别比如PM7或PM6D3,做个二面角旋转的柔性扫描,其中如果出现了极大点,尝试当过渡态初猜,照常搜索过渡态。半经验这种便宜级别找到过渡态之后,再用其结构当初始结构,用DFT进一步优化得到更精确的TS
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

15

帖子

0

威望

31

eV
积分
46

Level 2 能力者

11423#
发表于 Post on 2020-7-27 09:08:10 | 只看该作者 Only view this author
sobereva 发表于 2020-7-27 04:36
“Windows版本gaussian09只能设置200mw内存”这是严重的误解。要分清楚32bit和64bit,认真看此文
Gaussi ...

我按照链接所说,设置Default.Rou文件,它本是空白的,填写资源配置后,就报错,

559

帖子

0

威望

6210

eV
积分
6769

Level 6 (一方通行)

11424#
发表于 Post on 2020-7-27 10:05:19 | 只看该作者 Only view this author
sobereva 发表于 2020-7-27 04:40
当然得cd进入相应目录下再运行xtb

参考此文脚本里的folder_runall.sh

谢谢老师,平时用nohup xtb 1.xyz --chrg 0 --uhf 0 --gfn 2 --ohess verytight > 1.out 2>&1 &提交单个xtb任务,我将批量提交的脚本改成:
#!/bin/bash
shopt -s nullglob
for i in `ls -F |grep /`
do
cd $i
for inf in *.xyz
do
echo Running ${inf} ...
time xtb < ${inf} > ${inf//xyz/out} --chrg 0 --uhf 0 --gfn 2 --ohess verytight
echo ${inf} has finished
echo
done
cd ..
done
还是不能运行,请问老师是不是哪里改错了,谢谢老师!!!

6万

帖子

99

威望

6万

eV
积分
125148

管理员

公社社长

11425#
 楼主 Author| 发表于 Post on 2020-7-28 05:37:13 | 只看该作者 Only view this author
ljk 发表于 2020-7-27 09:08
我按照链接所说,设置Default.Rou文件,它本是空白的,填写资源配置后,就报错,

明显是你填的内容不当
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

6万

帖子

99

威望

6万

eV
积分
125148

管理员

公社社长

11426#
 楼主 Author| 发表于 Post on 2020-7-28 05:38:41 | 只看该作者 Only view this author
dingniu2 发表于 2020-7-27 10:05
谢谢老师,平时用nohup xtb 1.xyz --chrg 0 --uhf 0 --gfn 2 --ohess verytight > 1.out 2>&1 &提交单个x ...

我对远程提交的事不清楚。先确保脚本在单机上能正常运行
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

82

帖子

4

威望

2540

eV
积分
2702

Level 5 (御坂)

11427#
发表于 Post on 2020-7-28 06:15:11 | 只看该作者 Only view this author
dingniu2 发表于 2020-7-27 10:05
谢谢老师,平时用nohup xtb 1.xyz --chrg 0 --uhf 0 --gfn 2 --ohess verytight > 1.out 2>&1 &提交单个x ...

脚本的计算命令跟平时用的命令顺序是一致的。
把这一行time xtb < ${inf} > ${inf//xyz/out} --chrg 0 --uhf 0 --gfn 2 --ohess verytight
改成:
time xtb ${inf} --chrg 0 --uhf 0 --gfn 2 --ohess verytight > ${inf%.*}.out

另外for 循环之前可以考虑加环境变量如:
export OMP_NUM_THREADS=4
export MKL_NUM_THREADS=${OMP_NUM_THREADS}
ulimit -s unlimited
环境变量可以根据自己资源情况修改

559

帖子

0

威望

6210

eV
积分
6769

Level 6 (一方通行)

11428#
发表于 Post on 2020-7-28 11:46:50 | 只看该作者 Only view this author
sobereva 发表于 2020-7-28 05:38
我对远程提交的事不清楚。先确保脚本在单机上能正常运行

谢谢老师!

559

帖子

0

威望

6210

eV
积分
6769

Level 6 (一方通行)

11429#
发表于 Post on 2020-7-28 11:48:15 | 只看该作者 Only view this author
winnerwill 发表于 2020-7-28 06:15
脚本的计算命令跟平时用的命令顺序是一致的。
把这一行time xtb < ${inf} > ${inf//xyz/out} --chrg 0 - ...

非常感谢,真的可以了!!!谢谢!!!

273

帖子

0

威望

4430

eV
积分
4703

Level 6 (一方通行)

11430#
发表于 Post on 2020-7-28 17:04:39 | 只看该作者 Only view this author
Sob老师,对于有机小分子来说,由于np态跃迁禁阻,那么荧光/磷光的发射态是否一定是pp态呢?如果S1是np态,S2是pp态,那么荧光一定是从pp态往下发射吗?

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-23 02:47 , Processed in 0.246547 second(s), 22 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list