請問我這個(gè)化合物該用什么方法計(jì)算,急!重獎(jiǎng)!
有機(jī)單晶,DFT B3LYP with a 6-31G(d) basis set 計(jì)算優(yōu)化了得到了構(gòu)象和真實(shí)的單晶構(gòu)象做比較,投了篇文章。分子內(nèi)存在C-H…pi contact作用。
審稿人說分子如存在C-H…pi contact,就不能用這計(jì)算方法優(yōu)化計(jì)算,具體信息如下。 我該怎么辦呢?用什么方法計(jì)算?看了些文獻(xiàn)可自己還是不知道怎么計(jì)算啊。非本專業(yè)麻煩各位大家?guī)臀铱纯戳?br />
1 General methodology for QM calculations. Authors should give more insights in computational details. Was located minima confirmed by frequency calculations. Have they also done wave function stability calculations for optimized structures? If not those calculations have to be done in order to fully characterize located states.
Another suggestion on how to find the conformations closest to the X-ray conformations is to start from X-ray pdb structure and see how it will be optimized within QM.
2 Chosen level for calculations. Authors are paying attention to CH…pi , assuming that they are main forces stabilizing the compound ("The conformations of molecules investigated are determined mainly by hydrogen bonds and electrostatic interactions"). Indeed C-H…pi contact is mainly established by dispersion forces rather than electrostatic [1]. From this point of view chosen DFT method is not the best choice to study those compounds. DFT has its limitations to describe dispersion interaction satisfactory and requires introduction of dispersion correction. That's probably why simplistic DFT approach (B3LYP/6-31G (d)) gave poor results. Authors should approach the question of choosing right level of computations carefully.There are a couple of reviews available on this problem [2-3]. Additional calculations at higher level are reasonable and needed to complete the study.
我的化合物結(jié)構(gòu)式

[ Last edited by baiyuefei on 2011-6-16 at 20:32 ]
返回小木蟲查看更多
京公網(wǎng)安備 11010802022153號(hào)
1、審稿人問你是否對優(yōu)化的結(jié)構(gòu)做了頻率分析,確保沒有虛頻來證實(shí)這是個(gè)穩(wěn)定結(jié)構(gòu);是否做了波函數(shù)穩(wěn)定測試,也可以用自選污染S^2來確認(rèn)一下,如果自選污染比較大,嘗試用非限制開竅層計(jì)算雙自由自單重態(tài);嘗試用單晶結(jié)構(gòu)做初始結(jié)構(gòu)再優(yōu)化;
2、普通的DFT/6-31g*無論是方法還是基組都不足以表征C-H…pi 作用,可以嘗試用lc-b3lyp或者pbepbe或者dft-d等描述長程作用和dispersion forces ,基組連彌散都沒有,更不合適,可加大到6-311++g**再優(yōu)化。
估計(jì)就是這個(gè)意思吧。
另外評審的[1] [2-3]應(yīng)該是給你推薦的文獻(xiàn)吧。。那就更應(yīng)該認(rèn)真讀一讀,看看別人類似的體系是怎么做的了,而不是到這里來求助啊。這樣的評審真的不錯(cuò),不僅告訴你哪里不足,還給你指點(diǎn)怎么改正,相當(dāng)不錯(cuò)了。其他懶惰或者輕視你的,直接指出不足就斃了文章了。
是有兩篇文獻(xiàn),但看了下,不太明白。其實(shí)文章里的計(jì)算是我找好朋友幫我弄的。非本專業(yè),請諒解。
您能具體點(diǎn)我該用哪種計(jì)算方法優(yōu)化我的分子嗎?麻煩您了
,
try DFT uB3LYP with a 6-311++G(d,p) basis set