طريقة جوباك

تتنبأ طريقة جوباك ، والتي تُعرف غالبًا باسم طريقة جوباك-ريد ، بإحدى عشرة خاصية ديناميكية حرارية مهمة وشائعة الاستخدام للمكونات النقية انطلاقًا من التركيب الجزيئي فقط. سُميت هذه الطريقة نسبةً إلى كيفن ج. جوباك عام 1984 [ 1 ] ، ثم طُوّرت لاحقًا بالتعاون مع روبرت س. ريد. [ 2 ] تُعد طريقة جوباك امتدادًا لطريقة ليدرسن [ 3 ] ، وتستخدم مجموعات وصيغًا ومعاملات متشابهة جدًا للخصائص الثلاث التي دعمتها طريقة ليدرسن سابقًا ( درجة الحرارة الحرجة ، والضغط الحرج ، والحجم الحرج).

قام جوباك وريد بتوسيع نطاق الخصائص المدعومة، وإنشاء معلمات جديدة، وتعديل صيغ طريقة ليدرسن القديمة بشكل طفيف.

المبادئ الأساسية

أسلوب المساهمة الجماعية

مبدأ أسلوب المساهمة الجماعية

طريقة جوباك هي طريقة مساهمة جماعية . تستخدم هذه الأنواع من الطرق المعلومات الهيكلية الأساسية للجزيء الكيميائي، مثل قائمة بالمجموعات الوظيفية البسيطة، وتضيف معلمات إلى هذه المجموعات الوظيفية، وتحسب الخصائص الحرارية الفيزيائية وخصائص النقل كدالة لمجموع معلمات المجموعة.

يفترض جوباك عدم وجود تفاعلات بين المجموعات، ولذلك يستخدم فقط المساهمات الجمعية دون أي مساهمات للتفاعلات بينها. أما طرق مساهمات المجموعات الأخرى، وخاصةً طرق مثل UNIFAC التي تُقدّر خصائص الخليط كمعاملات النشاط، فتستخدم كلاً من معلمات المجموعات الجمعية البسيطة ومعلمات تفاعل المجموعات. وتكمن الميزة الكبرى لاستخدام معلمات المجموعات البسيطة فقط في قلة عدد المعلمات المطلوبة. بينما يزداد عدد معلمات تفاعل المجموعات المطلوبة بشكل كبير مع ازدياد عدد المجموعات (1 لمجموعتين، 3 لثلاث مجموعات، 6 لأربع مجموعات، 45 لعشر مجموعات، ويتضاعف هذا العدد إذا لم تكن التفاعلات متناظرة).

تسع من هذه الخصائص عبارة عن قيم مفردة لا تتأثر بدرجة الحرارة، ويتم تقديرها في الغالب بجمع بسيط لمساهمات المجموعات بالإضافة إلى حد إضافي. اثنتان من الخصائص المقدرة تعتمدان على درجة الحرارة: السعة الحرارية للغاز المثالي واللزوجة الديناميكية للسوائل. تستخدم معادلة السعة الحرارية أربعة معاملات، بينما تستخدم معادلة اللزوجة معاملين فقط. في كلتا الحالتين، يتم حساب معاملات المعادلة من خلال مساهمات المجموعات.

نقاط قوة وضعف النموذج

نقاط القوة

يعود الفضل في شعبية ونجاح طريقة جوباك بشكل أساسي إلى قائمة المجموعة الواحدة لجميع الخصائص. وهذا يسمح بالحصول على جميع الخصائص الإحدى عشرة المدعومة من خلال تحليل واحد للبنية الجزيئية.

تستخدم طريقة جوباك أيضًا مخططًا بسيطًا جدًا وسهل التحديد للمجموعات ، مما يجعل الطريقة قابلة للاستخدام للأشخاص الذين لديهم معرفة كيميائية أساسية فقط.

نقاط الضعف

الأخطاء المنهجية لطريقة جوباك (نقطة الغليان العادية)

أظهرت التطورات الحديثة في أساليب التقدير [ 4 ] [ 5 ] أن جودة طريقة جوباك محدودة. وقد ذكر المؤلفون الأصليون في ملخص المقالة الأصلية: "لا ندّعي دقة عالية، ولكن الطرق المقترحة غالبًا ما تكون بنفس دقة التقنيات الشائعة الاستخدام اليوم أو أكثر دقة منها".

