1. لماذا لا توجد مغانط سائلة في الطبيعة؟ — حاجز درجة حرارة كوري

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

ولكن لماذا لا يمكننا ببساطة صهر مغناطيس حديدي عادي في بوتقة للحصول على «مغناطيس سائل»؟ تكمن الإجابة في قانون أساسي من قوانين الفيزياء الكمومية والإحصائية: درجة حرارة كوري ($T_C$).

تنشأ المغناطيسية الحديدية (Ferromagnetism) من تفاعل التبادل الكمومي لهايزنبرغ (Exchange Interaction) الذي يصطف بموجبه اللف المغزلي (Spin) للإلكترونات المتجاورة $\mathbf{S}_i$ و$\mathbf{S}_j$ بشكل متوازٍ:

إن اقتران التبادل هذا حساس للغاية للمسافة الذرية البينية وللتهيج الحراري. فإذا سخّنا معدناً مغناطيسياً حديدياً فوق درجة حرارة كوري $T_C$ الخاصة به، تتغلب الطاقة الحرارية $k_B T$ على طاقة التبادل $J_{\text{ex}}$، فينهار الاصطفاف التلقائي للسبينات ويتحول المعدن إلى مادة مغناطيسية مسايرة (Paramagnetic) عشوائية.

وكما يوضح الجدول التالي، فإن درجة حرارة كوري $T_C$ لجميع العناصر والأكاسيد المغناطيسية الحديدية المعروفة تقع أقل بمئات الدرجات من نقطة انصهارها $T_m$:

يفقد الحديد المنصهر مغناطيسيته الحديدية بالكامل قبل أن ينصهر بـ $768\,\text{K}$ ($T_C \ll T_m$). لذلك لا يوجد في الطبيعة أي سائل متجانس أحادي الطور يمتلك مغناطيسية حديدية عند درجة حرارة الغرفة.

كيف تمكن إذن مهندس وكالة ناسا (NASA) ستيفن بابيل (Stephen Papell) في عام 1963 من ابتكار أول سائل مغناطيسي في العالم لضخ وقود الصواريخ السائل في انعدام الجاذبية باستخدام المجالات المغناطيسية؟ تمثّل الحل العبقري لبابيل في الالتفاف حول حاجز كوري عبر تصميم معلق غروي ثلاثي المكونات:

1. جسيمات نانوية مغناطيسية ($\sim 5\%$ حجماً): بلورات نانوية صلبة من المغنيتيت ($\text{Fe}_3\text{O}_4$) أو الماغهيميت ($\gamma\text{-Fe}_2\text{O}_3$) بقطر نموذجي $d \approx 10\,\text{nm}$، تعمل تحت درجة حرارة كوري بكثير وتحتفظ بتمغنطها القوي.
2. طبقة خافض للتوتر السطحي / Surfactant ($\sim 10\%$ حجماً): غلاف جزيئي بسمك $\delta \approx 2\,\text{nm}$ (مثل حمض الأوليك $\text{C}_{18}\text{H}_{34}\text{O}_2$)، يرتبط رأسه القطبي كيميائياً بسطح الجسيم النانوي بينما يمتد ذيله الهيدروكربوني في المذيب.
3. سائل حامل ($\sim 85\%$ حجماً): مذيب عضوي (كيروسين، زيت اصطناعي، إستر) أو ماء، حيث تحافظ الحركة البراونية الحرارية على الجسيمات النانوية معلقة إلى ما لا نهاية.

2. الفيزياء النانوية للمغناطيسية: حد النطاق الواحد لكيتل والمغناطيسية المسايرة الفائقة للانجفان

لماذا يجب أن يكون قطر جسيمات المغنيتيت قريباً تحديداً من $d \approx 10\,\text{nm}$؟ لو قمنا بمجرد خلط مسحوق حديد ميكرومتري ($d \approx 1\text{–}10\,\mu\text{m}$) في الزيت، لترسبت الحبيبات فوراً ولتكتلت في كتلة صلبة عند تقريب مغناطيس.

2.1 القطر الحرج لكيتل للجسيمات أحادية النطاق

ينقسم المغناطيس العياني إلى عدة نطاقات مغناطيسية لفايس (Weiss Domains) تفصل بينها جدران بلوخ لتقليل الطاقة المغناطيسية الساكنة للمجال المزيل للتمغنط ($\propto d^3$). غير أن إنشاء جدار نطاق يستهلك طاقة تبادل وتباين مناحٍ تتناسب مع مساحة المقطع ($\propto d^2$). وتحت قطر حرج $d_{\text{sd}}$ (حَسَبه تشارلز كيتل Charles Kittel)، يصبح تكوين الجدار غير مجدٍ طاقياً:

