عندما تغمس حلقة في الماء والصابون وتنفخ فيها بلطف، يتمدد غشاء سائل رقيق وينفصل ليسترخي في غضون أجزاء من الألف من الثانية مكوِّناً كرة متلألئة بألوان قوس قزح تسبح في الهواء. وإذا تأملت رغوة الصابون عبر عدسة مكبرة، ستكتشف أن ما يبدو فوضى عشوائية يخضع في الحقيقة لنظام هندسي صارم لا يقبل المساومة: فجدران الأغشية السائلة لا تلتقي أبداً أربعةً أو خمسةً معاً، بل تلتقي حصرياً ثلاثةً ثلاثةً على طول قنوات سائلة تُسمى حدود بلاتو بزاوية دقيقة قدرها $120^\circ$ ($2\pi/3$)، وتلتقي هذه القنوات بدورها حصرياً أربعاً أربعاً عند الرؤوس وفق الزاوية الرباعية الأوجه الدقيقة $\arccos(-1/3) \approx 109.4712^\circ$.

لماذا تتخذ فقاعة الصابون الحرة شكلاً كروياً دائماً؟ ولماذا عندما نصل فقاعة صغيرة بأخرى كبيرة عبر أنبوب تتقلص الفقاعة الصغيرة وتنفخ الفقاعة الكبيرة على عكس البديهة — ولماذا لا يؤدي هذا عدم الاستقرار نفسه إلى انهيار 300 مليون حويصلة هوائية في رئتينا مع كل زفير؟ ولماذا تنفجر فقاعات الماء النقي أو الكحول النقي في أجزاء من الألف من الثانية بينما تمنح كمية ضئيلة من المادة الخافضة للتوتر السطحي (الصابون) الغشاء مرونةً ذاتية الالتئام؟ وكيف سقطت حدسية اللورد كلفن التي صمدت 107 أعوام حول تقسيم الفضاء ثلاثي الأبعاد بأقل مساحة سطحية عام 1994 أمام هيكل ويري-فيلان؟

في هذا المقال نجمع بين الديناميكا الحرارية، والهندسة التفاضلية، وحساب التغيرات، ونظرية القياس الهندسي،وميكانيكا الموائع البينية لكشف الفيزياء الرياضية الكاملة لفقاعات الصابون والرغوة.

1. الديناميكا الحرارية للتوتر السطحي وعدم التماثل الجزيئي

في قلب السائل (Bulk) داخل الماء، يكون كل جزيء $\text{H}_2\text{O}$ محاطاً بشكل متماثل في جميع الاتجاهات بجزيئات مجاورة ترتبط معه عبر شبكة من الروابط الهيدروجينية وقوى فان دير فالس. ولذلك تتلاشى محصلة القوى المتوسطة زمنياً، ويستقر الجزيء في حالة طاقة وضع تماسكية منخفضة (عدد التناسق $z_b \approx 12$، وطاقة التماسك $-u_0$).

أما عند السطح البيني بين السائل والبخار، فينكسر هذا التماثل بشدة؛ إذ تقل كثافة الجزيئات في الطور الغازي بنحو ثلاث رتب أسية عنها في الطور السائل. ونتيجةً لذلك، يفقد الجزيء الموجود على السطح قرابة نصف جيرانه الجاذبين ($z_s < z_b$) ويمتلك طاقة وضع فائضة. ولنقل الجزيئات من داخل السائل إلى السطح لزيادة مساحة السطح البيني $A$ بمقدار $dA$، لا بد من بذل شغل عكوس موجب ضد قوى التماسك.

يُعطى التفاضل الكلي لطاقة غيبس الحرة $G(T, P, A, \{N_i\})$ بالعلاقة:

ومنها نحصل على التعريف الديناميكي الحراري الدقيق للتوتر السطحي $\gamma$:

ويحمل التوتر السطحي معنيين فيزيائيين متكافئين:

