计算化学公社

标题: 哪位有Q-chem破解版？ [打印本页]

作者
Author: lao7 时间: 2015-6-29 17:38
标题: 哪位有Q-chem破解版？
最近用Gaussian计算荧光发射谱，那叫一个慢呀！查了一下，高斯只能有限差分来计算。Q-chem能用解析梯度？哪位仁兄有这款软件？谢谢啦！

作者
Author: sobereva 时间: 2015-6-29 17:43
你误解了高斯。
高斯03就支持CIS的解析梯度，09开始支持TDDFT的解析梯度，荧光直接用高斯算就行了，速度比较理想，完全不需要Q-Chem。
而且Q-Chem的license是和机子绑定的，即便有程序也没法用。
Q-Chem这方面的优点也就是能支持到TDDFT的二阶解析导数，但二阶导数对于荧光计算的激发态优化并不是必须的。

作者
Author: lao7 时间: 2015-6-29 18:25

sobereva 发表于 2015-6-29 17:43
你误解了高斯。
高斯03就支持CIS的解析梯度，09开始支持TDDFT的解析梯度，荧光直接用高斯算就行了，速度比 ...

谢谢回复，我展示我的计算分子给你看看，请帮助看看为什么第五步开始就计算一直很慢了，前4步都是1天完成，第五步开始就一周左右（目前只计算到第6步）：
第一步：
%nprocshared=20
%mem=55GB
%chk=2NF1.chk
# Opt Freq b3lyp/6-31G** scrf=(solvent=water,pcm)

2NF1 ground state

2 1
C                -8.35184600 -4.37086300 -2.04627200
N                -9.69511900 -4.09121800 -1.53104700
C                -9.95150800 -2.64930000 -1.52529000
C                -8.94078600 -1.94531700 -0.61971000
N                -7.57146100 -2.24232700 -1.06487500
C                -7.31522600 -3.69885900 -1.15142700
C                -6.54551500 -1.40991700 -0.64025100
C                -6.71510900 -0.45090700 0.37305700
C                -5.68173700 0.42249700 0.74937700
C                -4.41526000 0.34187700 0.12484500
C                -4.23671300 -0.61190600 -0.89787000
C                -5.26542000 -1.43579800 -1.26374900
N                -5.87674000 1.37800100 1.76419000
C                -4.87448900 2.18817000 2.11448700
C                -3.60569300 2.16861800 1.54208900
C                -3.30831300 1.21612700 0.49572800
C                -7.21037600 1.46406200 2.43010700
C                -7.34226500 2.53819200 3.50019900
C                -2.60562700 3.11925500 2.02701500
O                -1.44188800 3.21253500 1.62592900
O                -2.19189700 1.11085400 -0.07938800
F                -5.06302900 -2.28251700 -2.31191700
H                -8.18044900 -5.45160600 -2.03783700
H                -8.21671800 -4.01100500 -3.08183100
H             -10.39192400 -4.56475000 -2.09981800
H             -10.95861300 -2.46294800 -1.14021900
H                -9.88634500 -2.20817400 -2.53585800
H                -9.10629100 -0.86703200 -0.67843000
H                -9.09504300 -2.26704100 0.42329600
H                -7.36910300 -4.13176100 -0.13982600
H                -6.31721900 -3.87538500 -1.54604800
H                -7.66470500 -0.41092800 0.88169900
H                -3.28548200 -0.66665900 -1.41378300
H                -5.08693700 2.89754900 2.89910500
H                -7.94571900 1.63944600 1.64061300
H                -7.40940100 0.48287000 2.86895200
H                -8.36393100 2.48236700 3.88565500
H                -6.66475600 2.37706400 4.34353000
H                -7.19586200 3.54651800 3.10307400
O                -3.05097200 3.94060100 2.99553200
H                -2.32113000 4.53000700 3.25398100
Ca                0.00000700 2.07663300 0.00000000
C                8.35182500 -4.37093100 2.04622200
N                9.69512000 -4.09123100 1.53108100
C                9.95146000 -2.64930400 1.52536400
C                8.94076800 -1.94533900 0.61973700
N                7.57142800 -2.24240200 1.06481900
C                7.31523400 -3.69894500 1.15132900
C                6.54547800 -1.41001400 0.64016400
C                6.71509500 -0.45096500 -0.37310400
C                5.68172400 0.42243700 -0.74943300
C                4.41522800 0.34178300 -0.12494100
C                4.23665400 -0.61204900 0.89772400
C                5.26535700 -1.43594300 1.26360600
N                5.87675500 1.37800900 -1.76417800
C                4.87452400 2.18821900 -2.11443500
C                3.60571600 2.16864300 -1.54206400
C                3.30829600 1.21606100 -0.49579700
C                7.21040900 1.46413000 -2.43004800
C                7.34235800 2.53846200 -3.49992800
C                2.60568200 3.11935600 -2.02691100
O                1.44194300 3.21263600 -1.62582400
O                2.19187700 1.11077800 0.07930900
F                5.06294300 -2.28269500 2.31174400
H                8.18046500 -5.45168000 2.03775500
H                8.21662400 -4.01109800 3.08177900
H                10.39190700 -4.56475000 2.09988400
H                10.95858100 -2.46291200 1.14035600
H                9.88622400 -2.20819800 2.53593600
H                9.10623400 -0.86705000 0.67848700
H                9.09509600 -2.26703900 -0.42326500
H                7.36918200 -4.13182600 0.13972400
H                6.31721100 -3.87551000 1.54589200
H                7.66471800 -0.41093800 -0.88169000
H                3.28540800 -0.66682800 1.41360600
H                5.08700400 2.89766100 -2.89898900
H                7.94575400 1.63931600 -1.64051100
H                7.40938000 0.48301600 -2.86908800
H                8.36401000 2.48263100 -3.88542100
H                6.66480900 2.37755500 -4.34326900
H                7.19605000 3.54672100 -3.10259800
O                3.05105500 3.94077000 -2.99535700
H                2.32123200 4.53021900 -3.25375900