حيث $A_{\text{ex}} \approx 10^{-11}\,\text{J/m}$ صلابة التبادل، $K \approx 1.35 \times 10^4\,\text{J/m}^3$ ثابت التباين المغناطيسي البلوري، و$M_d \approx 4.46 \times 10^5\,\text{A/m}$ تمغنط التشبع للمغنيتيت الكتلي $\text{Fe}_3\text{O}_4$. بالنسبة للمغنيتيت، يبلغ $d_{\text{sd}} \approx 25\text{–}30\,\text{nm}$. وبالتالي فإن جسيماً نانوياً بقطر $d = 10\,\text{nm}$ هو حتماً جسيم أحادي النطاق (Single-Domain Particle): تدور جميع سبيناته الذرية البالغ عددها $\sim 10^4$ معاً كثنائي قطب مغناطيسي عملاق واحد $m = M_d V = M_d \frac{\pi}{6}d^3 \approx 2.3 \times 10^{-19}\,\text{A}\cdot\text{m}^2$ ($\sim 25\,000$ مغنيطون بور $\mu_B$).

2.2 استرخاء نيل وبراون، والمغناطيسية المسايرة الفائقة

هل يظل هذا العزم العملاق $m$ مقيداً على طول المحور السهل للبلورة؟ عند $d = 10\,\text{nm}$، يكون الحجم $V \approx 5.24 \times 10^{-25}\,\text{m}^3$ ضئيلاً جداً بحيث يكون حاجز طاقة التباين $E_B = K V \approx 7 \times 10^{-21}\,\text{J}$ مقارباً للطاقة الحرارية $k_B T \approx 4.14 \times 10^{-21}\,\text{J}$ عند درجة حرارة الغرفة ($T = 300\,\text{K}$)!

ونتيجة لذلك، ينقلب العزم المغناطيسي تلقائياً وباستمرار بفعل التقلبات الحرارية عبر آليتين متوازيتين:

1. استرخاء نيل ($\tau_N$): يدور متجه العزم المغناطيسي داخلياً ضمن الشبكة البلورية الساكنة متجاوزاً الحاجز $K V$: $$\tau_N = \tau_0 \exp\left(\frac{K V}{k_B T}\right)$$ ومع $\tau_0 \approx 10^{-9}\,\text{s}$، يبلغ زمن الاسترخاء لجسيم مغنيتيت بقطر $10\,\text{nm}$ قيمة فائقة القصر $\tau_N \approx 5 \times 10^{-9}\,\text{s}$ (بضع نانوثوانٍ).
2. استرخاء براون ($\tau_B$): يدور الجسيم النانوي بأكمله مع غلافه الجزيئي (الحجم الهيدروديناميكي $V_H = \frac{\pi}{6}(d + 2\delta)^3$) فيزيائياً داخل المائع ذي اللزوجة $\eta$: $$\tau_B = \frac{3 \eta V_H}{k_B T}$$ ويخضع الزمن الفعال للعلاقة $\frac{1}{\tau_{\text{eff}}} = \frac{1}{\tau_N} + \frac{1}{\tau_B}$.

في غياب مجال خارجي ($\mathbf{H} = \mathbf{0}$)، تتبعثر اتجاهات العزوم حرارياً خلال نانوثوانٍ: فلا يظهر السائل المغناطيسي أي تمغنط متبقٍ ($M_r = 0$) ولا أي قوة قاهرة ($H_c = 0$) (انعدام التخلف المغناطيسي). ولكن عند تطبيق مجال $\mathbf{H}$، يستجيب كل جسيم بعزمه العملاق البالغ $25\,000\,\mu_B$، مولداً قابلية تمغنط ابتدائية $\chi_i$ أكبر بآلاف المرات من المواد البارامغناطيسية العادية. تُعرف هذه الظاهرة باسم المغناطيسية المسايرة الفائقة (Superparamagnetism)، وتُوصف بدقة بواسطة دالة لانجفان $\mathcal{L}(\xi)$:

