| 2 | 1/1 | 返回列表 |
| 查看: 2131 | 回復(fù): 1 | ||||
sfycong新蟲 (初入文壇)
|
[交流]
用Dmol3里的Cosmo模型計(jì)算sigma——profile的問題 已有1人參與
|
|
我現(xiàn)在用的是MS7.0,最近一直在用dmol3里的Cosmo計(jì)算分子的sigma--profile,但發(fā)現(xiàn)現(xiàn)在的新版本的Cosmo里的計(jì)算方法默認(rèn)為ibs,這和老方法有什么區(qū)別,和老方法相比算出的結(jié)果會(huì)不會(huì)有很大的差別? 1、我用MS7.0算水分子的sigma--profile,沒有修改input文件,計(jì)算結(jié)果和文獻(xiàn)上用guass和Cosmothermo軟件算出來的有很大的差別。那到底應(yīng)不應(yīng)該修改input文件,我列出我計(jì)算水的sigma--profile修改的input文件,如下: # Task parameters Calculate energy Symmetry off Max_memory 10480 File_usage smart # Electronic parameters Spin_polarization restricted Charge 0 Basis dnp basis_version basfile_v4.4 Pseudopotential none Functional vwn-bp Aux_density octupole Integration_grid fine Occupation fermi Cutoff_Global 5.5000 angstrom Scf_density_convergence 1.0000e-005 Scf_charge_mixing 2.0000e-001 Scf_iterations 500 Scf_diis 6 pulay Cosmo ibs COSMO_Dielectric 78.5400 Water # COSMO Environment Keywords Cosmo_Grid_Size 1082 Cosmo_Segments 92 Cosmo_Solvent_Radius 1.300000 Cosmo_A-Matrix_Cutoff 7.000000 Cosmo_RadCorr_Incr 0.000000 Cosmo_Radii(units) angstrom 1 1.300000 2 1.720000 3 1.300000 這樣的input文件有沒有什么問題?算出來的sigma--profile的圖片已上傳。 2.我要用Cosmo計(jì)算一些金屬離子的sigma--profile,是不是直接建單個(gè)離子的模型直接計(jì)算還是要考慮離子的水化問題建水化離子模型在來計(jì)算? 3.文獻(xiàn)上列出的Cosmo半徑,比如說:C:2.00A,H:1.30A,O:1.72A,Li+=1.750A,Na+=1.862A,這些半徑該怎么用? 我現(xiàn)在就是因?yàn)樗鉺igma--profile的問題犯愁啊,不知問題出在哪里,真心希望得到各位大神的指導(dǎo),本人將不勝感激。 |
專家顧問 (著名寫手)
Ab Initio Amateur
![]() |
專家經(jīng)驗(yàn): +224 |
|
1、有區(qū)別,兩種方法再劃分溶劑空洞的時(shí)候形狀應(yīng)該是有些不同的。除了溶劑化不同,4.0以前、4.1和5.5以后的版本計(jì)算方法也稍有差異,請(qǐng)注意help文檔中的MS40_Compat以及MS41_Compat關(guān)鍵詞。 2、一些金屬離子與水分子的作用是比較明顯的化學(xué)鍵,具有明確的方向性,這種通過“自洽反應(yīng)場”將水的構(gòu)型、物性等平均化卻不具備這樣的方向性,至少內(nèi)氛配位的幾個(gè)水分子要考慮顯式。 3、這個(gè)半徑要是沒記錯(cuò)的話應(yīng)該是決定孔洞面積的吧。實(shí)際溶劑化模型的操作無非是將溶質(zhì)外部罩上一個(gè)孔洞,計(jì)算這個(gè)孔洞和溶質(zhì)的電荷作用。Gaussian的溶劑化模型還是有有些不同的,雖然cpcm和cosmol是一樣的,但是很多內(nèi)置參數(shù)的區(qū)別因?yàn)槭诸^資料有限,我不是很了解。印象中Gaussian的溶劑參數(shù)是指定計(jì)算級(jí)別確定的,如uff、hf和pbe0水平的溶劑信息。 至于這些半徑你可以進(jìn)行指定,如你的input文件中這一行: Cosmo_Radii(units) angstrom 下面指定原子的編號(hào)及其對(duì)應(yīng)的溶劑化半徑即可。 另外想要精確的話 Cosmo_Grid_Size 1082 Cosmo_Segments 92 格點(diǎn)和面元片段也要充分。 |

| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085600材料與化工調(diào)劑 324分 +9 | llllkkkhh 2026-03-18 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 一志愿華中科技大學(xué),080502,354分求調(diào)劑 +4 | 守候夕陽CF 2026-03-18 | 4/200 |
|
|
[考研] 328求調(diào)劑,英語六級(jí)551,有科研經(jīng)歷 +3 | 生物工程調(diào)劑 2026-03-17 | 7/350 |
|
|
[考研] 一志愿中國海洋大學(xué),生物學(xué),301分,求調(diào)劑 +4 | 1孫悟空 2026-03-17 | 4/200 |
|
|
[考研] 085600材料與化工 +5 | 安全上岸! 2026-03-16 | 5/250 |
|
|
[教師之家] 焦慮 +8 | 水冰月月野兔 2026-03-13 | 12/600 |
|
|
[考研] 280求調(diào)劑 +6 | 咕嚕曉曉 2026-03-18 | 7/350 |
|
|
[考研] 0703化學(xué)求調(diào)劑 總分331 +3 | ZY-05 2026-03-13 | 3/150 |
|
|
[考研] 環(huán)境工程調(diào)劑 +8 | 大可digkids 2026-03-16 | 8/400 |
|
|
[基金申請(qǐng)] 被我言中:新模板不強(qiáng)調(diào)格式了,假專家開始管格式了 +4 | beefly 2026-03-14 | 4/200 |
|
|
[考研] 278求調(diào)劑 +3 | Yy7400 2026-03-13 | 3/150 |
|
|
[考研] 080500,材料學(xué)碩302分求調(diào)劑學(xué)校 +4 | 初識(shí)可樂 2026-03-14 | 5/250 |
|
|
[考研] 289求調(diào)劑 +4 | 這么名字咋樣 2026-03-14 | 6/300 |
|
|
[考研] 復(fù)試調(diào)劑 +3 | 呼呼?~+123456 2026-03-14 | 3/150 |
|
|
[考研] 308 085701 四六級(jí)已過求調(diào)劑 +7 | 溫喬喬喬喬 2026-03-12 | 14/700 |
|
|
[考研] 0703化學(xué)求調(diào)劑 +7 | 綠豆芹菜湯 2026-03-12 | 7/350 |
|
|
[考研] 295求調(diào)劑 +3 | 小匕仔汁 2026-03-12 | 3/150 |
|
|
[考研] 289求調(diào)劑 +3 | 李政瑩 2026-03-12 | 3/150 |
|
|
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +3 | d如愿上岸 2026-03-13 | 3/150 |
|
|
[考研] 290求調(diào)劑 +3 | ADT 2026-03-13 | 3/150 |
|