الحجز المغناطيسي للبلازما المتراكمة هو آلية شائعة للأجسام المدمجة، لكنه حتى الآن دُرس بالتفصيل بشكل شبه حصري على النجوم النيوترونية. الأقزام البيضاء في الأنظمة الثنائية القريبة، مثل القطبية (polars) والقطبية الوسيطة (intermediate polars)، توجه المادة أيضًا إلى أقطابها المغناطيسية، مُنشئة أعمدة تراكم موضعية. الرصدات الأخيرة للانفجارات الصغرى (micronovae) على أقزام بيضاء شديدة المغناطيسية أكدت أن الاحتراق النووي الحراري يمكن أن يحدث في مناطق محدودة عند الأقطاب. ومع ذلك، ظل السؤال حول ما إذا كانت الأقزام البيضاء قادرة على احتجاز بنى كبيرة مدعومة مغناطيسيًا مفتوحًا. ينقل هذا العمل الشكلية المطورة للنجوم النيوترونية إلى الظروف الفيزيائية للأقزام البيضاء.
يحل المؤلفون معادلة غراد-شافرانوف لتوازن متماثل محوري مع معادلة حالة متساوية الحرارة. تُطبق طريقة تكرارية متسقة ذاتيًا حيث تُحدد دالة المصدر من شرط تجمّد التدفق المغناطيسي: الكتلة لا تعبر الأسطح المغناطيسية. تُستخدم شبكة بحجم 256×256 مع تكثيف لوغاريتمي عند السطح. على عكس طريقة الاسترخاء الزائد المتعاقب التقليدية، تُحل المعادلة باستخدام حلال متفرق مباشر UMFPACK. للعثور على الكتلة القصوى المحتجزة، يُستخدم بحث ثنائي مع التحقق من التقارب عبر البواقي وخطأ الكتلة. تتباين المعلمات في النطاقات: المجال المغناطيسي من 10^5 إلى 3×10^8 غاوس، ودرجة الحرارة من 2×10^4 إلى 10^7 كلفن.
بالنسبة للمعلمات النموذجية (M=0.8 M⊙، R=7×10^8 سم، B=5×10^7 غاوس، T=1.59×10^5 كلفن)، يبلغ ارتفاع الجبل المميز حوالي 10^5 سم لكتلة 10^-3 M⊙. تصل الكثافة القصوى عند القاعدة إلى 2.48×10^7 غ/سم^3. تتشوه خطوط المجال المغناطيسي بشكل ملحوظ وتزاح نحو خط الاستواء، لكن التوازن يُحفظ. تخضع الكتلة القصوى التي يمكن احتجازها لقانون قوى: M_max = 4.07×10^-3 M⊙ (T/10^5 K)^-2 (B/5×10^7 غاوس)^2. وهذا يعادل حوالي 3% من الكتلة الحرجة التي حُصل عليها سابقًا تحليليًا للتراكمات الصغيرة. يقلل الدفن المغناطيسي العزم ثنائي القطب إلى 10% من قيمته الأولية، مع كون هذا الانخفاض النسبي مستقلًا عن شدة المجال. يبقى المجال ثنائي القطب في الغالب، لكن تُثار لحظات متعددة الأقطاب فردية الرتبة (l=3,5,7,9).
يقدم هذا العمل لأول مرة تقديرًا كميًا للجبال المحتجزة مغناطيسيًا على الأقزام البيضاء، مما يسد الفجوة بين نظرية النجوم النيوترونية ورصد المتغيرات الكارثية. على عكس التقديرات السابقة القائمة على نماذج كروية أو أسطوانية بسيطة، يُحل هنا معادلة التوازن الكاملة مع مراعاة متسقة ذاتيًا لتجمّد التدفق. وهذا يسمح بتحديد موثوق للكتلة القصوى ودرجة الدفن المغناطيسي. تُظهر النتائج أن الدفن المغناطيسي على الأقزام البيضاء أضعف منه على النجوم النيوترونية: إذ ينخفض المجال بنسبة 10% فقط، وليس بمقدار مراتب عشرية. يمكن أن يكون هذا النهج أساسًا لتفسير البيانات الرصدية حول المجالات المغناطيسية والتراكم في القطبية والقطبية الوسيطة.
يمكن توسيع النموذج بالتخلي عن معادلة الحالة متساوية الحرارة وإدراج ملفات درجة حرارة واقعية. إن مراعاة غوص الجبل في الغلاف الجوي للنجم واختلاط المادة، كما فعلت دراسات النجوم النيوترونية، سيسمح بتحسين البنية. إدخال الموصلية المحدودة والتأثيرات غير المثالية سيفتح الطريق لدراسة التبدد وإعادة الاتصال. المحاكاة الديناميكية المبنية على حالات التوازن التي حُصل عليها ستساعد في تقييم استقرار الجبال وتطورها. من المثير للاهتمام بشكل خاص مراعاة الدوران ونظام المروحة (propeller)، الذي قد يحد من تراكم الكتلة في الأنظمة غير المتزامنة.
النتائج قابلة للتطبيق على تفسير رصد المتغيرات الكارثية المغناطيسية، بما في ذلك القطبية والقطبية الوسيطة. كما أنها توفر أساسًا لتقدير إهليلجية الأقزام البيضاء واحتمال إصدار موجات جاذبية مستمرة يمكن أن تكتشفها المقاييس التداخلية الفضائية المستقبلية.
الخطوة التالية ستكون إدراج معادلة حالة واقعية ومراعاة غوص المادة المتراكمة في الغلاف الجوي للنجم. من الضروري أيضًا إجراء محاكاة هيدروديناميكية للتحقق من استقرار حالات التوازن التي تم العثور عليها ودراسة تطور عدم الاستقرار.
يربط هذا العمل نظرية الحجز المغناطيسي بالمظاهر الرصدية للتراكم على الأقزام البيضاء، مثل الانفجارات الصغرى وتغيرات المجالات المغناطيسية. يبقى السؤال مفتوحًا حول التطور طويل الأمد للمجال المغناطيسي تحت التراكم المطول ودور عدم الاستقرار في تدمير الجبال.
يمكن أن تمتلك الأقزام البيضاء مجالات مغناطيسية تتراوح من 10^3 إلى 10^9 غاوس، أي بفارق مليون ضعف. عند مجال 10^8 غاوس، تصل الكتلة القصوى المحتجزة إلى 0.067 من كتلة الشمس، وهو ما يعادل كتلة كوكب عملاق بأكمله.