حيث $\phi$ الكسر الحجمي للمادة المغناطيسية و$M_s = \phi M_d$ تمغنط التشبع للسائل. وفي المجالات الضعيفة ($\xi \ll 1$)، يعطي التقريب $\mathcal{L}(\xi) \approx \xi / 3$ القابلية الابتدائية $\chi_i = \frac{\mu_0 \phi M_d^2 V}{3 k_B T} \sim 1\text{–}5$.

3. الاستقرار الغروي: نظرية DLVO الموسعة والمقارنة مع الموائع المغناطيسية الريولوجية (MRF)

لمنع الجسيمات النانوية من التجلط والتلاصق غير العكوس، يجب أن يمتلك جهد التفاعل الكلي $U_{\text{total}}(l)$ كدالة للمسافة السطحية المعيارية $l = 2s / d$ حاجز تنافر مرتفعاً $\Delta U \gg 15\,k_B T$:

1. تجاذب لندن–فان دير فالس ($U_{\text{vdW}}$): تجاذب كمومي قصير المدى (ثابت هامكر $A_H \approx 10^{-19}\,\text{J}$) يتباعد إلى $-\infty$ عند التلامس ($l \to 0$): $$U_{\text{vdW}}(l) = -\frac{A_H}{6}\left[ \frac{2}{l^2 + 4l} + \frac{2}{(l+2)^2} + \ln\left(\frac{l^2 + 4l}{(l+2)^2}\right) \right]$$
2. تفاعل ثنائي القطب المغناطيسي ($U_{\text{dd}}$): في وضعية الرأس للذيل، يتجاذب ثنائيا قطب بطاقة $U_{\text{dd}}^{\text{min}}(s) = -\frac{\mu_0 m^2}{2\pi (d + s)^3}$.
3. التنافر الفراغي لماكور–روزينزويغ ($U_{\text{steric}}$): عندما تقترب الجسيمات لمسافة $s < 2\delta$، تنضغط سلاسل حمض الأوليك (كثافة التطعيم $\xi_s \approx 10^{18}\,\text{chains/m}^2$). يؤدي انخفاض إنتروبيا التشكيل للسلاسل البوليمرية إلى توليد تنافر إنتروبي قوي: $$U_{\text{steric}}(l) = \frac{\pi d^2 \xi_s k_B T}{2} \left[ 2 - \frac{l+2}{t}\ln\left(\frac{1+t}{1+l/2}\right) - \frac{l}{t} \right] \quad \left(\text{عند } l < 2t,\; t = \frac{2\delta}{d}\right)$$

هذا الدرع الفراغي بسمك $\delta \approx 2\,\text{nm}$ هو ما يميز جذرياً السائل المغناطيسي (FF) عن المائع المغناطيسي الريولوجي (MRF):

  flowchart TD
    A["جسيمات المغنيتيت النانوية في المعلق (d ≈ 10 nm)"] --> B{"توازن جهود التفاعل بين الجسيمات"}
    B -->|"قوى التجاذب"| C["تجاذب فان دير فالس U_vdW<br/>+ تجاذب ثنائي القطب المغناطيسي U_dd"]
    B -->|"قوى التنافر والحرارة"| D["تنافر فراغي إنتروبي U_steric (δ ≈ 2 nm)<br/>+ حركة براونية حرارية (k_B T)"]
    C --> E["توازن غروي وفق نظرية DLVO الموسعة:<br/>حاجز طاقة ΔU >> k_B T يمنع التكتل"]
    D --> E
    E --> F["سائل سوبر بارامغناطيسي مستقر لسنوات<br/>(لا يترسب ويحتفظ بسيولته تحت المجال القوي)"]

4. الهيدروديناميكا المغناطيسية (FHD) وعدم استقرار كاولي–روزينزويغ

لماذا ينتصب السطح الحر للسائل المغناطيسي فجأة على شكل شبكة منتظمة من القمم المخروطية عند تسليط مجال مغناطيسي عمودي؟ تأتي الإجابة من الهيدروديناميكا المغناطيسية للسوائل المغناطيسية (Ferrohydrodynamics - FHD) التي صاغها رونالد روزينزويغ.

4.1 قوة كيلفن الحجمية ومعادلة برنولي المعممة

في مائع قابل للتمغنط وغير موصل كهربائياً موضوع في مجال $\mathbf{H}$، يتأثر كل عنصر حجمي ذي تمغنط $\mathbf{M}$ بـ كثافة قوة كيلفن المغناطيسية:

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

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