لا تشمل قائمة المجموعات العديد من الجزيئات الشائعة بشكل كافٍ. فعلى وجه الخصوص، لا يتم التمييز بين المركبات العطرية والمكونات الحلقية العادية. وهذه مشكلة خطيرة لأن المكونات العطرية والأليفاتية تختلف اختلافًا كبيرًا.

كانت قاعدة البيانات التي استخدمها جوباك وريد للحصول على معلمات المجموعات صغيرة نسبيًا، ولم تغطِ سوى عدد محدود من الجزيئات المختلفة. وقد تحققت أفضل تغطية لنقاط الغليان العادية (438 مكونًا)، بينما كانت أسوأ تغطية لحرارة الانصهار (155 مكونًا). أما التطورات الحالية التي تستخدم بنوك البيانات، مثل بنك بيانات دورتموند أو قاعدة بيانات DIPPR، فتتمتع بتغطية أوسع بكثير.

تُظهر الصيغة المستخدمة للتنبؤ بنقطة الغليان العادية مشكلة أخرى. فقد افترض جوباك مساهمة ثابتة للمجموعات المضافة في السلاسل المتجانسة كالألكانات ، وهو ما لا يصف السلوك الحقيقي لنقاط الغليان العادية بدقة. [ 6 ] وبدلاً من المساهمة الثابتة، يجب تطبيق انخفاض في هذه المساهمة مع ازدياد عدد المجموعات. وتؤدي الصيغة المختارة في طريقة جوباك إلى انحرافات كبيرة للجزيئات الكبيرة والصغيرة، وتقدير جيد مقبول فقط للجزيئات متوسطة الحجم.

الصيغ

في الصيغ التالية، يرمز G i إلى مساهمة المجموعة. يتم حساب G i لكل مجموعة متاحة. إذا تكررت مجموعة ما عدة مرات، يتم حساب كل تكرار على حدة.

نقطة الغليان العادية

تيب[ك]=198.2+تيب،أنا.{\displaystyle T_{\text{b}}[{\text{K}}]=198.2+\sum T_{{\text{b}},i}.}

نقطة الانصهار

تيم[ك]=122.5+تيم،أنا.{\displaystyle T_{\text{m}}[{\text{K}}]=122.5+\sum T_{{\text{m}},i}.}

درجة الحرارة الحرجة

تيج[ك]=تيب[0.584+0.965تيج،أنا-(تيج،أنا)2]-1.{\displaystyle T_{\text{c}}[{\text{K}}]=T_{\text{b}}\left[0.584+0.965\sum T_{{\text{c}},i}-\left(\sum T_{{\text{c}},i}\right)^{2}\right]^{-1}.}

تتطلب معادلة درجة الحرارة الحرجة هذه نقطة غليان عادية ( T<sub> b</sub>) . يُنصح باستخدام هذه النقطة إذا توفرت قيمة تجريبية. من جهة أخرى، يمكن أيضًا إدخال نقطة الغليان العادية المُقدَّرة باستخدام طريقة جوباك، ولكن هذا سيؤدي إلى زيادة نسبة الخطأ.

ضغط حرج

Pج[حاجِز]=[0.113+0.0032شمالأ-Pج،أنا]-2،{\displaystyle P_{\text{c}}[{\text{bar}}]=\left[0.113+0.0032\,N_{\text{a}}-\sum P_{{\text{c}},i}\right]^{-2},}

حيث أن N a هو عدد الذرات في التركيب الجزيئي (بما في ذلك ذرات الهيدروجين).

المجلد الحرج

Vج[سم3/مول]=17.5+Vج،أنا.{\displaystyle V_{\text{c}}[{\text{cm}}^{3}/{\text{mol}}]=17.5+\sum V_{{\text{c}},i}.}

حرارة التكوين (غاز مثالي، 298 كلفن)

حتشكيل[كيلو جول/مول]=68.29+حاستمارة،أنا.{\displaystyle H_{\text{formation}}[{\text{kJ}}/{\text{mol}}]=68.29+\sum H_{{\text{form}},i}.}

طاقة جيبس ​​للتكوين (غاز مثالي، 298 كلفن)

جيتشكيل[كيلو جول/مول]=53.88+جياستمارة،أنا.{\displaystyle G_{\text{formation}}[{\text{kJ}}/{\text{mol}}]=53.88+\sum G_{{\text{form}},i}.}

السعة الحرارية (الغاز المثالي)