第二步：

%oldchk=2NFCa1
%chk=2NFCa2
%nprocshared=20
%mem=55GB
# B3LYP/6-31G(d,p) TD=NStates=6 SCRF=(Solvent=water) Geom=Check Guess=Read

NF: linear response vertical excited states

2 1

第三步：

%mem=55GB
%nprocshared=20
%oldchk=2NFCa1
%chk=2NFCa3
# B3LYP/6-31G(d,p) SCRF=(Solvent=Water,Read) Geom=Check Guess=Read

2NFCa: prepare for state-specific non-eq solvation
by saving the solvent reaction field from the ground state

2 1

NonEq=write

--link1--
%mem=55GB
%nprocshared=20
%chk=2NFCa3
# B3LYP/6-31G(d,p) TD(NStates=6,Root=1) SCRF=(Solvent=Water,ExternalIteration,Read) Geom=Check Guess=Read

2NFCa: read non-eq solvation from ground state and compute energy of the first excited with the state-specific method

2 1

NonEq=read

第四步：

%mem=55GB
%nprocshared=20
%oldchk=2NFCa2
%chk=2NFCa4
# B3LYP/6-31G(d,p) TD=(Read,NStates=6,Root=1) SCRF=(Solvent=water)
  Geom=Modify Guess=Read Opt=ReadFC

2NFCa: excited state opt
Modify geometry to break Cs symmetry
since first excited state is A"

2 1

13 9 10 11 -160.0
53 52 51 50 -188.0

第五步：

%mem=55GB
%nprocshared=20
%oldchk=2NFCa4
%chk=2NFCa5
# B3LYP/6-31G(d,p) TD=(Read,NStates=6,Root=1) Freq
  SCRF=(Solvent=water) Geom=Check Guess=Read

2NFCa excited state freq

2 1

第六步：

%mem=55GB
%nprocshared=20
%oldchk=2NFCa5
%chk=2NFCa6
# B3LYP/6-31G** TD=(read,NStates=6,Root=1) SCRF=(Solvent=Water,ExternalIteration,read) Geom=Check Guess=Read

2NFCa emission state-specific solvation at first excited state optimized geometry

2 1

NonEq=write

作者
Author: sobereva 时间: 2015-6-29 18:33
因为第五步做的是激发态的振动分析。

激发态优化只需要一阶导数，Gaussian的TDDFT支持一阶解析导数，所以快。
激发态振动分析需要二阶导数，Gaussian的TDDFT不支持二阶解析导数，所以很慢，耗时是3N步几何优化的耗时。

作者
Author: lao7 时间: 2015-7-1 09:10
谢谢回复，坚持用一下gaussian。

作者
Author: ZHANGZY 时间: 2016-11-20 22:15
你有q-chem的破解版吗？

作者
Author: lao7 时间: 2016-12-6 14:43

ZHANGZY 发表于 2016-11-20 22:15
你有q-chem的破解版吗？

没有，可以安装试用。

作者
Author: wangzhehyd 时间: 2017-10-23 19:28

lao7 发表于 2016-12-6 14:43
没有，可以安装试用。

能给个试用版的license不，我来crack

作者
Author: lao7 时间: 2017-10-26 11:04

wangzhehyd 发表于 2017-10-23 19:28
能给个试用版的license不，我来crack

试用版可以官方申请吧

作者
Author: wangzhehyd 时间: 2017-10-26 15:09

lao7 发表于 2017-10-26 11:04
试用版可以官方申请吧

http://bbs.keinsci.com/forum.php?mod=viewthread&tid=7074

作者
Author: work.teng 时间: 2018-11-7 12:49
最近也想用这个软件，老板又不同意购买！

作者
Author: sobereva 时间: 2018-11-8 00:13

work.teng 发表于 2018-11-7 12:49
最近也想用这个软件，老板又不同意购买！

一般来说这程序也没什么值得购买的，说服老板应该不容易，本来这程序也没有什么特别吸引人的功能

作者
Author: work.teng 时间: 2018-11-8 11:11

sobereva 发表于 2018-11-8 00:13
一般来说这程序也没什么值得购买的，说服老板应该不容易，本来这程序也没有什么特别吸引人的功能

哈哈，是的，准备换路子了

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