حيث $p_0$ الضغط الجوي، $\sigma$ التوتر السطحي، $\kappa \approx -\nabla_\perp^2 \zeta$ متوسط انحناء السطح، و$M_n = \mathbf{M} \cdot \mathbf{n}$ المركبة العمودية للتمغنط. يعمل الحد $-\frac{1}{2}\mu_0 M_n^2$ بمثابة سحب مغناطيسي نحو الخارج: فأي نتوء صغير $\zeta(x,y) > 0$ على السطح يُركّز خطوط المجال المغناطيسي (تأثير مانعة الصواعق)، مما يزيد $M_n$ محلياً ويسحب السائل إلى الأعلى بقوة أكبر!

4.2 تحليل الاستقرار الخطي وعلاقة التشتت (كاولي وروزينزويغ، 1967)

في عام 1967، درس م. كاولي (M. D. Cowley) ور. روزينزويغ (R. E. Rosensweig) طبقة مستوية من السائل المغناطيسي تحت مجال مغناطيسي عمودي منتظم $\mathbf{B}_0 = B_0 \hat{\mathbf{z}}$ تتعرض لموجة سطحية صغيرة $\zeta(\mathbf{r}_\perp, t) = \zeta_0 e^{i(\mathbf{k}\cdot\mathbf{r}_\perp - \omega t)}$ بعدد موجي $k = |\mathbf{k}| = 2\pi / \lambda$.

يقود الحل المشترك لمعادلات ماكسويل ($\nabla \cdot \mathbf{B} = 0$، $\nabla \times \mathbf{H} = \mathbf{0}$) ومعادلة برنولي FHD إلى علاقة التشتت لكاولي–روزينزويغ الشهيرة:

حيث $\chi_{\text{eff}}$ القابلية الهندسية الفعالة. تعبر الحدود الثلاثة عن صراع بين ثلاث قوى:

• $+g k$ (خطي في $k$): قوة الإرجاع التثاقلية (الجاذبية) التي تثبت الأمواج الطويلة ($k \to 0$).
• $+\frac{\sigma}{\rho}k^3$ (تكعيبي في $k$): قوة الإرجاع الشعرية (التوتر السطحي) التي تثبت الأمواج القصيرة ($k \to \infty$).
• $-\frac{\mu_0 M_0^2}{\rho}\frac{\chi_{\text{eff}}}{1+\chi_{\text{eff}}}k^2$ (تربيعي في $k$، سالب!): الشد المغناطيسي المزيل للاستقرار عند الأطوال الموجية المتوسطة!

عندما يتجاوز التمغنط $M_0$ عتبة حرجة $M_c$، يصبح المنحنى $\omega^2(k)$ سالباً ($\omega^2 < 0$) حول عدد موجي حرج $k_c$. يصبح التردد $\omega = \pm i \gamma$ تخيلياً صرفاً وتنمو السعة $\zeta(t) \propto e^{+\gamma t}$ أسياً!

ومن شرطي العتبة $\omega^2(k_c) = 0$ و$\left.\frac{\partial (\omega^2)}{\partial k}\right|_{k = k_c} = 0$، نستنتج المعاملات الحرجة لـ عدم استقرار روزينزويغ:

1. العدد الموجي الحرج $k_c$ والطول الموجي الحرج $\lambda_c$ (المسافة بين الأشواك): $$k_c = \sqrt{\frac{\rho g}{\sigma}} = \frac{1}{\ell_c} \implies \lambda_c = \frac{2\pi}{k_c} = 2\pi \sqrt{\frac{\sigma}{\rho g}}$$ لسائل مغناطيسي نموذجي ($\rho \approx 1200\,\text{kg/m}^3$، $\sigma \approx 0.026\,\text{N/m}$، $g = 9.81\,\text{m/s}^2$)، يبلغ الطول الشعري $\ell_c \approx 1.5\,\text{mm}$ والمسافة الفاصلة بين قمتين متجاورتين $\lambda_c \approx 9.3\,\text{mm}$ (حوالي $1\,\text{cm}$)!
2. التمغنط الحرج $M_c$ لظهور الأشواك: $$M_c^2 = \frac{2}{\mu_0}\left(1 + \frac{1}{\chi_{\text{eff}}}\right)\sqrt{\rho g \sigma}$$

4.3 لماذا تتشكل شبكة سداسية؟ (الرنين الثلاثي اللاخطي)

تحدد النظرية الخطية مقدار $k_c$، لكنها لا تحدد اتجاه المتجه $\mathbf{k}$. لماذا تنتظم الأشواك في شبكة سداسية ثنائية الأبعاد (Hexagonal Lattice) متقنة؟