جP[ج/(مولك)]=أأنا-37.93+[بأنا+0.210]تي+[جأنا-3.9110-4]تي2+[دأنا+2.0610-7]تي3.{\displaystyle C_{P}[{\text{J}}/({\text{mol}}\cdot {\text{K}})]=\sum a_{i}-37.93+\left[\sum b_{i}+0.210\right]T+\left[\sum c_{i}-3.91\cdot 10^{-4}\right]T^{2}+\left[\sum d_{i}+2.06\cdot 10^{-7}\right]T^{3}.}

تستخدم طريقة جوباك معادلة متعددة الحدود ذات أربعة معاملات لوصف اعتماد السعة الحرارية للغاز المثالي على درجة الحرارة. وتكون هذه المعاملات صالحة من 273  كلفن إلى حوالي 1000  كلفن، ويمكن توسيع نطاقها إلى 1500 كلفن مع هامش خطأ معين.

حرارة التبخر عند نقطة الغليان العادية

Δحبخار[كيلو جول/مول]=15.30+حبخار،أنا.{\displaystyle \Delta H_{\text{vap}}[{\text{kJ}}/{\text{mol}}]=15.30+\sum H_{{\text{vap}},i}.}

حرارة الانصهار

Δحفوس[كيلو جول/مول]=-0.88+حفوس،أنا.{\displaystyle \Delta H_{\text{fus}}[{\text{kJ}}/{\text{mol}}]=-0.88+\sum H_{{\text{fus}},i}.}

اللزوجة الديناميكية للسائل

ηل[باs]=مwهـxص[(ηأ-597.82)/تي+ηب-11.202]،{\displaystyle \eta _{\text{L}}[{\text{Pa}}\cdot {\text{s}}]=M_{\text{w}}exp{\left[\left(\sum \eta _{a}-597.82\right)/T+\sum \eta _{b}-11.202\right]},}

حيث M w هو الوزن الجزيئي .

تستخدم هذه الطريقة معادلة ذات معلَمين لوصف اعتماد اللزوجة الديناميكية على درجة الحرارة. ويذكر المؤلفون أن هذه المعلمات صالحة من درجة حرارة الانصهار وحتى 0.7 من درجة الحرارة الحرجة ( T <sub>r  </sub> <  0.7).

مساهمات المجموعة

