| 查看: 4024 | 回復(fù): 17 | ||||
SAYRG金蟲 (小有名氣)
|
[求助]
使用DSLF擬合吸附等溫線時(shí),怎麼樣才能讓得到的參數(shù)為identifiable(參數(shù)估計(jì)值唯一)
|
|||
|
想向大夥兒請(qǐng)教的是使用DSLF模型擬合等溫吸附曲線時(shí),有無其他條件需注意? 才會(huì)使得經(jīng)由擬合所得到的參數(shù)具有再現(xiàn)性、信度、一致性以及identifiable(參數(shù)估計(jì)值唯一)的性質(zhì)。 (即只得到一組參數(shù),而非有多組參數(shù)都可將該曲線擬合得很好的情況。) 目前只能確定限制條件之一就是所有的參數(shù)的數(shù)值都必須大於0。 在閱讀一些文獻(xiàn)後發(fā)現(xiàn),有透過DSLF模型擬合得到參數(shù)的文獻(xiàn)裡,其報(bào)導(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)) 其一,直接報(bào)導(dǎo)使用DSLF模型擬合不同溫度時(shí)該吸附材的某氣體等溫吸附曲線所得到的參數(shù)。 溫度1 q1 b1 q2 b2 n1 n2 溫度2 q1' b1' q2' b2' n1' n2' 溫度3 q1" b1" q2" b2" n1" n2" 其二,不像上頭有多組參數(shù),只會(huì)有一組參數(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 這兩種報(bào)導(dǎo)方式讓我想到兩個(gè)問題! 第一,為何會(huì)有這兩種報(bào)導(dǎo)方式。而第二種方式又是怎麼擬合與分析的呢? 第二,bi肯定會(huì)隨著溫度的變化而有不同,但是q1、q2就其物理意義而言, 應(yīng)該不會(huì)依溫度的高低而有所變動(dòng)。這兩種報(bào)導(dǎo)方式似乎表示著其看待參數(shù) 物理意義的角度有所不同。 而在使用DSLF擬合文獻(xiàn)*1(可在本樓附件一下載)中的化合物1於298K時(shí)吸附CO2的等溫吸附曲線時(shí),無論怎麼擬合,都無法得到其*1之Table S2.所報(bào)導(dǎo)的參數(shù)值。而我用?磉M(jìn)行擬合的軟體為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)) 對(duì)於同一個(gè)吸附劑在不同溫度下吸附某氣體時(shí),會(huì)隨溫度有變化的參數(shù)應(yīng)該就只有the affinity coefficients of site :bi (因?yàn)?bi = bi0 exp(Ei/RT) ; ln(bi) = ln(bi0)+Ei/RT )。因此將q1 q2 (或是q1 q2 n1 n2)看成共享參數(shù),使用DSLF模型以O(shè)rigin中的Global fitting功能進(jìn)行擬合,依舊得不到Table S2.中化合物1吸附CO2時(shí)的DSLF參數(shù)。 因此覺得在擬合時(shí)除了所有參數(shù)都大於0之外, 或許還有其他限制條件需要考慮進(jìn)?懟2艜?huì)蕼厦由擬合得到的參數(shù)具有再現(xiàn)性、一致性以及identifiable(參數(shù)估計(jì)值唯一)的性質(zhì)。 (即只得到一組參數(shù),而非有多組參數(shù)都可將該曲線擬合得很好的情況。) 在使用DSLF模型擬合曲線時(shí)需注意的事情與整個(gè)詳細(xì)過程想和大家請(qǐng)教! *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 |
興趣 |
金蟲 (小有名氣)
|
?????????????????????Щ???}?? ??????????????*1, *2, *3??M???????^???岻????? ???J???@?N?????????(?????·??????D**)???U????????????????????????x?????????????????????Ψ???(identifiable)?? ????????????S??????????Ⅲ???????N????????????????????Origin????1stOpt?e?????????????M??????????????????]?@?N???Ρ? ???????c??????????? **???*1?a???Y???е?Fig. S10.???Table. S6. *1 Mason, J.A., Sumida, K., Herm, Z.R., Krishna, R., Long, J.R., Evaluating Metal-Organic Frameworks for Post-Combustion Carbon Dioxide Capture via Temperature Swing Adsorption, Energy Environ. Sci. 4 (2011) 3030-3040. ???????W? http://pubs.rsc.org/en/Content/ArticleLanding/2011/EE/C1EE01720A ???????n?B?Y http://www.science.uva.nl/resear ... EnvironSci_2011.pdf ?a???Y?????dλ?? http://www.rsc.org/suppdata/ee/c1/c1ee01720a/c1ee01720a.pdf ??http://www.sendspace.com/file/ds79cv ??????????(isosteric heat of gas adsorption, Qst)???????t???S7??S12? *2 Xiang, S.C., He, Y., Zhang, Z., Wu, H., Zhou, W., Krishna, R., Chen, B., Microporous Metal-Organic Framework with Potential for Carbon Dioxide Capture at Ambient Conditions, Nature Communications 3: 954 http://dx.doi.org/doi:10.1038/ncomms1956 (2012) ?????c?a???Y??????????B?Y???d http://www.science.uva.nl/resear ... ons_3_954(2012).pdf *3 He, Y., Zhang, Z., Xiang, S., Fronczek, F.R., Krishna, R., Chen B., A robust doubly interpenetrated metal-organic framework constructed from a novel aromatic tricarboxylate for highly selective separation of small hydrocarbons, Chem. Commun. 48 (2012) 6493-6495. ?????c?a???Y??????????B?Y???d http://www.science.uva.nl/resear ... mmn_2012_withSI.pdf |
金蟲 (小有名氣)
|
??????x*1-1~3???Krishna Group*2?e?c???????P???īI???X??????????l??????????????????M??Global fitting(Global fitting with Parameter Sharing)*3????J??Dual site Langmuir-Freundlich (DSLF)??Dual site Langmuir (DSL)??Langmuir-Freundlich (LF) ......??????????е?Ei ??K??????????????(??Virial method??Clausius-Clapeyron method)??????????y????Global fitting??????S???@??? ?????н?D?ccode?^?K?? ????y?????M???^?????Global fitting?????????????0??E1>E2(in dual-site model)?? ??????в????????????????DSLF?r??ni?@????????????????????? *1-1?????е?D2 *1-1?????е?D3????D?D?f??t??????????? *1-1?????3034??P???????????????????t??????????? *1-1?a???Y???е??S6?cS2?????t???????????@??????c*1-1?????е?D4?t????^???????????? ?????t??????????? ?????t??????????? ?????t??????????? *4?????Y???е?DS8?????t????????????????????????????е?Ei ?K????Virial method ????? *1-1 Mason, J.A., Sumida, K., Herm, Z.R., Krishna, R., Long, J.R., Evaluating Metal-Organic Frameworks for Post-Combustion Carbon Dioxide Capture via Temperature Swing Adsorption, Energy Environ. Sci. 4 (2011) 3030-3040. ???????W? http://pubs.rsc.org/en/Content/ArticleLanding/2011/EE/C1EE01720A ???????n?B?Y http://www.science.uva.nl/research/cr/pdf/Mason_EnergyEnvironSci_2011.pdf ?a???Y?????dλ?? http://www.rsc.org/suppdata/ee/c1/c1ee01720a/c1ee01720a.pdf *1-2 Xiang, S.C., He, Y., Zhang, Z., Wu, H., Zhou, W., Krishna, R., Chen, B., Microporous Metal-Organic Framework with Potential for Carbon Dioxide Capture at Ambient Conditions, Nature Communications 3: 954 http://dx.doi.org/doi:10.1038/ncomms1956 (2012) ?????c?a???Y??????????B?Y???d http://www.science.uva.nl/research/cr/pdf/Xiang_NatureCommunications_3_954(2012).pdf *1-3 He, Y., Zhang, Z., Xiang, S., Fronczek, F.R., Krishna, R., Chen B., A robust doubly interpenetrated metal-organic framework constructed from a novel aromatic tricarboxylate for highly selective separation of small hydrocarbons, Chem. Commun. 48 (2012) 6493-6495. ?????c?a???Y??????????B?Y???d http://www.science.uva.nl/research/cr/pdf/He_ChemCommn_2012_withSI.pdf *2 Krishna Group http://www.science.uva.nl/research/cr/pdf/ *3 Global Fitting with Parameter Sharing In Origin http://www.originlab.com/index.aspx?go=Products/Origin/DataAnalysis/CurveFitting/GlobalFitting http://www.originlab.com/www/helponline/Origin/en/mergedProjects/Tutorial/Tutorial/Global_Fitting.html http://www.originlab.com/www/helponline/Origin/en/UserGuide/Global_fitting_with_parameter_sharing.html http://www.originlab.com/index.aspx?go=Products/Origin/DataAnalysis/CurveFitting/GlobalFitting&pid=1175 In 1stOpt ????????????27~28? 1stOpt?????? *4 Mcdonald, T.M., D'Alessandro, D.M., Krishna, R., Long, J.R., Enhanced carbon dioxide capture upon incorporation of N,N?- dimethylethylenediamine in the metal-organic framework CuBTTri, Chemical Science, 2 (2011) 2022-2028 ???? http://www.science.uva.nl/research/cr/pdf/McDonald_ChemicalScience_2001.pdf ?a???Y?? http://www.rsc.org/suppdata/sc/c1/c1sc00354b/c1sc00354b.pdf |
榮譽(yù)版主 (著名寫手)
彩色掛圖
![]() |
專家經(jīng)驗(yàn): +2 |
|
我上次還下了你的文獻(xiàn)看了一會(huì),我不是搞吸附的,但是沒有看明白。以下是我的一些理解: 對(duì)應(yīng)不同的溫度T下,你應(yīng)該有一組q關(guān)于p的數(shù)據(jù)(你的公式里是Q和x),根據(jù)這組數(shù)據(jù),以及公式Q = (q1*b1*x)/(1 + b1*x) + (q2*b2*x)/(1 + b2*x) 可以擬合得到該T溫度下對(duì)應(yīng)的q1,b1,q2,b2。 然后對(duì)于不同的溫度,q1和q2應(yīng)該是相同的,但是b1和b2是不同的,也就是說q1和q2與T無關(guān),b1和b2與T有關(guān)。 那么假如你有兩個(gè)溫度T1和T2,在擬合時(shí)應(yīng)該是q1=q1(T1)=q1(T2),q2=q2(T1)=q2(T2)作為公共參數(shù)擬合,b1(T1),b1(T2),b2(T1),b2(T2)則分別擬合。 這樣擬合完了之后應(yīng)該有1個(gè)q1和1個(gè)q2,2個(gè)b1,2個(gè)b2。 對(duì)應(yīng)n個(gè)溫度T是 1個(gè)q1和1個(gè)q2,n個(gè)b1,n個(gè)b2。 至此,q1和q2可以確定,另外還得到一組b1和一組b2(各分別與n個(gè)溫度T對(duì)應(yīng))。 然后將b1和b2分別對(duì)T擬合,根據(jù)公式b=b0 exp(E/RT)。R已知,那么b0和E未知,可以擬合出來b0和E。 但是,進(jìn)一步的理解,E應(yīng)該是不變的,它代表的是吸附所需要的一個(gè)特征能量,所以應(yīng)該是能從圖中看出來的。 |

