计算化学公社

 找回密码 Forget password
 注册 Register
Views: 11495|回复 Reply: 3
打印 Print 上一主题 Last thread 下一主题 Next thread

[新手求助] 色散矫正相关问题

[复制链接 Copy URL]

46

帖子

0

威望

405

eV
积分
451

Level 3 能力者

                        
一. 请问各位老师,我最近在对一个二聚体中间的两个具有面面堆积的分子进行势能面扫描,在两个中间具有面面堆积的环上引入虚原子,借助dimerscan 和 XYZ2QC生成高斯输入文件,做单点计算绘制势能面扫描曲线。我在优化这个二聚体的激发态S1时,采用与文献相同泛函 # m062x/6-31g(d,p) td opt 我优化这个二聚体有将近六七天时间,都没有得到结果。我这里应该考虑加入色散矫正项吗,改为# m062xD3/6-31g(d,p) td opt   任务7天后停掉了,请问老师该如何修改,谢谢老师。
二.采用QM/MM的方法,计算中间和最右边两个分子,优化S0和S1态利用四点法计算重组能。中间的分子单晶不含水,最右边分子单晶含水,水分子是普通水,用Mercury建模后,选择中间分子作为QM层,用0放开分子,-1冻结周围分子。我的问题是考虑分子中存在弱相互作用氢键 ,输入文件为# opt oniom(m062x/def2tzvp:uff=qeq)=embed 这里要不要在加上色散矫正项,改为m062xD3 ? 我这里基组选择有点大了优化激发态很慢,改为 Def-TZVP 或者 Def2-SVP 耗时会更少些吗?
谢谢老师

6万

帖子

99

威望

6万

eV
积分
125152

管理员

公社社长

2#
发表于 Post on 2019-5-9 14:13:29 | 只看该作者 Only view this author
“都没有得到结果”、“任务7天后停掉了”这种描述没法给回答者提供任何有效信息。自行看输出文件,把对于解决问题可能有效的信息呈献给回答者

这和加不加D3校正没有直接关系

优化用def-tzvp或6-311G**就足够了,def2-TZVP是给算能量用的

加不加D3仔细看此文
谈谈“计算时是否需要加DFT-D3色散校正?”
http://sobereva.com/413http://bbs.keinsci.com/thread-9772-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

46

帖子

0

威望

405

eV
积分
451

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2019-5-9 15:38:36 | 只看该作者 Only view this author
老师这个是输出文件,我看了下好像是我计算时间超出要求的时间了,您能帮我看一下具体是什么原因吗?谢谢老师了
Sender: LSF System <lsfadmin@node31>
Subject: Job 51442: <OPT> in cluster <njtech2> Exited

Job <OPT> was submitted from host <login01> by user <hlma> in cluster <njtech2>.
Job was executed on host(s) <14*node31>, in queue <high>, as user <hlma> in cluster <njtech2>.
</gpfs/users_home/hlma> was used as the home directory.
</gpfs/users_home/hlma/xueyan/saomiao/s1/.> was used as the working directory.
Started at Results reported on
Your job looked like:

------------------------------------------------------------
# LSBATCH: User input
#!/bin/sh

#BSUB -q high
#BSUB -n 14
#BSUB -R "span[hosts=1]"
#BSUB -cwd .
#BSUB -e %J.err
#BSUB -o %J.out
#BSUB -J OPT

jobname=78s

username=`whoami`
Current_WORKDIR=`pwd`

Inputfile=${jobname}.com

export g09root=/gpfs/soft/g16
#source $g09root/bsd/g09.profile

export GAUSS_EXEDIR=/gpfs/soft/g16
export LD_LIBRARY_PATH=$LD_LIBRARY_PATH:/gpfs/soft/g16

GAUSS_RUNDIR=/tmp/$username/$LSB_JOBID
if [ ! -a $GAUSS_RUNDIR ]; then
   echo "Scratch directory $GAUSS_RUNDIR created."
   mkdir -p $GAUSS_RUNDIR
fi

cp $Current_WORKDIR/$Inputfile  $GAUSS_RUNDIR

cd $GAUSS_RUNDIR

GAUSS_SCRDIR=$GAUSS_RUNDIR/${jobname}
if [ ! -a $GAUSS_SCRDIR ]; then
   echo "Scratch directory $GAUSS_SCRDIR created."
   mkdir -p $GAUSS_SCRDIR
fi
export GAUSS_SCRDIR
echo "Using $GAUSS_SCRDIR for temporary Gaussian 09 files."

echo "Starting Gaussian run at" `date`

time /gpfs/soft/g16/g16 < $GAUSS_RUNDIR/${jobname}.com >$GAUSS_RUNDIR/${jobname}.log

echo "Finished Gaussian run at" `date`

echo $Current_WORKDIR

mv $GAUSS_RUNDIR/*.chk $Current_WORKDIR

(... more ...)
------------------------------------------------------------

TERM_RUNLIMIT: job killed after reaching LSF run time limit.
Exited with exit code 140.

Resource usage summary:

    CPU time :                                   8428928.00 sec.
    Max Memory :                                 30546 MB
    Average Memory :                             13074.28 MB
    Total Requested Memory :                     -
    Delta Memory :                               -
    Max Swap :                                   128 MB
    Max Processes :                              5
    Max Threads :                                19
    Run time :                                   604808 sec.
    Turnaround time :                            638611 sec.

The output (if any) follows:

Scratch directory /tmp/hlma/51442 created.
Scratch directory /tmp/hlma/51442/78s created.
Using /tmp/hlma/51442/78s for temporary Gaussian 09 files.
Starting Gaussian run at Wed May 1 04:42:35 CST 2019


PS:

Read file <51442.err> for stderr output of this job.

46

帖子

0

威望

405

eV
积分
451

Level 3 能力者

4#
 楼主 Author| 发表于 Post on 2019-5-9 15:41:59 | 只看该作者 Only view this author
这个二聚体我优化过T1态用UDFT 出来了,S1态的总也优化不出来卡了半个月了

本版积分规则 Credits rule

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

GMT+8, 2026-2-23 15:23 , Processed in 0.166657 second(s), 23 queries , Gzip On.

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