 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 4039 | 回復(fù): 17 | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | |||
SAYRG金蟲 (小有名氣)
|
[求助]
使用DSLF擬合吸附等溫線時，怎麼樣才能讓得到的參數(shù)為identifiable(參數(shù)估計值唯一)
|
||
|
想向大夥兒請教的是使用DSLF模型擬合等溫吸附曲線時，有無其他條件需注意？ 才會使得經(jīng)由擬合所得到的參數(shù)具有再現(xiàn)性、信度、一致性以及identifiable(參數(shù)估計值唯一)的性質(zhì)。 (即只得到一組參數(shù)，而非有多組參數(shù)都可將該曲線擬合得很好的情況。) 目前只能確定限制條件之一就是所有的參數(shù)的數(shù)值都必須大於0。 在閱讀一些文獻後發(fā)現(xiàn)，有透過DSLF模型擬合得到參數(shù)的文獻裡，其報導(dǎo)參數(shù)的方式大致上有兩種。以使用DSLF模型的為例。 DSLF Q = (q1*b1*x^n1)/(1 + b1*x^n1) + (q2*b2*x^n2)/(1 + b2*x^n2) 或 Q = (q1*b1*x^(1/n1))/(1 + b1*x^(1/n1)) + (q2*b2*x^(1/n2))/(1 + b2*x^(1/n2)) 其一，直接報導(dǎo)使用DSLF模型擬合不同溫度時該吸附材的某氣體等溫吸附曲線所得到的參數(shù)。 溫度1 q1 b1 q2 b2 n1 n2 溫度2 q1' b1' q2' b2' n1' n2' 溫度3 q1" b1" q2" b2" n1" n2" 其二，不像上頭有多組參數(shù)，只會有一組參數(shù)。 q1 b10 E1 n1 q2 b20 E2 n2 b1=b10 exp(E1/RT), ln(b1) = ln(b10)+E/RT b2=b20 exp(E2/RT), ln(b2) = ln(b20)+E/RT 這兩種報導(dǎo)方式讓我想到兩個問題！ 第一，為何會有這兩種報導(dǎo)方式。而第二種方式又是怎麼擬合與分析的呢？ 第二，bi肯定會隨著溫度的變化而有不同，但是q1、q2就其物理意義而言， 應(yīng)該不會依溫度的高低而有所變動。這兩種報導(dǎo)方式似乎表示著其看待參數(shù) 物理意義的角度有所不同。 而在使用DSLF擬合文獻*1(可在本樓附件一下載)中的化合物1於298K時吸附CO2的等溫吸附曲線時，無論怎麼擬合，都無法得到其*1之Table S2.所報導(dǎo)的參數(shù)值。而我用?磉M行擬合的軟體為Origin 8.6版。(化合物1之CO2以及N2的單組成吸附資料在本樓附件二或是可取用附件一的Table S3.，而我自己擬合的結(jié)果在本樓末) 另外有想過以下事情。 Q = (q1*b1*x^n1)/(1 + b1*x^n1) + (q2*b2*x^n2)/(1 + b2*x^n2) 或 Q = (q1*b1*x^(1/n1))/(1 + b1*x^(1/n1)) + (q2*b2*x^(1/n2))/(1 + b2*x^(1/n2)) 對於同一個吸附劑在不同溫度下吸附某氣體時，會隨溫度有變化的參數(shù)應(yīng)該就只有the affinity coefficients of site ：bi (因為 bi = bi0 exp(Ei/RT) ; ln(bi) = ln(bi0)+Ei/RT )。因此將q1 q2 (或是q1 q2 n1 n2)看成共享參數(shù)，使用DSLF模型以O(shè)rigin中的Global fitting功能進行擬合，依舊得不到Table S2.中化合物1吸附CO2時的DSLF參數(shù)。 因此覺得在擬合時除了所有參數(shù)都大於0之外， 或許還有其他限制條件需要考慮進?懟２艜溝糜蓴M合得到的參數(shù)具有再現(xiàn)性、一致性以及identifiable(參數(shù)估計值唯一)的性質(zhì)。 (即只得到一組參數(shù)，而非有多組參數(shù)都可將該曲線擬合得很好的情況。) 在使用DSLF模型擬合曲線時需注意的事情與整個詳細過程想和大家請教！ *1 Enhanced CO2 Binding Affinity of a High-Uptake rht-Type Metal−Organic Framework Decorated with Acylamide Groups Supporting Information中的Figure S10與Table S2 http://pubs.acs.org/doi/suppl/10.1021/ja110042b |
回復(fù)此樓» 本帖附件資源列表
-
歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺，不對該內(nèi)容負責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識產(chǎn)權(quán)問題，其責(zé)任在于用戶本人，如對版權(quán)有異議，請聯(lián)系郵箱：xiaomuchong@tal.com - 附件 1 : ja110042b_si_001.pdf
- 附件 2 : ja110042b_SI_table-S3.xls
2012-11-19 20:42:51, 708.53 K
2012-11-19 21:00:46, 53.5 K
» 收錄本帖的淘帖專輯推薦
 興趣 |