في المفكوك اللاخطي للطاقة الحرة السطحية $\mathcal{F}[\zeta]$، يكسر المجال المغناطيسي العمودي تماثل الانعكاس بين القمة والقاع ($\zeta \to -\zeta$)، مولداً حداً تكعيبياً للتفاعل الرنيني بين ثلاث متجهات موجية:

ولكي يكون مجموع ثلاثة متجهات متساوية الطول $|\mathbf{k}_1| = |\mathbf{k}_2| = |\mathbf{k}_3| = k_c$ مساوياً للصفر ($\mathbf{k}_1 + \mathbf{k}_2 + \mathbf{k}_3 = \mathbf{0}$)، يجب أن تشكل مثلثاً متساوي الأضلاع، أي أن تكون الزاوية بينها $120^\circ$ بالضبط! إن تراكب ثلاث موجات مستوية بزاوية $120^\circ$ يولّد في الفضاء الحقيقي النمط السداسي الشهير لأشواك روزينزويغ.

  flowchart TD
    A["سطح مستوٍ للسائل المغناطيسي تحت مجال عمودي B_0"] --> B["اضطراب شعري حراري ζ(x,y) بعدد موجي k"]
    B --> C["التنافس في علاقة التشتت ω²(k):<br/>+ g·k (الجاذبية) + (σ/ρ)·k³ (التوتر السطحي)<br/>- (μ₀·M₀²/ρ)·[χ/(1+χ)]·k² (الشد المغناطيسي)"]
    C --> D{"هل يتجاوز التمغنط M_0 العتبة الحرجة M_c؟"}
    D -->|"لا (M_0 < M_c)"| E["ω²(k) > 0 لجميع قيم k:<br/>يبقى السطح الحر أملس ومستقراً"]
    D -->|"نعم (M_0 > M_c)"| F["ω²(k_c) < 0 عند k_c = √(ρg/σ):<br/>نمو أسي لقمم الموجات"]
    F --> G["رنين ثلاثي لاخطي (k₁ + k₂ + k₃ = 0 بزاوية 120°)<br/>تشكل الشبكة السداسية للأشواك (λ_c ≈ 1 cm)"]

5. التطبيقات المتقدمة في الصناعة والطب والفن الحركي

تؤدي السوائل المغناطيسية دوراً حيوياً في العديد من التقنيات الحديثة:

1. أختام الفراغ الدوارة في صناعة أشباه الموصلات والأقراص الصلبة: في غرف الفراغ الفائق (UHV) لتصنيع الرقائق الإلكترونية، يثبت مغناطيس حلقي حلقة سائلة من السائل المغناطيسي في الفجوة الميكرومترية حول الأعمدة الدوارة بسرعة $10\,000\,\text{rpm}$، مما يحقق إحكاماً تاماً دون احتكاك صلب أو تآكل.
2. تبريد وتخميد مكبرات الصوت عالية الدقة (Hi-Fi): يوضع السائل المغناطيسي في الفجوة المغناطيسية حول ملف الصوت (Voice Coil)، فينقل الحرارة بكفاءة تفوق الهواء بـ $5$ مرات ويخمد بلزوجته الرنين الميكانيكي غير المرغوب فيه للغشاء.
3. فرط الحرارة المغناطيسي وتوجيه الأدوية في علاج الأورام: تُوجه جسيمات $\text{Fe}_3\text{O}_4$ المتوافقة حيوياً نحو الورم وتُعرض لمجال مغناطيسي متناوب عالي التردد ($f \approx 100\text{–}500\,\text{kHz}$)، فتبدد حرارة عبر استرخاء نيل وبراون ترفع حرارة الورم موضعياً إلى $42\text{–}45\,^\circ\text{C}$ لتدمير الخلايا السرطانية انتقائياً.
4. الفن الحركي (ساتشيكو كوداما Sachiko Kodama): تقوم الفيزيائية والفنانة اليابانية ساتشيكو كوداما (Protrude, Flow) بمزامنة مغانط كهربائية مع الموسيقى لنحت السوائل المغناطيسية في أبراج عضوية من الأشواك تتراقص في الوقت الفعلي.

إن كل شوكة تنتصب على سطح السائل المغناطيسي هي تجسيد مرئي للتوازن الدقيق بين تفاعل التبادل الكمومي لهايزنبرغ، وإنتروبيا الغرويات، ومعادلات ماكسويل ونافييه-ستوكس المقترنة.