1. الطاقة الحرة الفائضة لكل وحدة مساحة ($\text{J/m}^2$): عند ثبات الحجم $V$، يتجه النظام تلقائياً نحو تقليل طاقته الحرة $\gamma A$ عبر تقليل مساحته السطحية $A$. وبموجب المتباينة المتساوية المحيط (Isoperimetric Inequality) في الفضاء ثلاثي الأبعاد $A^3 \ge 36\pi V^2$، فإن الكرة هي الشكل الهندسي الوحيد الذي يحقق أصغر مساحة سطحية لحجم معين.
2. قوة الشد المماسية لكل وحدة طول ($\text{N/m}$): وفق صيغة كيركوود-باف (1949) في الميكانيكا الإحصائية، يُعبر عن $\gamma$ بتكامل عدم تماثل موتر الضغط بين المركبة العمودية $P_N$ والمركبة المماسية $P_T(z)$ عبر طبقة الانتقال البينية: $$\gamma = \int_{-\infty}^{+\infty} \bigl[ P_N - P_T(z) \bigr]\,dz$$

2. الهندسة التفاضلية للأسطح المنحنية ومعادلة يونغ-لابلاس

عندما يكون السطح البيني منحنياً، تولد قوى الشد المماسية $\gamma$ مركبة عمودية صافية تتجه نحو مركز الانحناء (الجانب المقعر)، مما يفرض قفزة في الضغط $\Delta P = P_{\text{in}} - P_{\text{out}} > 0$. ليكن $R_1, R_2$ نصفي قطر الانحناء الرئيسيين ($\kappa_1 = 1/R_1, \kappa_2 = 1/R_2$) لسطح أملس $\Sigma \subset \mathbb{R}^3$. يُعرَّف متوسط الانحناء $H$ وانحناء غاوس $K$ بالعلاقتين:

عند إزاحة رقعة متناهية الصغر $dA = R_1 R_2\,d\theta_1\,d\theta_2$ إزاحةً افتراضية عمودية $\delta\xi$، يكون تغير المساحة من الرتبة الأولى $\delta(dA) = dA\left(\frac{1}{R_1} + \frac{1}{R_2}\right)\delta\xi$ والحجم الممسوح $\delta(dV) = dA\,\delta\xi$. ومن شرط الاتزان الترموديناميكي $\delta F = \gamma\,\delta(dA) - \Delta P\,\delta(dV) = 0$ نشتق معادلة يونغ-لابلاس (1805–1806):

بينما تمتلك قطرة الماء في الهواء سطحاً بينياً واحداً فقط ($\Delta P_{\text{droplet}} = 2\gamma/R$)، فإن فقاعة الصابون في الهواء عبارة عن قشرة سائلة رقيقة بسمك $h \ll R$ تمتلك سطحين بينيين (غاز-سائل) (داخلي بنصف قطر $R$ وخارجي بنصف قطر $R+h$)، ومن هنا يظهر المعامل $4$:

3. الأسطح الدنيا ($H = 0$)، وكارثة السطح السلسلي (Catenoid)، ومسألة بلاتو

عندما يمتد غشاء صابون مفتوح عبر إطار سلكي دون حبس حجم من الهواء، يتعرض كلا وجهيه لنفس الضغط الجوي ($\Delta P = 0$). وتفرض معادلة يونغ-لابلاس عندئذٍ أن ينعدم متوسط الانحناء في كل نقطة:

تُعرف الأسطح التي تحقق $H \equiv 0$ باسم الأسطح الدنيا (Minimal Surfaces). وفي التمثيل الديكارتي $z = u(x,y)$، يؤدي تصغير دالي المساحة $A[u] = \iint \sqrt{1+u_x^2+u_y^2}\,dx\,dy$ إلى معادلة أويلر-لاغرانج التفاضلية الجزئية غير الخطية (لاغرانج، 1760):

كل نقطة غير مستوية على غشاء الصابون المفتوح هي محلياً نقطة سرج زائدية ($K < 0$). والسطح الأدنى الدوراني الممتد بين حلقتين دائريتين متحدتي المحور بنصف قطر $R$ والمسافة بينهما $2h$ هو السطح السلسلي (الكاتينويد، Catenoid) الذي اكتشفه أويلر عام 1744:

عند إبعاد الحلقتين تدريجياً لما يتجاوز النسبة الحرجة ($\coth u_c = u_c \approx 1.19968$):

يختفي الحل الحقيقي عبر تشعب العقدة والسرج، فينقبض عنق الغشاء الصابوني وينقطع قافزاً إلى قرصين مستويين منفصلين (حل غولدشميت غير المستمر). وقد حُلت مسألة بلاتو العامة عام 1930 على يد جيسي دوغلاس (الحائز على أول ميدالية فيلدز عام 1936) وتيبور رادو.

4. مفارقة الفقاعتين، ونضج أوستفالد، ومادة السرفكتانت الرئوية (DPPC)

إذا وصلنا فقاعة صابون صغيرة ($R_{\text{small}}$) بفقاعة كبيرة ($R_{\text{large}}$) عبر أنبوب مزود بصمام، فإنه بمجرد فتح الصمام تتقلص الفقاعة الصغيرة بسرعة وتفرغ هواءها داخل الفقاعة الكبيرة فتزداد كِبَراً! فبحسب القانون $\Delta P = 4\gamma/R$، يتناسب الضغط الداخلي عكسياً مع نصف القطر:

  flowchart TD
    A["فقاعتا صابون متصلتان بأنبوب ذي صمام:<br/>فقاعة صغيرة (R_small) وفقاعة كبيرة (R_large)"] --> B["معادلة يونغ-لابلاس: ΔP = 4γ / R<br/>الضغط الفائض يتناسب عكسياً مع نصف القطر R"]
    B --> C["ضغط الفقاعة الصغيرة: ΔP_small = 4γ / R_small<br/>ضغط الفقاعة الكبيرة: ΔP_large = 4γ / R_large"]
    C --> D["اختلال الضغط: ΔP_small > ΔP_large<br/>(الانحناء الأشد يولد ضغطاً داخلياً أعلى!)"]
    D --> E["فتح الصمام: يتدفق الغاز من الفقاعة الصغيرة ذات الضغط العالي<br/>إلى الفقاعة الكبيرة ذات الضغط المنخفض"]
    E --> F["تغذية راجعة إيجابية:<br/>تتقلص الفقاعة الصغيرة أكثر فيرتفع ضغطها وتنهار تماماً"]
    D --> G["أزمة فسيولوجية في رئات الثدييات:<br/>300 مليون حويصلة هوائية بأنصاف أقطار متفاوتة (50–200 μm)"]
    G --> H["الخلايا السنخية من النوع الثاني تفرز السرفكتانت الرئوي (DPPC)"]
    H --> I["عند تقلص الحويصلة تنضغط طبقة DPPC الأحادية بشدة،<br/>مما يخفض التوتر السطحي γ إلى 0.5–2 mN/m"]
    I --> J["تتعادل ضغوط لابلاس (2γ_small / R_small ≈ 2γ_large / R_large)<br/>مما يمنع انخماص الرئة وانهيار الحويصلات!"]

وفي الرغوة الحقيقية، يؤدي انتشار الغاز عبر الأغشية السائلة من الفقاعات الصغيرة إلى الكبيرة إلى خشونة الرغوة أو نضج أوستفالد (Ostwald Ripening). وفي عام 1952، أثبت جون فون نيومان باستخدام مبرهنة غاوس-بونيه أنه في الرغوة الجافة ثنائية الأبعاد، تتغير مساحة الخلية المضلعة ذات $n$ أضلاع وفق قانون فون نيومان:

فالخلايا التي لها $n < 6$ أضلاع تنكمش وتختفي، والمسدسات ($n = 6$) تظل ثابتة تماماً، والخلايا التي لها $n > 6$ تنمو باستمرار (وتم تعميم القانون لثلاثة أبعاد بواسطة ماكفرسون وسرولوفيتز عام 2007).