» 猜你喜歡
請問還有沒有用Latex寫文章的小伙伴們？
已經(jīng)有0人回復(fù)
-
光學(xué)工程學(xué)碩調(diào)劑信息
已經(jīng)有32人回復(fù)
-
物理學(xué)I論文潤色/翻譯怎么收費?
已經(jīng)有153人回復(fù)
-
歡迎加入課題組
已經(jīng)有0人回復(fù)
-
散金幣，求好運，祝面上順利！
已經(jīng)有78人回復(fù)
-
華東師范大學(xué)芯片設(shè)計徐瓏真實水平咋樣
已經(jīng)有1人回復(fù)
-
調(diào)劑
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復(fù)
-
中德博士后2026項目5月30日截止- 國內(nèi)博士后入站、德國兩年全職
已經(jīng)有0人回復(fù)
-
面向全球申請人的中德博士后2026項目5月30日截止- 國內(nèi)博后入站德國兩年全職
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有2人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 一些關(guān)於孔洞材料吸附氣體時，吸附量轉(zhuǎn)換、吸附選擇性計算、等溫吸附曲線擬合的問題
已經(jīng)有1人回復(fù)
- 如何在Origin等擬合軟件上以DSLF等模型擬合多條吸附等溫線(見帖內(nèi)文檔)
已經(jīng)有2人回復(fù)
SAYRG
金蟲 (小有名氣)
- 應(yīng)助: 8 (幼兒園)
- 金幣: 1456.7
- 紅花: 1
- 帖子: 230
- 在線: 181.1小時
- 蟲號: 1565042
- 注冊: 2012-01-05
- 專業(yè): 理論和計算化學(xué)
★
送鮮花一朵
csgt0: 金幣+1, 鼓勵交流 2012-11-30 09:29:19
送鮮花一朵csgt0: 金幣+1, 鼓勵交流 2012-11-30 09:29:19
|
不知道能不能再向你請教一些問題？ 1. " 擬合的時候不同溫度的q1，q2應(yīng)該設(shè)為常數(shù) " 我個人非常同意這個看法，q1,q2應(yīng)該不會隨著溫度而有變化。 a. q1, q2是先以其他方式求得的嗎？那方式又是什麼呢？ b. 還是在擬合過程當(dāng)中將q1, q2視作共享參數(shù)，進而透過Global fitting得到？ c. 我想了解怎麼得到q1, q2的。 2. " 比較好的方法是通過分子模擬估計吸附熱等參數(shù)，確定數(shù)據(jù)可能確定范圍再代入擬合。 " a. 目前有辦法透過維里(Virial)方程式計算零覆蓋率(zero coverage)與不同覆蓋率時的等量吸附熱(Qst, Isosteric heat or absorption of heat)。但是不知道如何推估E1和E2。可否和我說說？ |
4樓2012-11-21 18:29:12
csgt0
榮譽版主 (著名寫手)
彩色掛圖
-
專家經(jīng)驗: +2 - 計算強帖: 2
- 應(yīng)助: 367 (碩士)
- 貴賓: 0.457
- 金幣: 7807
- 散金: 8897
- 紅花: 20
- 沙發(fā): 7
- 帖子: 2809
- 在線: 515.8小時
- 蟲號: 449046
- 注冊: 2007-11-02
- 專業(yè): 金屬非晶態(tài)、準(zhǔn)晶和納米晶
- 管轄: 計算模擬
2樓2012-11-20 11:59:07
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
csgt0: 金幣+2, 多謝應(yīng)助 2012-11-21 11:27:15
SAYRG: 金幣+10, ★★★很有幫助, 很有收穫！非常謝謝！ 2012-11-21 14:25:19
感謝參與，應(yīng)助指數(shù) +1
csgt0: 金幣+2, 多謝應(yīng)助 2012-11-21 11:27:15
SAYRG: 金幣+10, ★★★很有幫助, 很有收穫！非常謝謝！ 2012-11-21 14:25:19
|
本帖內(nèi)容被屏蔽 |
» 本帖已獲得的紅花（最新10朵）
3樓2012-11-21 00:17:30
SAYRG
金蟲 (小有名氣)
- 應(yīng)助: 8 (幼兒園)
- 金幣: 1456.7
- 紅花: 1
- 帖子: 230
- 在線: 181.1小時
- 蟲號: 1565042
- 注冊: 2012-01-05
- 專業(yè): 理論和計算化學(xué)
|
我將目前知道一些線索就寫在下方了！ 1. 所有的參數(shù)都大於0 2. 在 Single site的模型裡，b=b0 exp(E/RT) ; ln(b)=ln(b0)+E/RT E是否就是吸附劑吸附某氣體於零覆蓋率時的等量吸附熱(zero-coverage isosteric heat of gas adsorption)？ 在 Dual site的模型裡， b1=b10 exp(E1/RT) ; ln(b1)=ln(b10)+E1/RT b2=b20 exp(E1/RT) ; ln(b2)=ln(b20)+E2/RT b1b2=b10b20 exp((E1+E2)/RT) ; ln(b1b2)=ln(b10b20)+(E1+E2)/RT E1+E2是否就是吸附劑吸附某氣體於零覆蓋率時的等量吸附熱(zero-coverage isosteric heat of gas adsorption)？ E1和E2又是怎麼得到呢？ 3. q 或是 q1, q2根據(jù)其物理意義應(yīng)該不會隨著溫度變動而改變。 a. q1, q2是先以其他方式求得的嗎？那方式又是什麼呢？ b. 還是在擬合過程當(dāng)中將q1, q2視作共享參數(shù)，進而透過Global fitting得到？ LF Q = (q*b*x^n)/(1 + b*x^n) 或 Q = (q*b*x^(1/n))/(1 + b*x^(1/n)) b=b0 exp(E/RT) ; ln(b)=ln(b0)+E/RT DSLF Q = (qm1*b1*x^n1)/(1 + b1*x^n1) + (qm2*b2*x^n2)/(1 + b2*x^n2) 或 Q = (qm1*b1*x^(1/n1))/(1 + b1*x^(1/n1)) + (qm2*b2*x^(1/n2))/(1 + b2*x^(1/n2)) b1=b10 exp(E1/RT) ; ln(b1)=ln(b10)+E1/RT b2=b20 exp(E1/RT) ; ln(b2)=ln(b20)+E2/RT b1b2=b10b20 exp((E1+E2)/RT) ; ln(b1b2)=ln(b10b20)+(E1+E2)/RT DSL Q = (qm1*b1*x)/(1 + b1*x) + (qm2*b2*x)/(1 + b2*x) b1=b10 exp(E1/RT) ; ln(b1)=ln(b10)+E1/RT b2=b20 exp(E1/RT) ; ln(b2)=ln(b20)+E2/RT b1b2=b10b20 exp((E1+E2)/RT) ; ln(b1b2)=ln(b10b20)+(E1+E2)/RT |
5樓2012-11-22 17:13:13
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 335求調(diào)劑 +4 | yuyu宇 2026-03-23 | 5/250 |
|
|---|---|---|---|---|
|
|
[考研] 299求調(diào)劑 +5 | 某某某某位 2026-03-21 | 5/250 |
|
|
[考研] 材料專碩找調(diào)劑 +3 | 哈哈哈吼吼吼哈 2026-03-23 | 3/150 |
|
|
|
[考研] 291 求調(diào)劑 +4 | 化工2026屆畢業(yè)?/a> 2026-03-21 | 5/250 |
|
|
|
[考研] 276求調(diào)劑 +3 | YNRYG 2026-03-21 | 4/200 |
|
|
|
[考研] 280分求調(diào)劑 一志愿085802 +4 | PUMPT 2026-03-22 | 7/350 |
|
|
|
[考研] 311求調(diào)劑 +6 | 冬十三 2026-03-18 | 6/300 |
|
|
|
[考研] 求調(diào)劑一志愿海大，0703化學(xué)學(xué)碩304分，有大創(chuàng)項目，四級已過 +6 | 幸運哩哩 2026-03-22 | 10/500 |
|
|
|
[考研] 生物學(xué)調(diào)劑 +5 | Surekei 2026-03-21 | 5/250 |
|
|
|
[考研] 318求調(diào)劑 +4 | plum李子 2026-03-21 | 7/350 |
|
|
|
[考博] 招收博士1-2人 +3 | QGZDSYS 2026-03-18 | 4/200 |
|
|
|
[考研] 考研調(diào)劑 +3 | 呼呼？~+123456 2026-03-21 | 3/150 |
|
|
|
[考研] 0805 316求調(diào)劑 +3 | 大雪深藏 2026-03-18 | 3/150 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | 妮妮ninicgb 2026-03-21 | 4/200 |
|
|
|
[考研] 302求調(diào)劑 +12 | 呼呼呼。。。。 2026-03-17 | 12/600 |
|
|
|
[考研] 296求調(diào)劑 +6 | www_q 2026-03-18 | 10/500 |
|
|
|
[考研] 295材料求調(diào)劑，一志愿武漢理工085601專碩 +5 | Charlieyq 2026-03-19 | 5/250 |
|
|
|
[考研] 一志愿南理工085701環(huán)境302求調(diào)劑院校 +3 | 葵梓衛(wèi)隊 2026-03-20 | 3/150 |
|
|
|
[考研] 086500 325 求調(diào)劑 +3 | 領(lǐng)帶小熊 2026-03-19 | 3/150 |
|
|
|
[考研] 材料學(xué)碩318求調(diào)劑 +5 | February_Feb 2026-03-19 | 5/250 |
|