الجزيئات والبلازما: كيف قام العلماء بتعليم الإلكترونات احترام التبعية

الهيدروجين، أبسط عنصر في الكون، يحمل المفتاح لفهم المادة في الحالات القصوى. في باطن النجوم، في منشآت الاندماج النووي الحراري الخاضعة للرقابة، أو في تجارب الضغط بالليزر، توجد المادة على شكل بلازما – “حساء” ساخن من الأيونات والإلكترونات. للتنبؤ بسلوك مثل هذه البلازما، من الضروري معرفة خصائصها الأساسية: معادلة الحالة، ودرجة التأين، وآلية تكوين الذرات والجزيئات. من الصعب للغاية حساب هذه المعلمات نظريًا للبلازما غير المثالية نظرًا للحاجة إلى مراعاة التأثيرات الكمية والتفاعلات الجماعية بين عدد كبير من الجسيمات. المعيار الذهبي هنا هو طريقة مونت كارلو الكمومية مع تكاملات المسار، والتي تعطي نتائج دقيقة للغاية. ومع ذلك، فإن سعره يمثل تكاليف حسابية هائلة عند محاولة وصف الأنظمة التي تتكون من عدد كبير من الجسيمات. ولذلك، لدراسة البلازما في نظام يكون فيه الانحطاط الكمي للإلكترونات منخفضًا، غالبًا ما يستخدم العلماء طرقًا شبه كلاسيكية أسرع. في عمليات المحاكاة هذه، تتحرك الجسيمات على طول مسارات كلاسيكية، لكن تفاعلها لا يوصف بإمكانات كولوم المعتادة، ولكن بواسطة “كمونات زائفة” خاصة تأخذ في الاعتبار مبدأ عدم اليقين الكمي والتفاعل المتبادل في النموذج. واحدة من أشهر هذه الإمكانات الزائفة تم اقتراحها من قبل غونتر كيلبج في ستينيات القرن الماضي، وتم استخدامها عدة مرات لوصف البلازما الساخنة. تم تحسينه لاحقًا ليأخذ في الاعتبار سلوك درجة الحرارة المنخفضة للنظام بشكل صحيح. لكن إمكانات كيلبج المحسنة كانت بها أيضًا عيبًا خطيرًا. عند محاكاة بلازما الهيدروجين عند درجات حرارة أقل من 50 ألف درجة، بدأت تتشكل في محاكاة النظام مجموعات مستقرة غير موجودة في الطبيعة، حيث تم ربط العديد من الإلكترونات والبروتونات في بنية واحدة مدمجة. وأدى ذلك إلى التقليل الكارثي من قيم الطاقة في عمليات المحاكاة مقارنة بالحسابات الأكثر دقة. حدد الباحثون الروس لأنفسهم مهمة العثور على مصدر هذه المشكلة وتعديل طريقة الديناميكيات الجزيئية بطريقة تتخلص من القطع الأثرية العنقودية، مع الحفاظ على الكفاءة الحسابية لهذا النهج. تم نشر العمل في مجلة فيزياء البلازما، وتم دعم البحث من قبل مؤسسة العلوم الروسية (المنحة رقم 24-19-00746). اقترح العلماء أن المشكلة تكمن في عدم مراعاة مبدأ باولي بدقة كافية للإلكترونات التي لها نفس إسقاط الدوران. في الصيغة الحالية لقوة تنافرها، كان من المفترض ضمنا أن الإلكترونات هي جسيمات نقطية. ومع ذلك، على المستوى الكمي، يتم “تلطخ” الإلكترون في الفضاء بحجم مميز يُعرف باسم طول موجة دي برولي الحراري. عندما يقترب إلكترونين “ملطخين” بنفس إسقاط الدوران من بعضهما البعض، فإن المسافة الفعالة بينهما لتفاعل كولوم تصبح أصغر من المسافة بين مراكز حزمهما الموجية، مما ينبغي أن يعزز التنافر. لتصحيح السلوك غير الصحيح، أدخل الفيزيائيون عاملًا تجريبيًا في صيغة القوة التنافرية بين هذه الإلكترونات. المعنى المادي للتعديل بسيط: يتم حساب القوة كما لو أن الإلكترونات ليست على مسافة r من بعضها البعض، ولكن على مسافة (r − αα)، حيث α هو طول موجتها الحرارية و α هي معلمة ترتيب مختارة للوحدة. هذا التصحيح، على الرغم من أنه لا يحتوي على اشتقاق نظري صارم، إلا أنه يقدم مساهمة كبيرة فقط في منطقة المسافات الصغيرة بترتيب ، حيث تكون التأثيرات الكمومية كبيرة. وعلى مسافات كبيرة، تتحول القوة بسلاسة إلى تنافر كولوم عادي. بمساعدتها، يتم استقرار سلوك النظام. في الشكل، يتم تصوير الحزم الموجية لإلكترونين لهما نفس إسقاط الدوران بشكل تخطيطي على أنها مناطق ذات حجم محدود، وهذه الدقة هي التي تجعلنا نتوقع تنافرًا أقوى من تنافر الجسيمات “النقطية” على مسافات قصيرة / © فيزياء البلازما تأثير طريقة تعديل القوة على تكوين العناقيد. تم الحصول على المنحنى الأسود دون مراعاة الطول الموجي المحدود للإلكترون ويظهر ذروة شاذة – وهي علامة على تكوين حالة غير مادية مرتبطة وانتهاك لمبدأ باولي. المنحنى الأحمر الذي تم الحصول عليه مع تعديل القوى ليس له مثل هذه الذروة، مما يشير إلى منع التجمع / © فيزياء البلازما وأوضح جورجي ديميانوف، كبير المحاضرين في قسم فيزياء العمليات ذات درجة الحرارة العالية في MIPT، جوهر النهج الجديد: “يمكننا القول أن نهجنا يأخذ في الاعتبار الحجم المميز للإلكترونات أثناء تفاعلها. ونتيجة لذلك، يبدأ التنافر الفعال بين الإلكترونات بنفس إسقاط الدوران في الظهور على مسافات أكبر. بين الجسيمات عما كانت عليه في النموذج الأصلي، ويمنع هذا التعديل تكوين حالات مرتبطة اصطناعية تنتهك مبدأ باولي في المحاكاة. أجرى العلماء سلسلة واسعة النطاق من حسابات الديناميكيات الجزيئية باستخدام إمكانات كاذبة معدلة لبلازما الهيدروجين غير المتحللة مع معامل انحطاط صغير. كما هو متوقع، عند درجات حرارة أعلى من 50 ألف درجة، تبين أن الاختلافات بين طريقة مونت كارلو الكمومية المقترحة والأكثر دقة كانت ضئيلة – أقل من 1٪ للطاقة والضغط. ومع ذلك، مع انخفاض درجة الحرارة، زادت التناقضات، وفي أبرد نقطة تم أخذها بعين الاعتبار، تبين أن الطاقة من نموذج الديناميكيات الجزيئية تم الاستهانة بها بعشرات بالمائة. في الوقت نفسه، أظهر تحليل وظائف التوزيع الشعاعي بوضوح أن هذا لم يعد بسبب تكوين مجموعات عملاقة – وكانت ذروتها المميزة غائبة عن الرسوم البيانية. المشكلة، وفقًا للمؤلفين، هي أنه على الرغم من أن إمكانات كيلبج الكاذبة المحسّنة تقارب بشكل صحيح طاقة التفاعل، إلا أن القوة المقابلة يمكن أن تنحرف بشكل ملحوظ عن القوة الدقيقة، الأمر الذي يؤدي في درجات الحرارة المنخفضة إلى زيادة الجذب بين الإلكترونات والبروتونات، ونتيجة لذلك، إلى المبالغة في تقدير تكوين الجزيئات في عمليات المحاكاة. اعتماد الطاقة الإجمالية لبلازما الهيدروجين على نسبة الكتلة mp/me في محاكاة الديناميكيات الجزيئية في ظل شرطين: (يسار) بدون جزيئات و (يمين) مع جزيئات / © فيزياء البلازما علق بافيل ليفاشوف، رئيس قسم فيزياء عمليات درجات الحرارة المرتفعة في MIPT، على هذه النتيجة: “كشف التحليل عن قيود منهجية مهمة للنموذج المستخدم. أتاح تعديل تفاعل القوة القضاء على النهج غير المادي للإلكترونات مع نفس إسقاط الدوران، لكن المزيد من المقارنة مع الحسابات الإحصائية الكمومية أظهرت أنه من أجل إعادة إنتاج خصائص النظام بشكل أكثر دقة، من الضروري توضيح شكل زوج الإمكانات الفعال، أي تدرجه. يعمل التعديل المقترح للقوى على استقرار المحاكاة ويجعل من الممكن دراسة بنية بلازما الهيدروجين في نظام درجات الحرارة المنخفضة الذي لم يكن من الممكن الوصول إليه سابقًا؛ علاوة على ذلك، فإن هذه الطريقة أقل تكلفة حسابيًا بشكل ملحوظ من طريقة مونت كارلو الكمومية، ويمكن استخدامها في الحسابات البارامترية واسعة النطاق لمعادلة حالة الهيدروجين وبلازما الديوتيريوم، وهو أمر مهم لمشاكل الفيزياء الفلكية والاندماج النووي الحراري بالقصور الذاتي، بينما يهدف المزيد من البحث إلى اختبار الفرضية حول اختلاف التدرج الكاذب مع التفاعل الدقيق وتعميم التعديل المقترح للقوى على أنظمة البلازما متعددة المكونات الأكثر تعقيدًا، بما في ذلك عمليات المحاكاة مع عدد متغير من الجسيمات لتحديد أكثر موثوقية للحد الديناميكي الحراري.