تحديد ثوابت الاتزان

يتم تحديد ثوابت الاتزان من أجل تحديد كمية الاتزان الكيميائي . عندما يُعبَّر عن ثابت الاتزان K كنسبة تركيز،

ك=[S]σ[تي]τ[أ]α[ب]β{\displaystyle K={\frac {\mathrm {[S]} ^{\sigma }\mathrm {[T]} ^{\tau }\cdots {\mathrm {[A]} ^{\alpha }\mathrm {[B]} ^{\beta }\cdots }}}

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

إن معادلة التوازن المذكورة أعلاه هي دالة لتركيزات [ A] و[B]، إلخ، للمواد الكيميائية المتوازنة. ويمكن تحديد قيمة ثابت التوازن إذا أمكن قياس أي من هذه التركيزات. وتتمثل الطريقة العامة في قياس التركيز المطلوب لسلسلة من المحاليل ذات التركيزات التحليلية المعروفة للمتفاعلات. عادةً، تُجرى عملية المعايرة بوضع متفاعل واحد أو أكثر في وعاء المعايرة ومتفاعل واحد أو أكثر في السحاحة. بمعرفة التركيزات التحليلية للمتفاعلات في وعاء التفاعل وفي السحاحة، يمكن اشتقاق جميع التركيزات التحليلية كدالة لحجم (أو كتلة) المعاير المضاف.

يمكن اشتقاق ثوابت التوازن عن طريق أفضل مطابقة للبيانات التجريبية مع نموذج كيميائي لنظام التوازن.

الأساليب التجريبية

توجد أربع طرق تجريبية رئيسية. للاطلاع على الطرق الأقل شيوعًا، انظر روسوتي وروسوتي. [ 1 ] في جميع الحالات، يمكن توسيع النطاق باستخدام طريقة التنافس . مثال على تطبيق هذه الطريقة موجود في سيانيد البلاديوم (II) .

القياسات الكهروكيميائية

يُقاس التركيز الحر [A] أو النشاط {A} لنوع كيميائي A باستخدام قطب انتقائي للأيونات ، مثل القطب الزجاجي . إذا تمت معايرة القطب باستخدام معايير النشاط، يُفترض أن معادلة نرنست تنطبق بالشكل التالي:

هـ=هـ0+RتينFln{أ}{\displaystyle E=E^{0}+{\frac {RT}{nF}}\ln \mathrm {\{A\}} }

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

صح=نFRتي(هـ0-هـ){\displaystyle \mathrm {pH} ={\frac {nF}{RT}}\left(E^{0}-E\right)}

عند درجة حرارة 298 كلفن، تساوي وحدة الرقم الهيدروجيني  الواحدة تقريبًا 59 ملي فولت. [ 2 ]  

عند معايرة القطب الكهربائي بمحاليل ذات تركيز معروف، عن طريق معايرة حمض قوي - قاعدة قوية، على سبيل المثال، يتم افتراض معادلة نرنست المعدلة.

هـ=هـ0+sسجل10[أ]{\displaystyle E=E^{0}+s\log _{10}\mathrm {[A]} }

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

النطاق والقيود

القطب الزجاجي هو الأكثر استخدامًا، وهو انتقائي لأيون الهيدروجين. يناسب هذا القطب جميع تفاعلات التوازن الحمضي القاعدي . يمكن قياس قيم log β التي تتراوح بين 2 و11 تقريبًا مباشرةً باستخدام المعايرة الجهدية بواسطة القطب الزجاجي . هذا النطاق الواسع من قيم ثابت الاستقرار (من 100 إلى 10^ 11 تقريبًا ) ممكن بفضل الاستجابة اللوغاريتمية للقطب. تكمن القيود في أن معادلة نرنست لا تنطبق عند قيم الأس الهيدروجيني المنخفضة جدًا أو المرتفعة جدًا.

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

القياسات الطيفية الضوئية

الامتصاص

يفترض أن قانون بير-لامبرت ينطبق.

أ=لεج{\displaystyle A=l\sum {\varepsilon c}}

حيث يُمثل l طول المسار البصري، وε الامتصاصية المولية عند وحدة طول المسار، و c التركيز. قد تُساهم أكثر من مادة في الامتصاصية. من حيث المبدأ، يُمكن قياس الامتصاصية عند طول موجي واحد فقط، ولكن في الممارسة العملية الحالية، من الشائع تسجيل الأطياف الكاملة.

النطاق والقيود

يُذكر عادةً حدٌّ أعلى لقيمة log 10 β يبلغ 4، وهو ما يتوافق مع دقة القياسات، ولكنه يعتمد أيضًا على شدة التأثير. يجب أن تكون أطياف الأنواع المساهمة متميزة بوضوح عن بعضها البعض.

شدة التألق (اللمعان)

يفترض أن شدة الضوء المتناثر هي دالة خطية لتركيزات الأنواع.

أنا=φج{\displaystyle I=\sum \varphi c}

حيث φ هو ثابت التناسب.

النطاق والقيود

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

قياسات الإزاحة الكيميائية بالرنين النووي المغناطيسي

يُفترض أن التبادل الكيميائي سريع على المقياس الزمني للرنين النووي المغناطيسي. يمثل الانزياح الكيميائي الفردي δ المتوسط ​​المرجح بالكسر المولي لانزياحات δ للنوى في الأنواع المساهمة.

دلتا¯=xأنادلتاأناxأنا{\displaystyle {\bar {\delta }}={\frac {\sum x_{i}\delta _{i}}{\sum x_{i}}}}

مثال: تم تحديد قيمة pKa لمجموعة الهيدروكسيل في حمض الستريك من بيانات الإزاحة الكيميائية للكربون -13 لتكون 14.4 . لم يكن من الممكن استخدام قياس الجهد الكهربائي أو مطيافية الأشعة فوق البنفسجية المرئية لهذا التحديد. [ 3 ]

النطاق والقيود

تُفرض دقة قياسات الإزاحة الكيميائية المحدودة حدًا أقصى يبلغ حوالي 4 على لوغاريتم β (10 ) . يقتصر ذلك على الأنظمة غير المغناطيسية. لا يمكن استخدام الرنين النووي المغناطيسي للبروتون (1H NMR) مع محاليل المركبات في  الماء ( 1H₂O ) .

القياسات الحرارية

يُجرى قياس ثابت التوازن (K) وتغير المحتوى الحراري (ΔH ) في آنٍ واحد لمركبات الإضافة بنسبة 1:1 بشكل روتيني باستخدام قياس السعرات الحرارية بالمعايرة متساوية الحرارة . ويُحدّ من إمكانية تطبيق هذه الطريقة على أنظمة أكثر تعقيدًا توفر البرامج المناسبة.

النطاق والقيود

لا تتوفر أدلة كافية في الوقت الحالي.

أسلوب المنافسة

يمكن استخدام طريقة التنافس عندما تكون قيمة ثابت الاستقرار كبيرة جدًا بحيث لا يمكن تحديدها بطريقة مباشرة. وقد استخدمها شوارزنباخ لأول مرة في تحديد ثوابت استقرار معقدات EDTA مع أيونات المعادن.

لتبسيط الأمر، لنفترض تحديد ثابت الاستقراركأب{\displaystyle K_{AB}}من مركب ثنائي، AB ، من كاشف A مع كاشف آخر B.

كأب=[أب][أ][ب]{\displaystyle K_{AB}={\frac {[AB]}{[A][B]}}}

حيث يمثل [X] تركيز النوع X في محلول ذي تركيب معين عند التوازن.

يتم اختيار الرابط C الذي يشكل معقدًا أضعف مع A. ثابت الاستقرار، K<sub> AC </sub>، صغير بما يكفي لتحديده بطريقة مباشرة. على سبيل المثال، في حالة معقدات EDTA، يكون A أيونًا فلزيًا، وقد يكون C بولي أمين مثل ثنائي إيثيلين تريامين .

كأج=[أج][أ][ج]{\displaystyle K_{AC}={\frac {[AC]}{[A][C]}}}

ثابت الاستقرار، K ، لتفاعل المنافسة

أج+بأب+ج{\displaystyle AC+B\leftrightharpoons AB+C}

يمكن التعبير عنها على النحو التالي

ك=[أب][ج][أج][ب]{\displaystyle K={\frac {[AB][C]}{[AC][B]}}}

ويترتب على ذلك أن

كأب=ك×كأج{\displaystyle K_{AB}=K\times K_{AC}}

حيث K هو ثابت الاستقرار لتفاعل المنافسة. وبالتالي، فإن قيمة ثابت الاستقراركأب{\displaystyle K_{AB}}يمكن اشتقاقها من القيم المحددة تجريبياً لـ K وكأج{\displaystyle K_{AC}}.

الأساليب الحسابية

يُفترض أن البيانات التجريبية المُجمّعة تتألف من مجموعة من نقاط البيانات. عند كل نقطة بيانات رقم i ، تُعرف التراكيز التحليلية للمتفاعلات، T A ( i ) ، TB ( i ) وما إلى ذلك، بالإضافة إلى كمية مُقاسة، y i ، تعتمد على واحد أو أكثر من هذه التراكيز التحليلية. يتكون الإجراء الحسابي العام من أربعة مكونات رئيسية :

  1. تعريف نموذج كيميائي للتوازنات
  2. حساب تركيزات جميع الأنواع الكيميائية في كل محلول
  3. تحسين ثوابت التوازن
  4. اختيار الطراز

يمكن حساب قيمة ثابت التوازن لتكوين مركب 1:1، مثل نوع المضيف والضيف، باستخدام تطبيق جدولة مخصص، Bindfit: [ 4 ] في هذه الحالة، يمكن تنفيذ الخطوة 2 بإجراء غير تكراري ويمكن استخدام روتين Solver المبرمج مسبقًا للخطوة 3.

النموذج الكيميائي

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

صأ+qب--أصبq{\displaystyle {\ce {{{\mathit {p}}A}+{\mathit {q}}B\cdots <=>A_{\mathit {p}}B_{\mathit {q}}\cdots }}}:βصq=[أصبq][أ]ص[ب]q{\displaystyle \beta _{pq\cdots }={\frac {[{\ce {A}}_{p}{\ce {B}}_{q}\cdots ]}{[{\ce {A}}]^{p}[{\ce {B}}]^{q}\cdots }}}

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

في المحاليل المائية، يتم تقييد تركيزات البروتون (أيون الهيدرونيوم) وأيون الهيدروكسيد عن طريق التفكك الذاتي للماء.

ح2يا--ح++أوه-{\displaystyle {\ce {H2O <=> H+ + OH-}}}:كدبليو=[ح+][ياح-][ح2يا]{\displaystyle K_{\mathrm {W} }^{'}={\frac {[H^{+}][OH^{-}]}{[H_{2}O]}}}

في المحاليل المخففة، يُفترض أن تركيز الماء ثابت، لذلك يُكتب تعبير التوازن في شكل حاصل الضرب الأيوني للماء.

كدبليو=[ح+][أوه-]{\displaystyle K_{\mathrm {W} }={\ce {[H+]}}[{\ce {OH-}}]\,}

عندما يجب اعتبار كل من H + وOH- كمتفاعلات ، يُستبعد أحدهما من النموذج بتحديد أن تركيزه يُستنتج من تركيز الآخر. عادةً ما يُعطى تركيز أيون الهيدروكسيد بالعلاقة التالية:

[أوه-]=كدبليو[ح+]{\displaystyle [{\ce {OH-}}]={\frac {K_{{\ce {W}}}}{[{\ce {H+}}]}}\,}

في هذه الحالة، يكون ثابت الاتزان لتكوين أيون الهيدروكسيد ذا معاملات قياسية تساوي -1 بالنسبة للبروتون، وصفر بالنسبة للمتفاعلات الأخرى. وهذا له آثار مهمة على جميع تفاعلات اتزان البروتنة في المحلول المائي، وعلى ثوابت التحلل المائي على وجه الخصوص.

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

عادة ما يتم تقدير قيم ثابت التوازن مبدئيًا بالرجوع إلى مصادر البيانات .

حسابات تحديد الأنواع الكيميائية

حساب التوزيع الكيميائي هو عملية يتم فيها حساب تركيزات جميع الأنواع الكيميائية في نظام متوازن، بمعرفة التركيزات التحليلية، T A و T B وما إلى ذلك، للمتفاعلات A و B وما إلى ذلك. وهذا يعني حل مجموعة من المعادلات غير الخطية لتوازن الكتلة.

تيأ=[أ]+1،نكصβصq[أ]ص[ب]qتيب=[ب]+1،نكqβصq[أ]ص[ب]qهـتج.{\displaystyle {\begin{aligned}{\ce {T_{A}}}&=[{\ce {A}}]+\sum _{1,nk}p\beta _{pq\cdots }[{\ce {A}}]^{p}[{\ce {B}}]^{q}\cdots \\{\ce {T_{B}}}&=[{\ce {B}}]+\sum _{1,nk}q\beta _{pq\cdots }[{\ce {A}}]^{p}[{\ce {B}}]^{q}\cdots \\etc.\end{aligned}}}

بالنسبة للتركيزات الحرة [A]، [B]، إلخ. عند قياس الرقم الهيدروجيني (أو القوة الدافعة الكهربائية المكافئة، E)، يتم الحصول على التركيز الحر لأيونات الهيدروجين، [H]، من القيمة المقاسة كما يلي:

[ح]=10-صح{\displaystyle [\mathrm {H} ]=10^{-\mathrm {pH} }}أو[ح]=هـ-نFRتي(هـ-هـ0){\displaystyle [\mathrm {H} ]=e^{\mathrm {{-{\frac {nF}{RT}}}(E-E^{0})} }}

ويتم حساب التركيزات الحرة للمتفاعلات الأخرى فقط. أما تركيزات المركبات فتُستنتج من التركيزات الحرة عبر النموذج الكيميائي.

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

تيأ=0،نكصβصq[أ]ص[ب]qتيب=0،نكqβصq[أ]ص[ب]q{\displaystyle {\begin{aligned}T_{\ce {A}}&=\sum _{0,nk}p\beta _{pq}[{\ce {A}}]^{p}[{\ce {B}}]^{q}\\[4pt]T_{\ce {B}}&=\sum _{0,nk}q\beta _{pq}[{\ce {A}}]^{p}[{\ce {B}}]^{q}\\\end{aligned}}}
β10=β01=1{\displaystyle \beta _{10}=\beta _{01}=1}

وبهذه الطريقة، يتم التعامل مع جميع الأنواع الكيميائية، بما في ذلك المواد المتفاعلة الحرة ، بنفس الطريقة، حيث تم تكوينها من مزيج المواد المتفاعلة المحدد بواسطة المعاملات القياسية.

في نظام المعايرة، تُحسب التراكيز التحليلية للمتفاعلات عند كل نقطة معايرة من الظروف الابتدائية، وتراكيز وحجوم السحاحة. ويُعطى التركيز التحليلي (الكلي) للمتفاعل R عند نقطة المعايرة رقم i بالعلاقة التالية:

تيR=R0+vأنا[R]v0+vأنا{\displaystyle T_{{\ce {R}}}={\frac {{\ce {R}}_{0}+v_{i}{\ce {[R]}}}{v_{0}+v_{i}}}}

حيث R₀ هي الكمية الابتدائية للمادة R في وعاء المعايرة، و v₀ هو الحجم الابتدائي، و[R] هو تركيز المادة R في السحاحة، و vᵢ هو الحجم المضاف. يُعتبر تركيز المادة المتفاعلة غير الموجودة في السحاحة صفرًا.

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

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

تحسين ثابت التوازن

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

رأنا=yأناملاحظة-yأناحساب{\displaystyle r_{i}=y_{i}^{\text{obs}}-y_{i}^{\text{calc}}}

ثم تُعطى دالة الهدف الأكثر عمومية بالصيغة التالية

يو=أناجرأنادبليوأناجرج{\displaystyle U=\sum _{i}\sum _{j}r_{i}W_{ij}r_{j}\,}

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

يو=أنادبليوأناأنارأنا2{\displaystyle U=\sum _{i}W_{ii}r_{i}^{2}}

مع W ij = 0 عندما ji . غالبًا ما يتم استخدام أوزان الوحدة، W ii = 1 ، ولكن في هذه الحالة، تكون القيمة المتوقعة لـ U هي الجذر التربيعي لمتوسط ​​مربعات الأخطاء التجريبية.

يمكن إجراء عملية التصغير باستخدام طريقة جاوس-نيوتن . أولاً، يتم تبسيط دالة الهدف عن طريق تقريبها كتوسيع متسلسلة تايلور من الدرجة الأولى حول مجموعة معلمات أولية، p .

يو=يو0+أنايوصأنادلتاصأنا{\displaystyle U=U^{0}+\sum _{i}{\frac {\partial U}{\partial p_{i}}}\delta p_{i}}

تُضاف الزيادات δ p i إلى المعاملات الأولية المناظرة بحيث تكون U أقل من U 0. عند الحد الأدنى ، تكون المشتقات U / p i مرتبطة ببساطة بعناصر مصفوفة جاكوبي ، J

ججك=yججألجصك{\displaystyle J_{jk}={\frac {\partial y_{j}^{\mathrm {calc} }}{\partial p_{k}}}}

حيث يمثل p k المعامل k في عملية التحسين، وتكون قيمها مساوية للصفر. قد يكون واحد أو أكثر من ثوابت الاتزان من معاملات التحسين. مع ذلك، فإن الكميات المقاسة (انظر أعلاه) المُمثلة بـ y لا تُعبر عنها بدلالة ثوابت الاتزان، بل بدلالة تراكيز الأنواع، وهي دوال ضمنية لهذه المعاملات. لذا، يجب الحصول على عناصر جاكوبيان باستخدام التفاضل الضمني .

يتم حساب الزيادات في المعلمة δ p عن طريق حل المعادلات العادية ، المستمدة من الشروط التي U / p = 0 عند الحد الأدنى.

(جتيدبليوج)دلتاص=جتيدبليور{\displaystyle {\left(J^{\mathrm {T} }WJ\right)\delta p=J^{\mathrm {T} }Wr}}

تُضاف الزيادات δ p بشكل متكرر إلى المعاملات

صن+1=صن+دلتاص{\displaystyle \mathbf {p} ^{n+1}=\mathbf {p} ^{n}+\delta \mathbf {p} }

حيث n هو رقم التكرار. تُعاد حسابات تركيزات الأنواع وقيم y عند كل نقطة بيانات. تستمر التكرارات حتى لا يتحقق انخفاض ملحوظ في U ، أي حتى يتم استيفاء معيار التقارب. مع ذلك، إذا لم تُؤدِّ المعلمات المُحدَّثة إلى انخفاض دالة الهدف، أي إذا حدث تباعد، فيجب تعديل حساب الزيادة. أبسط تعديل هو استخدام جزء، f ، من الزيادة المحسوبة، وهو ما يُعرف بتقطيع الإزاحة.

صن+1=صن+ودلتاص{\displaystyle \mathbf {p} ^{n+1}=\mathbf {p} ^{n}+f\mathbf {\delta p} }

في هذه الحالة، يبقى اتجاه متجه الإزاحة، δ p ، دون تغيير. أما في خوارزمية ليفنبرغ-ماركوارت الأكثر فعالية ، فيتم تدوير متجه الإزاحة باتجاه اتجاه الانحدار الأسرع ، وذلك بتعديل المعادلات العمودية.

(جتيدبليوج+λأنا)دلتاص=جتيدبليور{\displaystyle \mathbf {\left(J^{\mathrm {T} }WJ+\lambda I\right)\delta p=J^{\mathrm {T} }Wr} }

حيث λ هو معامل ماركوارت و I هي مصفوفة الوحدة. وقد تم اقتراح طرق أخرى للتعامل مع التباعد. [ 6 ]

تُثار مشكلة خاصة مع بيانات الرنين النووي المغناطيسي وبيانات قياس الطيف الضوئي. ففي حالة الأخيرة، تكون الكمية المرصودة هي الامتصاصية، A ، ويمكن كتابة قانون بير-لامبرت على النحو التالي:

أλ=ل(εصq..)λجصq..{\displaystyle A_{\lambda }=l\sum (\varepsilon _{pq..})_{\lambda }c_{pq..}}

يمكن ملاحظة أنه بافتراض أن التركيزات، c، معروفة، فإن الامتصاص، A ، عند طول موجي معين،λ{\displaystyle \lambda }وطول المسارل{\displaystyle l}، هي دالة خطية لمعاملات الامتصاص المولية، ε . مع  طول مسار 1 سم، في تدوين المصفوفة

أ=εج{\displaystyle \mathbf {A} ={\boldsymbol {\varepsilon }}\mathbf {C} \,}

هناك طريقتان لحساب معاملات الامتصاص المولية المجهولة.

(1) تُعتبر قيم ε معلماتٍ لعملية التصغير، ويتم بناء مصفوفة جاكوبي على هذا الأساس. ومع ذلك، يتم حساب قيم ε نفسها في كل خطوة من خطوات التحسين باستخدام طريقة المربعات الصغرى الخطية.
ε=(جتيج)-1جتيأ{\displaystyle {\boldsymbol {\varepsilon }}=\mathbf {\left(C^{\mathrm {T} }C\right)^{-1}C^{\mathrm {T} }A} }
باستخدام القيم المُحسَّنة لثوابت التوازن للحصول على التوزيع الكيميائي. المصفوفة
(جتيج)-1جتي{\displaystyle \mathbf {\left(C^{T}C\right)^{-1}C^{T}} }
هذا مثال على المعكوس الزائف .
أظهر غولوب وبيريرا [ 7 ] كيف يمكن اشتقاق المعكوس الزائف بحيث يمكن حساب الزيادات في المعلمات لكل من الامتصاصية المولية وثوابت التوازن عن طريق حل المعادلات العادية.
(2) يُكتب قانون بير-لامبرت على النحو التالي
ελ=أλ-1ج{\displaystyle \mathbf {\boldsymbol {\varepsilon }} _{\lambda }=\mathbf {A} _{\lambda }^{-1}\mathbf {C} \,}
يتم إيجاد قيم الامتصاص المولي المجهولة لجميع الأنواع "الملونة" باستخدام طريقة المربعات الصغرى الخطية غير التكرارية، وذلك بقياس طول موجي واحد في كل مرة. تُجرى الحسابات مرة واحدة في كل دورة تحسين، باستخدام قيم ثابت الاستقرار التي تم الحصول عليها في تلك الدورة لحساب قيم تركيز الأنواع في المصفوفة.ج{\displaystyle \mathbf {C} }.

أخطاء المعلمات والارتباط

في المنطقة القريبة من الحد الأدنى لدالة الهدف، U ، يقترب النظام من نظام المربعات الصغرى الخطي، والذي

ص=(جتيدبليوج)-1جتيدبليوyoبs{\displaystyle \mathbf {p=\left(J^{\mathrm {T} }WJ\right)^{-1}J^{\mathrm {T} }Wy^{\mathrm {obs} }} }

لذا، فإن قيم المعلمات هي (تقريبًا) توليفات خطية من قيم البيانات المرصودة ، ويمكن الحصول على أخطاء المعلمات، p ، من خلال انتشار الخطأ من المشاهدات، y obs ، باستخدام الصيغة الخطية. لنفترض أن مصفوفة التباين والتباين المشترك للمشاهدات يُرمز لها بـ Σ ولمصفوفة المعلمات بـ Σ p . عندئذٍ،

Σص=(جتيدبليوج)-1جتيدبليوΣyدبليوتيج(جتيدبليوج)-1{\displaystyle \mathbf {\Sigma ^{p}=\left(J^{\mathrm {T} }WJ\right)^{-1}J^{\mathrm {T} }W\Sigma ^{y}W^{\mathrm {T} }J(J^{\mathrm {T} }WJ)^{-1}} }

عندما يكون W = ( Σ y ) − 1 ، فإن هذا يتبسط إلى

Σص=(جتيدبليوج)-1{\displaystyle \mathbf {\Sigma ^{p}=\left(J^{\mathrm {T} }WJ\right)^{-1}} }

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

دبليوك=1σهـ2+(هـv)ك2σv2{\displaystyle W_{k}={\frac {1}{\sigma _{E}^{2}+\left({\frac {\partial E}{\partial v}}\right)_{k}^{2}\sigma _{v}^{2}}}}

حيث σ E هو الخطأ في جهد القطب أو الرقم الهيدروجيني ، ( E / v)k هو ميلمنحنى المعايرةوσvهو الخطأ في الحجم المضاف.

عند استخدام أوزان الوحدة ( W = I ، p = ( JTJ ) ⁻¹JTy ) ، يُفترض أن الأخطاء التجريبية غير مترابطة ومتساوية: Σy = σ²I ، حيث σ² هو تباين مشاهدة ذات وزن وحدة، و I هي مصفوفة الوحدة . في هذه الحالة ، يُقارب σ² بـ

σ2=يوند-نص{\displaystyle \sigma ^{2}={\frac {U}{n_{\mathrm {d} }-n_{\mathrm {p} }}}}

حيث U هي القيمة الدنيا لدالة الهدف و n d و n p هما عدد البيانات والمعلمات على التوالي.

Σص=يوند-نص(جتيج)-1{\displaystyle \mathbf {\Sigma ^{p}} ={\frac {U}{n_{\mathrm {d} }-n_{\mathrm {p} }}}\left(\mathbf {J} ^{\mathrm {T} }\mathbf {J} \right)^{-1}}

في جميع الحالات، يُعطى تباين المعلمة pᵢ بالصيغة Σ pᵢᵢ ، ويُعطى التغاير بين المعلمتين pᵢ و pⱼ بالصيغة Σ pⱼᵢ . الانحراف المعياري هو الجذر التربيعي للتباين. تعكس تقديرات الخطأ هذه الأخطاء العشوائية فقط في القياسات . أما عدم اليقين الحقيقي في المعلمات فهو أكبر بسبب وجود أخطاء منهجية ، والتي لا يمكن تحديدها كميًا بحكم تعريفها.

لاحظ أنه على الرغم من أن الملاحظات قد تكون غير مترابطة، إلا أن المعلمات تكون مترابطة دائمًا .

الثوابت المشتقة

بعد تحسين الثوابت التراكمية، يُفيد غالبًا اشتقاق ثوابت متدرجة منها. وتتمثل الطريقة العامة في كتابة الصيغ التعريفية لجميع الثوابت المعنية، ثم مساواة التراكيز. على سبيل المثال، لنفترض أننا نرغب في اشتقاق قيمة pKa لإزالة بروتون واحد من حمض ثلاثي القاعدة، LH₃ ، مثل حمض الستريك .

ل3-+ح+-- LH2-: [LH2-]=β11[ل3-][ح+]ل3-+2ح+-- LH2-: [LH2-]=β12[ل3-][ح+]2ل3-+3ح+-- LH3: [LH3]=β13[ل3-][ح+]3{\displaystyle {\begin{aligned}{\ce {L^3-}}+{\ce {H+ <=>}}\ {\ce {LH^2-}}&:\ [{\ce {LH^2-}}]=\beta _{11}[{\ce {L^3-}}][{\ce {H+}}]\\{\ce {L^3-}}+{\ce {2H+ <=>}}\ {\ce {LH2^-}}&:\ [{\ce {LH2^-}}]=\beta _{12}[{\ce {L^3-}}][{\ce {H+}}]^{2}\\{\ce {L^3-}}+{\ce {3H+ <=>}}\ {\ce {LH3}}&:\ [{\ce {LH3}}]=\beta _{13}[{\ce {L^3-}}][{\ce {H+}}]^{3}\end{aligned}}}

ثابت الارتباط التدريجي لتكوين LH 3 يُعطى بالعلاقة التالية

LH2-+ح+--LH3 ؛ [LH3]=ك[LH2-][ح+]{\displaystyle {\ce {{LH2^{-}}+H+<=>LH3\ ;\quad \ [LH3]}}=K[{\ce {LH2^{-}}}][{\ce {H+}}]}

استبدل تعبيرات تركيزي LH 3 و LH 2 في هذه المعادلة

β13[ل3-][ح+]3=كβ12[ل3-][ح+]2[ح+]{\displaystyle \beta _{13}[{\ce {L^3-}}][{\ce {H+}}]^{3}=K\beta _{12}[{\ce {L^3-}}][{\ce {H+}}]^{2}[{\ce {H+}}]}

ومن ثم

β13=كβ12؛ك=β13β12{\displaystyle \beta _{13}=K\beta _{12};K={\frac {\beta _{13}}{\beta _{12}}}\,}

وبما أن pKa = −log 10 1 / K ⁠ فإن قيمتها تُعطى بواسطة

صكأ1=سجل10β13-سجل10β12{\displaystyle {\ce {p}}K_{{\ce {a1}}}=\log _{10}\beta _{13}-\log _{10}\beta _{12}\,}
صكأ2=سجل10β12-سجل10β11{\displaystyle {\ce {p}}K_{{\ce {a2}}}=\log _{10}\beta _{12}-\log _{10}\beta _{11}\,}
صكأ3=سجل10β11{\displaystyle {\ce {p}}K_{{\ce {a3}}}=\log _{10}\beta _{11}\,}

لاحظ الترقيم العكسي لـ pK و log β. عند حساب الخطأ في الثابت التدريجي، يجب مراعاة حقيقة أن الثوابت التراكمية مترابطة . عن طريق انتشار الخطأ

σك2=σβ122+σβ132-2σβ12σβ13ρ12،13{\displaystyle \sigma _{K}^{2}=\sigma _{\beta _{12}}^{2}+\sigma _{\beta _{13}}^{2}-2\sigma _{\beta _{12}}\sigma _{\beta _{13}}\rho _{12,13}\,}

و

σسجل10ك=σكك{\displaystyle \sigma _{\log _{10}K}={\frac {\sigma _{K}}{K}}}

اختيار الطراز

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

دالة الهدف

عندما يتم اشتقاق الأوزان بشكل صحيح من تقديرات الخطأ التجريبي، فإن القيمة المتوقعة لـ U / n dn p هي 1. [ 8 ] لذلك من المفيد جدًا تقدير الأخطاء التجريبية واشتقاق بعض الأوزان المعقولة منها لأن هذا مؤشر مطلق على جودة المطابقة.

عند استخدام الأوزان الموحدة ، يُفترض أن جميع المشاهدات لها نفس التباين. ومن المتوقع أن يكون U / n dn p مساوياً لهذا التباين.

أخطاء المعلمات

يُفضّل أن تكون الأخطاء في ثوابت الاستقرار متناسبة تقريبًا مع الخطأ التجريبي. على سبيل المثال، في بيانات معايرة الرقم الهيدروجيني، إذا تم قياس الرقم الهيدروجيني بدقة منزلتين عشريتين، فلا ينبغي أن تتجاوز أخطاء لوغاريتم β (log 10 β ) 0.01. في الدراسات الاستكشافية التي لا تُعرف فيها طبيعة الأنواع الموجودة مسبقًا، يمكن اختبار ومقارنة العديد من النماذج الكيميائية المختلفة. ستكون هناك نماذج قد تكون فيها الشكوك في أفضل تقدير لثابت التوازن أكبر قليلًا أو حتى بشكل ملحوظ من σ pH ، خاصةً مع تلك الثوابت التي تحكم تكوين الأنواع الثانوية نسبيًا، ولكن يبقى القرار بشأن حجم الشكوك المقبول أمرًا شخصيًا. يمكن جعل عملية اتخاذ القرار بشأن تضمين التوازنات غير المؤكدة نسبيًا في النموذج، ومقارنة النماذج المتنافسة بشكل عام، موضوعية، وقد أوضحها هاميلتون. [ 8 ]

توزيع البواقي

عند الحد الأدنى في يمكن تقريب النظام إلى نظام خطي، وترتبط البواقي في حالة الأوزان الموحدة بالملاحظات بواسطة

ر=yoبs-ج(جتيتي)-1جتيyoبs{\displaystyle \mathbf {r=y^{\mathrm {obs} }-J\left(J^{\mathrm {T} }T\right)^{-1}J^{\mathrm {T} }y^{\mathrm {obs} }} }

تُعرف المصفوفة المتناظرة والمتساوية القوة J ( J T T ) −1 J في أدبيات الإحصاء باسم مصفوفة القبعة ، H. وبالتالي،

ر=(أنا-ح)yoبs{\displaystyle \mathbf {r=\left(I-H\right)y^{\mathrm {obs} }} }

و

مر=(أنا-ح)مy(أنا-ح){\displaystyle \mathbf {M^{r}=\left(I-H\right)M^{y}\left(I-H\right)} }

حيث I هي مصفوفة الوحدة ، و Mr و My هما مصفوفتا التباين والتباين المشترك للبواقي والمشاهدات على التوالي. وهذا يُظهر أنه على الرغم من أن المشاهدات قد تكون غير مترابطة، فإن البواقي تكون مترابطة دائمًا .

يُظهر الرسم البياني على اليمين نتيجة تحسين ثوابت استقرار Ni (Gly) و Ni(Gly) و Ni(Gly) ⁻³ (حيث GlyH = جليسين ). تُعرض القيم المرصودة على شكل معينات زرقاء، وتُضاف إليها تراكيز الأنواع كنسبة مئوية من إجمالي النيكل. تظهر القيم المتبقية في المربع السفلي. لا تتوزع هذه القيم المتبقية عشوائيًا كما هو متوقع، ويعود ذلك إلى تباين جهود وصلات السائل وتأثيرات أخرى عند أسطح التماس بين الزجاج والسائل. هذه التأثيرات بطيئة جدًا مقارنةً بمعدل الوصول إلى حالة التوازن.

القيود المادية

تُدمج بعض القيود الفيزيائية عادةً في الحسابات. على سبيل المثال، يجب أن تكون جميع تراكيز المواد المتفاعلة الحرة والأنواع الكيميائية ذات قيم موجبة، ويجب أن تكون ثوابت الارتباط ذات قيم موجبة أيضاً.

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

القيود الكيميائية

عند تحديد ثوابت استقرار معقدات الفلزات والروابط، من الممارسات الشائعة تثبيت ثوابت بروتنة الروابط عند قيم تم تحديدها باستخدام بيانات مستقاة من محاليل خالية من الفلزات. كما تُثبّت ثوابت تحلل أيونات الفلزات عادةً عند قيم تم الحصول عليها باستخدام محاليل خالية من الروابط. عند تحديد ثوابت استقرار المعقدات الثلاثية، M <sub>p</sub> A<sub> q</sub> B <sub> r </sub>، من الممارسات الشائعة تثبيت قيم المعقدات الثنائية المقابلة M <sub>p′ </sub> A<sub> q′</sub> و M <sub>p′′</sub> B <sub>q′′</sub> ، عند قيم تم تحديدها في تجارب منفصلة. يقلل استخدام هذه القيود من عدد المعاملات المطلوب تحديدها، ولكنه قد يؤدي إلى التقليل من تقدير الأخطاء المحسوبة في قيم ثوابت الاستقرار المُحسّنة.

نماذج أخرى

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

المبادئ الديناميكية الحرارية لتفاعلات المضيف والضيف

يمكن تقييم الديناميكا الحرارية لتفاعل المضيف والضيف باستخدام مطيافية الرنين النووي المغناطيسي، ومطيافية الأشعة فوق البنفسجية/المرئية، ومعايرة السعرات الحرارية متساوية الحرارة. [ 9 ] يوفر التحليل الكمي لقيم ثابت الارتباط معلومات ديناميكية حرارية مفيدة.

ثابت الارتباط ،كأ{\displaystyle K_{a}^{\ominus }}يمكن تعريفها من خلال التعبير

كأ={حجي}{ح}{جي}=[حجي][ح][جي]×Γ{\displaystyle K_{a}^{\ominus }={\frac {\{HG\}}{\{H\}\{G\}}}={\frac {[HG]}{[H][G]}}\times \Gamma }

حيث يمثل {HG} النشاط الديناميكي الحراري للمركب عند الاتزان. يمثل {H} نشاط المضيف، ويمثل {G} نشاط الضيف. الكميات[حجي]{\displaystyle [HG]}،[ح]{\displaystyle [H]}و[جي]{\displaystyle [G]}وهي التراكيز المقابلة وΓ{\displaystyle \Gamma }هو ناتج قسمة معاملات النشاط .

في الممارسة العملية، يتم تعريف ثابت التوازن عادةً بدلالة التركيزات.

كأ=[حجي][ح][جي]{\displaystyle K_{a}={\frac {[HG]}{[H][G]}}}

عند استخدام هذا التعريف، يُفترض ضمنيًا أن حاصل قسمة معاملات النشاط له قيمة عددية تساوي واحدًا. عندئذٍ يبدو أن ثابت التوازن،كأ{\displaystyle K_{A}}له بُعد 1/التركيز، لكن هذا لا يمكن أن يكون صحيحًا نظرًا لتغير طاقة جيبس ​​الحرة القياسية.Δجي{\displaystyle \Delta G^{\ominus }}يتناسب مع لوغاريتمكأ{\displaystyle K_{A}}.

Δجي=-Rتيlnكأ{\displaystyle \Delta G^{\ominus }=-RT\ln {K_{A}^{\ominus }}}

يتم حل هذه المفارقة الظاهرية عندما يكون بُعدΓ{\displaystyle \Gamma }يُعرَّف بأنه مقلوب بُعد حاصل قسمة التركيزات. ويترتب على ذلك أنΓ{\displaystyle \Gamma }يُعتبر ثابتًا في جميع الظروف التجريبية ذات الصلة. ومع ذلك، من الممارسات الشائعة إضفاء بُعد، مثل ملي مول لكل لتر أو ميكرو مول لكل لتر، على قيمة K التي تم تحديدها تجريبيًا.

كبيركأ{\displaystyle K_{a}}تشير القيمة إلى أن جزيئات المضيف والضيف تتفاعل بقوة لتشكيل مركب المضيف والضيف.

تحديد قيم ثابت الارتباط وثابت الحركة

تكوين معقدات بسيطة بين المضيف والضيف

عندما تتحد جزيئات المضيف والضيف لتشكيل مركب واحد، يتم تمثيل حالة التوازن على النحو التالي:

ح+جيحجي{\displaystyle H+G\leftrightharpoons HG}

ويُعرَّف ثابت الاتزان، K، على النحو التالي:

ك=[حجي][ح][جي]{\displaystyle K={\frac {[HG]}{[H][G]}}}

حيث يرمز [X] إلى تركيز النوع الكيميائي X (يُفترض أن جميع معاملات النشاط لها قيمة عددية تساوي 1). معادلات توازن الكتلة، عند أي نقطة بيانات،

تيح=[ح]+ك[ح][جي]{\displaystyle T_{H}=[H]+K[H][G]}
تيجي=[جي]+ك[ح][جي]{\displaystyle T_{G}=[G]+K[H][G]}

أينتيجي{\displaystyle T_{G}}وتيح{\displaystyle T_{H}}يمكن اختزال التركيزات الكلية للمضيف والضيف إلى معادلة تربيعية واحدة، على سبيل المثال، [G]، وبالتالي يمكن حلها تحليليًا لأي قيمة معينة لـ K. ومن ثم يمكن اشتقاق التركيزات [H] و[HG].

[ح]=تيح-تيجي+[جي]{\displaystyle [H]=T_{H}-T_{G}+[G]}
[حجي]=ك[ح][جي]{\displaystyle [HG]=K[H][G]}

الخطوة التالية في الحساب هي حساب القيمة،Xأناجألج{\displaystyle X_{i}^{calc}}، بكمية تتوافق مع الكمية المرصودةXأناoبs{\displaystyle X_{i}^{obs}}ثم، يمكن تعريف مجموع المربعات، U، على جميع نقاط البيانات، np، على النحو التالي:

يو=أنا=1،نص(Xأناoبs-Xأناجألج)2{\displaystyle U=\sum _{i=1,np}(X_{i}^{obs}-X_{i}^{calc})^{2}}

ويمكن تقليل هذا التأثير بالنسبة لقيمة ثابت الاستقرار K، ومعامل مثل الإزاحة الكيميائية للنوع HG (بيانات الرنين النووي المغناطيسي) أو امتصاصه المولي (بيانات الأشعة فوق البنفسجية/المرئية). ويمكن إجراء عملية التقليل في برنامج جداول بيانات مثل إكسل باستخدام أداة SOLVER المدمجة.

ينطبق هذا الإجراء على نواتج الإضافة بنسبة 1:1.

تفاعل التكوين المعقد العام

لكل حالة توازن تتضمن مضيفًا، H، وضيفًا، G

صح+qجيحصجيq{\displaystyle pH+qG\leftrightharpoons H_{p}G_{q}}

ثابت الاتزان،βصq{\displaystyle \beta _{pq}}، ويُعرَّف على النحو التالي

βصq=[حصجيq][ح]ص[جي]q{\displaystyle \beta _{pq}={\frac {[H_{p}G_{q}]}{[H]^{p}[G]^{q}}}}

قيم التركيزات الحرة،[ح]{\displaystyle [H]}و[جي]{\displaystyle [G]}يتم الحصول عليها عن طريق حل معادلات توازن الكتلة بقيم معروفة أو مقدرة لثوابت الاستقرار.

تيح=[ح]+صβصq[ح]ص[جي]q{\displaystyle T_{H}=[H]+\sum p\beta _{pq}[H]^{p}[G]^{q}}
تيجي=[جي]+qβصq[ح]ص[جي]q{\displaystyle T_{G}=[G]+\sum q\beta _{pq}[H]^{p}[G]^{q}}

بعد ذلك، يمكن أيضًا حساب تركيزات كل نوع من الأنواع المعقدة على النحو التالي:[حصجيq]=βصq[ح]ص[جي]q{\displaystyle [H_{p}G_{q}]=\beta _{pq}[H]^{p}[G]^{q}}تختلف العلاقة بين تركيز نوع معين والكمية المقاسة باختلاف تقنية القياس، كما هو موضح في كل قسم أعلاه. وباستخدام هذه العلاقة، يمكن تحسين مجموعة المعاملات، وقيم ثابت الاستقرار، وقيم خصائص مثل معامل الامتصاص المولي أو الإزاحات الكيميائية المحددة، من خلال عملية تحسين المربعات الصغرى غير الخطية . لمزيد من التفاصيل حول النظرية، انظر قسم "تحديد ثوابت الاتزان". تجدون بعض البرامج الحاسوبية المتخصصة في قسم " التطبيقات" .

التعاون

في التعاونية، يؤثر ارتباط الليجاند الأولي على ألفة الخلية المضيفة لليجاندات اللاحقة. في التعاونية الإيجابية، يعزز حدث الارتباط الأول ألفة الخلية المضيفة لليجاندات الأخرى. ومن أمثلة التعاونية الإيجابية والسلبية الهيموجلوبين ومستقبل الأسبارتات، على التوالي. [ 10 ]

الارتباط العام بين المضيف والضيف. (1) ارتباط الضيف أ. (2) ارتباط الضيف ب. (3) ارتباط الضيف أ-ب بتعاون إيجابي. (4) ارتباط الضيف أ-ب بتعاون سلبي.

تمت دراسة الخصائص الديناميكية الحرارية للتعاونية بهدف تحديد المعايير الرياضية التي تميز التعاونية الإيجابية عن التعاونية السلبية. وتنص معادلة طاقة جيبس ​​الحرة التقليدية على ما يلي:Δجي=Δح-تيΔS {\displaystyle \Delta G=\Delta H-T\Delta S\ }مع ذلك، لتحديد كمية التعاونية في نظام المضيف والضيف، يجب مراعاة طاقة الربط. يوضح الرسم التخطيطي على اليمين ربط A، وربط B، والربط التعاوني الإيجابي بين A وB، وأخيرًا، الربط التعاوني السلبي بين A وB. لذلك، يمكن صياغة معادلة طاقة جيبس ​​الحرة بصيغة بديلة كما يلي:

ΔجيS=Δجيأ+Δجيب-Δجيأب{\displaystyle \Delta G_{S}^{\circ }=\Delta G_{A}^{\circ }+\Delta G_{B}^{\circ }-\Delta G_{AB}^{\circ }}
ΔحS=Δحأ+Δحب-Δحأب{\displaystyle \Delta H_{S}^{\circ }=\Delta H_{A}^{\circ }+\Delta H_{B}^{\circ }-\Delta H_{AB}^{\circ }}
 تيΔجيS=تيΔحأ+تيΔحب-تيΔSأب{\displaystyle \ T\Delta G_{S}^{\circ }=T\Delta H_{A}^{\circ }+T\Delta H_{B}^{\circ }-T\Delta S_{AB}^{\circ }}

أين:

Δجيأ{\displaystyle \Delta G_{A}^{\circ }}= الطاقة الحرة للربط أ
Δجيب{\displaystyle \Delta G_{B}^{\circ }}= الطاقة الحرة للربط B
ΔجيS{\displaystyle \Delta G_{S}^{\circ }}= الطاقة الحرة للربط بين A و B المرتبطين
Δجيأب{\displaystyle \Delta G_{AB}^{\circ }}= مجموع طاقات الربط الحرة

يُعتبر أنه إذاΔجيS{\displaystyle \Delta G_{S}^{\circ }}أكثر من مجموعΔجيأ{\displaystyle \Delta G_{A}^{\circ }}وΔجيب{\displaystyle \Delta G_{B}^{\circ }}إنه متعاون بشكل إيجابي. إذاΔجيS{\displaystyle \Delta G_{S}^{\circ }}كلما كانت قيمة أقل، كان التفاعل تعاونيًا سلبيًا. [ 11 ] لا تقتصر كيمياء المضيف والضيف على تفاعلات المستقبلات والروابط، بل تتجلى أيضًا في أنظمة الاقتران الأيوني. تُدرس هذه التفاعلات في وسط مائي باستخدام مضيفات عضوية فلزية اصطناعية وجزيئات ضيف عضوية. على سبيل المثال، يرتبط مستقبل متعدد الكاتيونات يحتوي على النحاس (المضيف) بجزيئات مثل رباعي الكربوكسيلات، وثلاثي الكربوكسيلات، والأسبارتات، والأسيتات (الضيوف). توضح هذه الدراسة أن الإنتروبيا ، وليس الإنثالبي، هي التي تحدد طاقة الربط للنظام، مما يؤدي إلى تعاونية سلبية. ينشأ التغير الكبير في الإنتروبيا من إزاحة جزيئات المذيب المحيطة بالرابط والمستقبل. عندما ترتبط عدة جزيئات أسيتات بالمستقبل، فإنها تُطلق جزيئات ماء أكثر إلى البيئة المحيطة مقارنةً بجزيء رباعي الكربوكسيلات. أدى ذلك إلى انخفاض في الطاقة الحرة، مما يشير إلى أن النظام يتعاون بشكل سلبي. [ 12 ] في دراسة مماثلة، باستخدام غوانيدينيوم و Cu(II) وضيوف متعدد الكربوكسيلات، تبين أن التفاعل التعاوني الإيجابي يتحدد إلى حد كبير بالمحتوى الحراري. [ 13 ] بالإضافة إلى الدراسات الديناميكية الحرارية، فإن كيمياء المضيف والضيف لها تطبيقات بيولوجية أيضًا.

التطبيقات

بعض الأنظمة البسيطة قابلة للحساب باستخدام جداول البيانات. [ 4 ] [ 14 ]

نُشر عدد كبير من برامج الحاسوب العامة لحساب ثوابت الاتزان. انظر [ 15 ] للاطلاع على قائمة المراجع. البرامج الأكثر استخدامًا هي:

مراجع

  1. روسوتي، ف.  ج.  س.؛ روسوتي، هـ. (1961). تحديد ثوابت الاستقرار . ماكجرو هيل.
  2. "تعريفات مقاييس الرقم الهيدروجيني، والقيم المرجعية القياسية، وقياس الرقم الهيدروجيني، والمصطلحات ذات الصلة" (ملف PDF) . مجلة الكيمياء البحتة والتطبيقية. 57 : 531-542 . 1985. doi : 10.1351/pac198557030531 . S2CID 14182410 . 
  3. سيلفا، أندريه إم إن؛ كونغ، شياولي؛ هايدر، روبرت سي. (2009). "تحديد قيمة pKa لمجموعة الهيدروكسيل في مركبات ألفا-هيدروكسي كاربوكسيلات، السترات والمالات واللاكتات، باستخدام الرنين المغناطيسي النووي للكربون -13 : دلالات على تنسيق المعادن في الأنظمة البيولوجية". المعادن الحيوية . 22 (5): 771-778 . doi : 10.1007/s10534-009-9224-5 . PMID 19288211. S2CID 11615864 .  
  4. 1 2 هيبرت، د.ب.؛ ثوردارسون، ب. (2017). "موت مخطط أيوب، والشفافية، والعلوم المفتوحة، والأدوات الإلكترونية، وأساليب تقدير عدم اليقين، وتطورات أخرى في تحليل بيانات الكيمياء فوق الجزيئية" . مجلة الاتصالات الكيميائية . 52 (87): 12792-12805 . doi : 10.1039/c6cc03888c . PMID 27779264 . 
  5. موتيكايتيس، ر. ج.؛ مارتيل، أ. إ. (1982). "بيست - برنامج جديد لحساب دقيق لمعاملات التوازن للأنظمة المعقدة متعددة المكونات". المجلة الكندية للكيمياء 60 ( 19): 2403-2409 . doi : 10.1139/v82-347 .  
  6. 1 2 بوتفين، ب. ج. (1990). "نمذجة توازن المحاليل المعقدة. الجزء الأول: تحسين سريع وخالٍ من المتاعب لثوابت التوازن باستخدام طريقة المربعات الصغرى" . المجلة الكندية للكيمياء . 68 (12): 2198-2207 . doi : 10.1139/v90-337 . 
  7. غولوب، ج. هـ.؛ بيريرا، ف. (1973). "تفاضل الدوال المعكوسة الزائفة ومسائل المربعات الصغرى غير الخطية التي تنفصل متغيراتها". مجلة SIAM للتحليل العددي 10 ( 2): 413-432 . Bibcode : 1973SJNA...10..413G . doi : 10.1137/0710036 . 
  8. 1 2 هاميلتون، دبليو سي (1964). الإحصاء في العلوم الفيزيائية . نيويورك، نيويورك: رونالد برس. 
  9. بينيرو، أ.؛ بانكي، إكس.؛ بيريز-كاساس، إس.؛ توفار، إي.؛ غارسيا، أ.؛ فيلا، أ.؛ أميغو، أ.؛ مارك، إيه إي؛ كوستا، إم. (2007). "حول توصيف معقدات المضيف والضيف: التوتر السطحي، والقياس الحراري، والديناميكا الجزيئية للسايكلودكسترينات مع عامل سطحي غير أيوني". مجلة الكيمياء الفيزيائية ب . 111 (17): 4383-92 . doi : 10.1021/jp0688815 . PMID 17428087 . 
  10. كوشلاند، د. (1996). "الأساس البنيوي للتعاونية السلبية: المستقبلات والإنزيمات". الرأي الحالي في البيولوجيا البنيوية . 6 (6): 757-761 . doi : 10.1016/S0959-440X(96)80004-2 . PMID 8994875 . 
  11. جينكس، دبليو بي (1981). "حول إسناد وجمع طاقات الربط" . وقائع الأكاديمية الوطنية للعلوم، الولايات المتحدة الأمريكية . 78 (7): 4046-4050 . Bibcode : 1981PNAS...78.4046J . doi : 10.1073 / pnas.78.7.4046 . PMC 319722. PMID 16593049 .  
  12. دوبرزانسكا، ل؛ لويد، ج؛ إسترهويزن، س؛ باربور، ل (2003). "دراسات حول الأصل الديناميكي الحراري للتعاونية السلبية في التعرف الجزيئي عن طريق اقتران الأيونات". مجلة الجمعية الكيميائية الأمريكية . 125 (36): 10963-10970 . doi : 10.1021/ja030265o . PMID 12952478 . 
  13. هيوز، أ.؛ أنسلين، إي. (2007). "مضيف كاتيون يُظهر تعاونًا إيجابيًا في الماء" . وقائع الأكاديمية الوطنية للعلوم، الولايات المتحدة الأمريكية . 104 (16): 6538-6543 . Bibcode : 2007PNAS..104.6538H . doi : 10.1073 /pnas.0609144104 . PMC 1871821. PMID 17420472 .  
  14. بيلو، إي. جوزيف (2011). إكسل للكيميائيين: دليل شامل ( الطبعة الثالثة). وايلي-في سي إتش. رقم ISBN  978-0-470-38123-6.
  15. غانز، ب.؛ ساباتيني، أ.؛ فاكا، أ. (1996). "دراسة التوازنات في المحلول. تحديد ثوابت التوازن باستخدام مجموعة برامج HYPERQUAD". تالانتا . 43 (10): 1739-1753 . doi : 10.1016/0039-9140(96)01958-3 . PMID 18966661 . 
  16. مارتيل، أ. إي.؛ موتيكايتيس، ر. ج. (1992). تحديد واستخدام ثوابت الاستقرار . وايلي-في سي إتش. رقم ISBN   0471188174.
  17. 1 2 ليجيت، د. ج.، محرر. (1985). الطرق الحسابية لتحديد ثوابت التكوين . مطبعة بلينوم. ISBN  978-0-306-41957-7.
  18. غامب، هـ.؛ مايدر، م.؛ ماير، س. ج.؛ زوبربوهلر، أ. (1985). "حساب ثوابت التوازن من بيانات التحليل الطيفي متعدد الأطوال الموجية - اعتبارات رياضية". تالانتا . 32 (95): 95-101 . doi : 10.1016/0039-9140(85)80035-7 . PMID 18963802 .