مجموعةتي سيP cV cتي بT mHcaseنموذج Gأبجداندماج الهيدروجينبخار الهيدروجينη aη b
بيانات الحالة الحرجةدرجات حرارة التحولات الطوريةالخواص الحرارية الكيميائيةالسعات الحرارية للغاز المثاليإنثالبيات التحولات الطوريةاللزوجة الديناميكية
مجموعات غير حلقية
−CH 30.0141-0.00126523.58-5.10-76.45-43.961.95E+1-8.08E−31.53E−4−9.67E−80.9082.373548.29-1.719
−CH 20.01890.00005622.8811.27-20.648.42-9.09E−19.50E−2-5.44E-51.19E−82.5902.22694.16-0.199
>CH−0.01640.00204121.7412.6429.8958.36-2.30E+12.04E−1-2.65E−41.20E−70.7491.691-322.151.187
>ج<0.00670.00432718.2546.4382.23116.02-6.62E+14.27E−1-6.41E−43.01E−7-1.4600.636-573.562.307
=CH 20.0113-0.00285618.18-4.32-9.6303.772.36E+1-3.81E−21.72E−4-1.03E-7-0.4731.724495.01-1.539
=CH−0.0129-0.00064624.968.7337.9748.53-8.001.05E−1-9.63E−53.56E−82.6912.20582.28-0.242
=C<0.01170.00113824.1411.1483.9992.36-2.81E+12.08E−1-3.06E−41.46E−73.0632.138غير  متوفرغير  متوفر
=ج=0.00260.00283626.1517.78142.14136.702.74E+1-5.57E−21.01E−4-5.02E−84.7202.661غير  متوفرغير  متوفر
≡CH0.0027-0.0008469.20-11.1879.3077.712.45E+1-2.71E−21.11E−4-6.78E−82.3221.155غير  متوفرغير  متوفر
≡C−0.00200.00163727.3864.32115.51109.827.872.01E−2−8.33E−61.39E-94.1513.302غير  متوفرغير  متوفر
مجموعات الحلقات
−CH 20.01000.00254827.157.75-26.80-3.68-6.038.54E−2-8.00E−6-1.80E−80.4902.398307.53-0.798
>CH−0.01220.00043821.7819.888.6740.99-2.05E+11.62E−1-1.60E−46.24E−83.2431.942-394.291.251
>ج<0.00420.00612721.3260.1579.7287.88-9.09E+15.57E−1-9.00E-44.69E−7-1.3730.644غير  متوفرغير  متوفر
=CH−0.00820.00114126.738.132.0911.30-2.145.74E−2-1.64E−6-1.59E-81.1012.544259.65-0.702
=C<0.01430.00083231.0137.0246.4354.05-8.251.01E−1-1.42E−46.78E−82.3943.059-245.740.912
مجموعات الهالوجين
-F0.0111-0.005727-0.03-15.78-251.92-247.192.65E+1-9.13E−21.91E−4-1.03E-71.398-0.670غير  متوفرغير  متوفر
-Cl0.0105-0.00495838.1313.55-71.55-64.313.33E+1-9.63E−21.87E−4−9.96E−82.5154.532625.45-1.814
-Br0.01330.00577166.8643.43-29.48-38.062.86E+1-6.49E−21.36E−4-7.45E-83.6036.582738.91-2.038
-I0.0068-0.00349793.8441.6921.065.743.21E+1-6.41E−21.26E−4-6.87E−82.7249.520809.55-2.224
مجموعات الأكسجين
-OH (كحول)0.07410.01122892.8844.45-208.04-189.202.57E+1-6.91E−21.77E−4−9.88E−82.40616.8262173.72-5.057
-OH (فينول)0.02400.0184-2576.3482.83-221.65-197.37-2.811.11E−1-1.16E−44.94E−84.49012.4993018.17-7.314
-O- (غير حلقي)0.01680.00151822.4222.23-132.22-105.002.55E+1-6.32E−21.11E−4−5.48E−81.1882.410122.09-0.386
−O− (حلقة)0.00980.00481331.2223.05-138.16-98.221.22E+1-1.26E−26.03E−5−3.86E−85.8794.682440.24-0.953
>C=O (غير حلقي)0.03800.00316276.7561.20-133.22-120.506.456.70E−2-3.57E-52.86E−94.1898.972340.35-0.350
C=O (حلقة)0.02840.00285594.9775.97-164.50-126.273.04E+1-8.29E−22.36E−4-1.31E−70.6.645غير  متوفرغير  متوفر
O=CH− (ألديهيد)0.03790.00308272.2436.90-162.03-143.483.09E+1-3.36E−21.60E−4−9.88E−83.1979.093740.92-1.713
-COOH (حمض)0.07910.007789169.09155.50-426.72-387.872.41E+14.27E−28.04E−5-6.87E−811.05119.5371317.23-2.578
−COO− (إستر)0.04810.00058281.1053.60-337.92-301.952.45E+14.02E−24.02E−5-4.52E-86.9599.633483.88-0.966
=O (بخلاف ما سبق)0.01430.010136-10.502.08-247.61-250.836.821.96E−21.27E−5-1.78E-83.6245.909675.24-1.340
مجموعات النيتروجين
−NH 20.02430.01093873.2366.89-22.0214.072.69E+1-4.12E−21.64E−4−9.76E−83.51510.788غير  متوفرغير  متوفر
>NH (غير حلقي)0.02950.00773550.1752.6653.4789.39-1.217.62E−2-4.86E−51.05E−85.0996.436غير  متوفرغير  متوفر
>NH (حلقة)0.01300.01142952.82101.5131.6575.611.18E+1-2.30E−21.07E−4-6.28E−87.4906.930غير  متوفرغير  متوفر
>N− (غير حلقي)0.01690.0074911.7448.84123.34163.16-3.11E+12.27E−1-3.20E-41.46E−74.7031.896غير  متوفرغير  متوفر
-N = (غير حلقي)0.0255-0.0099غير  متوفر74.60غير  متوفر23.61غير  متوفرغير  متوفرغير  متوفرغير  متوفرغير  متوفرغير  متوفر3.335غير  متوفرغير  متوفر
-N= (حلقة)0.00850.00763457.5568.4055.5279.938.83-3.84E-34.35E−5-2.60E−83.6496.528غير  متوفرغير  متوفر
=NHغير  متوفرغير  متوفرغير  متوفر83.0868.9193.70119.665.69-4.12E−31.28E−4−8.88E−8غير  متوفر12.169غير  متوفرغير  متوفر
−CN0.0496-0.010191125.6659.8988.4389.223.65E+1-7.33E−21.84E−4-1.03E-72.41412.851غير  متوفرغير  متوفر
-NO 20.04370.006491152.54127.24-66.57-16.832.59E+1-3.74E−31.29E−4−8.88E−89.67916.738غير  متوفرغير  متوفر
مجموعات الكبريت
-SH0.00310.00846363.5620.09-17.33-22.993.53E+1-7.58E−21.85E−4-1.03E-72.3606.884غير  متوفرغير  متوفر
-S- (غير حلقي)0.01190.00495468.7834.4041.8733.121.96E+1-5.61E−34.02E−5-2.76E−84.1306.817غير  متوفرغير  متوفر
−S− (حلقة)0.00190.00513852.1079.9339.1027.761.67E+14.81E−32.77E−5-2.11E−81.5575.984غير  متوفرغير  متوفر