تحتوي رئتا الإنسان على نحو 300 مليون حويصلة هوائية (Alveoli) ($R \approx 50\text{–}200\,\mu\text{m}$). ولو كان التوتر السطحي $\gamma$ ثابتاً، لانهارت الحويصلات الصغيرة فوراً داخل الكبيرة (الانخماص الرئوي، Atelectasis). غير أن الخلايا السنخية من النوع الثاني تمنع ذلك بإفراز السرفكتانت الرئوي الغني بمركب ثنائي بالميتويل فوسفاتيديل كولين (DPPC)؛ فعند الزفير وتقلص الحويصلة، تنضغط طبقة DPPC غير الذائبة في الماء بشدة فينخفض $\gamma(R)$ إلى نحو $0.5\text{–}2\,\text{mN/m}$، مما يحفظ استقرار جميع الحويصلات.

5. قوانين جوزيف بلاتو (1873)، وبرهان جين تايلور (1976)، ومبرهنة الفقاعة المزدوجة

صاغ الفيزيائي البلجيكي الكفيف جوزيف بلاتو عام 1873 القوانين الهندسية لاتزان الرغوة، وبرهنتها رياضياً جين تايلور (Jean E. Taylor) عام 1976 باستخدام نظرية القياس الهندسي:

1. قانون بلاتو الأول (3 أغشية بزاوية $120^\circ$): تلتقي أغشية الصابون حصرياً ثلاثةً ثلاثةً على طول حدود بلاتو بزوايا متساوية قدرها $120^\circ$ ($2\pi/3$) ($\sum_{i=1}^3 \hat{\mathbf{t}}_i = \mathbf{0} \implies \cos\theta = -1/2$). أما تقاطع 4 أغشية بزاوية $90^\circ$ ($L_{\text{cross}} = 4\sqrt{2} \approx 5.657$) فهو غير مستقر وينشطر فوراً إلى عقدتين ثلاثيتين بزاوية $120^\circ$ وفق شجرة شتاينر ($L_{\text{Steiner}} = 2(1+\sqrt{3}) \approx 5.464$، انتقال طوبولوجي T1).
2. قانون بلاتو الثاني (4 حواف بزاوية $109.47^\circ$): تلتقي حدود بلاتو حصرياً أربعاً أربعاً عند الرؤوس وفق زاوية مارالدي الرباعية الأوجه: $$\theta_{\text{Maraldi}} = \arccos\!\left(-\frac{1}{3}\right) \approx 109.47122^\circ \quad \left(109^\circ 28' 16''\right)$$
3. مبرهنة الفقاعة المزدوجة (Hutchings, Morgan, Ritoré, Ros 2000): عند التحام فقاعتين بنصفي قطر $R_1 > R_2$، ينحني الحاجز المشترك (Septum) نحو داخل الفقاعة الكبرى بنصف قطر انحناء $R_{\text{septum}}$ يحقق: $$\frac{1}{R_{\text{septum}}} = \frac{1}{R_2} - \frac{1}{R_1} \quad \Longleftrightarrow \quad R_{\text{septum}} = \frac{R_1 R_2}{R_1 - R_2}$$

6. حدسية كلفن (1887) في مواجهة هيكل ويري-فيلان (1994)

في عام 1887، طرح اللورد كلفن مسألة تقسيم الفضاء $\mathbb{R}^3$ إلى خلايا متساوية الحجم ($V=1$) بأقل مساحة سطحية. وبسبب الإحباط الهندسي، اقترح كلفن ثماني الأوجه المقطوع (ذو 14 وجهاً: 6 مربعات و8 مسدسات، خلية فورونوي لشبكة BCC) المنحني قليلاً بمساحة معيارية $A_{\text{Kelvin}} \approx 5.30628$.

وبعد 107 أعوام، نقض الفيزيائيان دينيس ويري وروبرت فيلان (كلية ترينيتي في دبلن) عام 1994 حدسية كلفن باكتشاف هيكل ويري-فيلان (Weaire-Phelan Structure) المستوحى من هيدرات الكلاثرات من النوع الأول، وتتكون خليته الأولية المكعبة من 8 كثيرات أوجه متساوية الحجم من نوعين: اثنين من اثني عشري الأوجه الخماسي غير المنتظم ($5^{12}$) وستة من أربعة عشري الأوجه ($5^{12}6^2$ ذي وجهين مسدسين و12 وجهاً خماسياً). ولأن $88.9\%$ من الأوجه خماسية (زاويتها الداخلية $108^\circ$ قريبة جداً من $109.47^\circ$)، تنخفض المساحة بنسبة $0.34\%$ إلى $A_{\text{WP}} \approx 5.28825$، وهو الهيكل الذي شُيِّد به المركز الوطني للرياضات المائية في بكين («مكعب الماء»).

7. مرونة غيبس-مارانغوني، والتداخل الضوئي، وضغط الانفصال DLVO

إن انخفاض التوتر السطحي وحده لا يكفي لاستقرار الفقاعة (فالإيثانول النقي ينفجر فوراً). بل يكمن سر بقاء غشاء الصابون في مرونة غيبس-مارانغوني (Gibbs-Marangoni Elasticity):

  flowchart TD
    A["اضطراب خارجي (تيار هواء، تصريف جاذبي، موجة حرارية)<br/>يشد ويرقق غشاء الصابون محلياً"] --> B["تتوسع المساحة المحلية ΔA؛<br/>فتنخفض كثافة امتزاز السرفكتانت السطحية Γ محلياً"]
    B --> C["انخفاض Γ يؤدي إلى ارتفاع التوتر السطحي المحلي<br/>γ_thin فوق γ_thick المحيط به!"]
    C --> D["ينشأ تدرج مماسي في التوتر السطحي ∇_s γ<br/>(إجهاد قص مارانغوني τ_M = ∇_s γ يتجه نحو البقعة الرقيقة)"]
    D --> E["حمل مارانغوني (Marangoni Flow):<br/>المنطقة ذات γ المرتفع تسحب الطبقة الأحادية وتجر معها الماء السائل باللزوجة"]
    E --> F["يلتئم سمك الغشاء h ذاتياً في أجزاء من الألف من الثانية،<br/>مما يمنع التمزق الشعري!"]
    F --> G["يستمر التصريف الجاذبي في ترقيق الغشاء إلى أطوال موجية مرئية (100–700 nm)"]
    G --> H["التداخل الضوئي للأغشية الرقيقة يولد أشرطة قوس قزح المتلألئة"]
    H --> I["عند h < 30 nm، يوازن ضغط الانفصال DLVO Π(h) تجاذب فان دير فالس؛<br/>والانقلاب الطوري π يكوّن 'غشاء نيوتن الأسود' (h ≈ 4 nm)"]

عندما يتمدد الغشاء ويرق محلياً، تنخفض كثافة السرفكتانت السطحية $\Gamma$ فيرتفع التوتر السطحي المحلي $\gamma_{\text{thin}}$. ويعمل إجهاد مارانغوني $\boldsymbol{\tau}_M = \nabla_s \gamma$ ومعامل مرونة غيبس

على سحب السائل فوراً نحو المنطقة المرققة لإصلاحها (وهو نفس التأثير المسؤول عن «دموع النبيذ» على جدران الكأس). وعندما يتناقص السمك $h$ إلى نطاق الضوء المرئي، يحدث تداخل الأغشية الرقيقة ($2n_{\text{film}}h\cos\theta_t + \lambda/2 = m\lambda$). وتحت $h \ll \lambda/4 \approx 30\,\text{nm}$، يؤدي انقلاب الطور بمقدار $\pi$ عند السطح الأمامي إلى تداخل هدام لجميع الأطوال الموجية المرئية فيبدو الغشاء أسود تماماً (غشاء نيوتن الأسود، $h \approx 4\,\text{nm}$)، محمياً من الانهيار بفضل ضغط الانفصال (Disjoining Pressure) وفق نظرية DLVO $\Pi(h) = \Pi_{\text{vdW}}(h) + \Pi_{\text{el}}(h) + \Pi_{\text{sr}}(h)$.