|
本帖內(nèi)容被屏蔽 |
榮譽(yù)版主 (著名寫手)
彩色掛圖
![]() |
專家經(jīng)驗(yàn): +2 |

|
本帖內(nèi)容被屏蔽 |
金蟲 (小有名氣)
送鮮花一朵|
不知道能不能再向你請(qǐng)教一些問題? 1. " 擬合的時(shí)候不同溫度的q1,q2應(yīng)該設(shè)為常數(shù) " 我個(gè)人非常同意這個(gè)看法,q1,q2應(yīng)該不會(huì)隨著溫度而有變化。 a. q1, q2是先以其他方式求得的嗎?那方式又是什麼呢? b. 還是在擬合過程當(dāng)中將q1, q2視作共享參數(shù),進(jìn)而透過Global fitting得到? c. 我想了解怎麼得到q1, q2的。 2. " 比較好的方法是通過分子模擬估計(jì)吸附熱等參數(shù),確定數(shù)據(jù)可能確定范圍再代入擬合。 " a. 目前有辦法透過維里(Virial)方程式計(jì)算零覆蓋率(zero coverage)與不同覆蓋率時(shí)的等量吸附熱(Qst, Isosteric heat or absorption of heat)。但是不知道如何推估E1和E2?煞窈臀艺f說? |
金蟲 (小有名氣)
|
我將目前知道一些線索就寫在下方了! 1. 所有的參數(shù)都大於0 2. 在 Single site的模型裡,b=b0 exp(E/RT) ; ln(b)=ln(b0)+E/RT E是否就是吸附劑吸附某氣體於零覆蓋率時(shí)的等量吸附熱(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是否就是吸附劑吸附某氣體於零覆蓋率時(shí)的等量吸附熱(zero-coverage isosteric heat of gas adsorption)? E1和E2又是怎麼得到呢? 3. q 或是 q1, q2根據(jù)其物理意義應(yīng)該不會(huì)隨著溫度變動(dòng)而改變。 a. q1, q2是先以其他方式求得的嗎?那方式又是什麼呢? b. 還是在擬合過程當(dāng)中將q1, q2視作共享參數(shù),進(jìn)而透過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 |
|
本帖內(nèi)容被屏蔽 |
金蟲 (小有名氣)
|
?L???Origin ?? Global fitting????????У??????F?o?????(?????)??????? ???????J??Global fitting????????]?e?????????Щ????????????? ??????????*1??????o?????c?c???h?????x??? *1 http://www.gaoyang168.com/bbs/viewthread.php?tid=5201919&page=1#pid7 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0703化學(xué)調(diào)劑 +3 | 18889395102 2026-03-18 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +6 | Ymlll 2026-03-18 | 9/450 |
|
|
[考研] 一志愿985,本科211,0817化學(xué)工程與技術(shù)319求調(diào)劑 +10 | Liwangman 2026-03-15 | 10/500 |
|
|
[考研] 267一志愿南京工業(yè)大學(xué)0817化工求調(diào)劑 +10 | SUICHILD 2026-03-12 | 10/500 |
|
|
[考研] 330求調(diào)劑 +3 | 小材化本科 2026-03-18 | 3/150 |
|
|
[考研] 344求調(diào)劑 +6 | knight344 2026-03-16 | 7/350 |
|
|
[考研] 一志愿中國(guó)海洋大學(xué),生物學(xué),301分,求調(diào)劑 +4 | 1孫悟空 2026-03-17 | 4/200 |
|
|
[考研] 295求調(diào)劑 +3 | 一志愿京區(qū)211 2026-03-18 | 5/250 |
|
|
[考研] 0854可跨調(diào)劑,一作一項(xiàng)核心論文五項(xiàng)專利,省、國(guó)級(jí)證書40+數(shù)一英一287 +8 | 小李0854 2026-03-16 | 8/400 |
|
|
[考研] 302求調(diào)劑 +10 | 呼呼呼。。。。 2026-03-17 | 10/500 |
|
|
[考研] 303求調(diào)劑 +4 | 睿08 2026-03-17 | 6/300 |
|
|
[考研] 環(huán)境工程調(diào)劑 +8 | 大可digkids 2026-03-16 | 8/400 |
|
|
[考研] 268求調(diào)劑 +7 | 好運(yùn)連綿不絕 2026-03-12 | 8/400 |
|
|
[考研] 268求調(diào)劑 +8 | 一定有學(xué)上- 2026-03-14 | 9/450 |
|
|
[考研] 308求調(diào)劑 +4 | 是Lupa啊 2026-03-16 | 4/200 |
|
|
[考研] 085601求調(diào)劑 +4 | Du.11 2026-03-16 | 4/200 |
|
|
[考研] 26考研求調(diào)劑 +6 | 丶宏Sir 2026-03-13 | 6/300 |
|
|
[考研] 材料301分求調(diào)劑 +5 | Liyouyumairs 2026-03-12 | 5/250 |
|
|
[考研] 308求調(diào)劑 +3 | 是Lupa啊 2026-03-12 | 3/150 |
|
|
[考博] 福州大學(xué)楊黃浩課題組招收2026年專業(yè)學(xué)位博士研究生,2026.03.20截止 +3 | Xiangyu_ou 2026-03-12 | 3/150 |
|