寻找过渡态的方法tsqst2qst3.doc

《寻找过渡态的方法tsqst2qst3.doc》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、过渡态简单点说就是化学反应中瞬间形成的高自由能的不稳定复合物，形象点说就是就是反应势能面上的一个鞍点，所谓的找过渡态也就是要找到这个反应势能面上的这个鞍点，这里以DA反应作为例子简单说下我用gaussian找过渡态的经验。Gaussian里面找过渡态主要有3种方法：TS、QST2以及QST3。这里套用gaussian手册里的说明。TS是进行过渡态而不是局部最小值的优化计算。QST2是使用STQN方法寻找过渡态结构。这个选项需要输入反应物和产物分子结构，先后通过两组连续的标题和分子说明部分定义。注意在这两个结构中的原子顺序必需一致。TS不能

2、和QST2合用。QST3是使用STQN方法寻找过渡态结构。这个选项需要输入反应物，产物和最初的过渡态结构，先后由三组连续的标题和分子说明部分定义。注意在这三个结构中的原子顺序必需一致。TS不能和QST3合用。接下来贴出几个例子来说明下TS、QST2和QST3的输入文件该怎么写%chk=c2h4_c4h6_ts%nproc=1%mem=200mb#opt=(ts,calcfc)freq=noramanb3lyp/6-31g* TitleCardRequired01C            -0.37960200  1.41026100  0

3、.50972400H            -0.06409400  1.04025400  1.48028400H            -0.26589500  2.48067200  0.40089100C            -1.26020300  0.70558300  -0.28508700H            -1.84635800  1.22289600  -1.04415700C            -1.26023800  -0.70552500  -0.28508700H            -1.846

4、41800  -1.22280900  -1.04415900C            -0.37967300  -1.41024700  0.50972300H            -0.06414500  -1.04025600  1.48028200H            -0.26601500  -2.48066300  0.40088700C            1.45648200  0.69084700  -0.25401700H            1.98402400  1.24692000  0.5107910

5、0H            1.29283800  1.24312300  -1.17152100C            1.45645300  -0.69091000  -0.25401800H            1.29278400  -1.24317800  -1.17152200H            1.98396500  -1.24700700  0.51079200qst2%chk=c2h4_c4h6_qst2%nproc=1%mem=200mb#opt=(calcfc,qst2)freq=noramanb3lyp/

6、6-31g*TitleCardRequired01C            -3.17791906  1.13289796  -1.01187682H            -2.64475531  0.20519304  -1.01187682H            -4.24791906  1.13289796  -1.01187682C            -2.50264475  2.30787525  -1.01187682H            -3.03580850  3.23558018  -1.01187682C 

7、           -0.96264475  2.30787525  -1.01187682H            -0.42948100  3.23558018  -1.01187682C            -0.28040135  1.13693077  -1.01187682H            -0.80804382  0.20607449  -1.01187682H            0.78957977  1.14328807  -1.01187682C            -2.59776413  0.98

8、894279  1.10133233H            -3.14031063  0.06692559  1.12200992H            -3.12699887  1.91