مثال على الحساب

الأسيتون (البروبانون) هو أبسط الكيتونات ، ويُفصل إلى ثلاث مجموعات في طريقة جوباك: مجموعتان ميثيل (-CH₃ ) ومجموعة كيتون واحدة (C=O). ولأن مجموعة الميثيل موجودة مرتين، يجب جمع مساهماتها مرتين.

−CH 3>C=O (غير حلقي)
ملكيةعدد المجموعاتقيمة المجموعةعدد المجموعاتقيمة المجموعةجيأنا{\displaystyle \sum G_{i}}القيمة التقديريةوحدة
تي سي
2
0.0141
1
0.0380
0.0662
500.5590
ك
P c
2
-1.20E-03
1
3.10E−03
7.00E−04
48.0250
حاجِز
V c
2
65.0000
1
62.0000
192.0000
209.5000
مل/مول
تي ب
2
23.5800
1
76.7500
123.9100
322.1100
ك
T m
2
-5.1000
1
61.2000
51.0000
173.5000
ك
تكوين الهيدروجين
2
-76.4500
1
-133.2200
-286.1200
-217.8300
كيلوجول/مول
تشكيل G
2
-43.9600
1
-120.5000
-208.4200
-154.5400
كيلوجول/مول
ج ص : أ
2
1.95E+01
1
6.45E+00
4.55E+01
ج ص : ب
2
-8.08E-03
1
6.70E−02
5.08E−02
ج ب : ج
2
1.53E−04
1
-3.57E-05
2.70E−04
C p : d
2
−9.67E−08
1
2.86E−09
-1.91E-07
ج ب
عند درجة حرارة  300  كلفن
75.3264
جول/(مول·كلفن)
اندماج الهيدروجين
2
0.9080
1
4.1890
6.0050
5.1250
كيلوجول/مول
بخار الهيدروجين
2
2.3730
1
8.9720
13.7180
29.0180
كيلوجول/مول
η a
2
548.2900
1
340.3500
1436.9300
η b
2
-1.7190
1
-0.3500
-3.7880
η
عند درجة حرارة  300  كلفن
0.0002942
باس

مراجع

  1. جوباك، كيه جي (1984). منهج موحد لتقدير الخصائص الفيزيائية باستخدام تقنيات إحصائية متعددة المتغيرات (PDF) (MS). معهد ماساتشوستس للتكنولوجيا.
  2. Joback K. G., Reid R. C., "Estimation of Pure-Component Properties from Group-Contributions", Chem. Eng. Commun. , 57, 233–243, 1987.
  3. Lydersen A. L., "Estimation of Critical Properties of Organic Compounds", University of Wisconsin College Engineering, Eng. Exp. Stn. Rep. 3, Madison, Wisconsin, 1955.
  4. Constantinou L., Gani R., "طريقة جديدة لمساهمة المجموعة لتقدير خصائص المركبات النقية"، AIChE J. ، 40(10)، 1697-1710، 1994.
  5. Nannoolal Y., Rarey J., Ramjugernath J., "تقدير خصائص المكونات النقية الجزء 2. تقدير بيانات الخصائص الحرجة عن طريق مساهمة المجموعة"، توازن الطور السائل ، 252(1-2)، 1-27، 2007.
  6. Stein S. E., Brown R. L., "Estimation of Normal Boiling Points from Group Contributions", J. Chem. Inf. Comput. Sci. 34, 581–587 (1994).