2025-11-13T11:31:10.165583

Localized fermions on the triangular lattice with Ising-like interactions

Regeciová, Kapcia
The model of localized fermions on the triangular lattice is analyzed in means of the Monte Carlo simulations in the grand canonical ensemble. The Hamiltonian of the system has a form of the extended Hubbard model (at the atomic limit) with nearest-neighbor Ising-like magnetic $J$ interactions and onsite Coulomb $U$ interactions. The model is investigated for both signs of $J$, arbitrary $U$ interaction and arbitrary chemical potential $μ$ (or, equivalently, arbitrary particle concentration $n$). Based on the specific heat capacity and sublattice magnetization analyses, the phase diagrams of the model are determined. For ferromagnetic case ($J<0$), the transition from the ordered phase (which is a standard ferromagnet and can be stable up to $k_{B}T/|J| \approx 0.61$) is found to be second-order (for sufficiently large temperatures $k_{B}T/|J| \gtrsim 0.2$) or first-order (for $-1<U/|J|<-0.65$ at the half-filling, i.e., $n=1$). In the case of $J>0$, the ordered phase occurs in a range of $-1/2<U/|J|<0$ (for $n=1$), while for larger $U$ the state with short-range order is also found (also for $n \neq 1$). The ordered phase is characterized by an antiferromagnetic arrangement of magnetic moments in two sublattices forming the hexagonal lattice. The transition from this ordered phase, which is found also for $μ\neq 0$ ($n \neq 1$) and $U/|J|>-1/2$ is always second-order for any model parameters. The ordered phase for $J>0$ can be stable up to $k_{B}T/|J| \approx 0.06$.
academic

الفرميونات المحلية على الشبكة الثلاثية مع تفاعلات تشبه إيزينج

المعلومات الأساسية

  • معرّف الورقة: 2412.10808
  • العنوان: الفرميونات المحلية على الشبكة الثلاثية مع تفاعلات تشبه إيزينج
  • المؤلفون: Lubomíra Regeciová (أكاديمية سلوفاكيا للعلوم)، Konrad Jerzy Kapcia (جامعة آدم ميتسكيفيتش)
  • التصنيف: cond-mat.stat-mech cond-mat.other cond-mat.quant-gas cond-mat.str-el physics.comp-ph
  • تاريخ النشر: 17 أكتوبر 2025
  • رابط الورقة: https://arxiv.org/abs/2412.10808

الملخص

تحلل هذه الدراسة نموذج الفرميونات المحلية على الشبكة الثلاثية من خلال محاكاة مونت كارلو في المجموعة الكنسية. يتخذ هاميلتونيان النظام شكل نموذج هوبارد الممتد (الحد الذري)، يتضمن تفاعلات إيزينج المغناطيسية بين أقرب الجيران J وتفاعلات كولوم الموضعية U. تفحص الدراسة كلا الإشارتين لـ J وأي تفاعل U وأي جهد كيميائي μ (أو بشكل مكافئ، أي كثافة جزيئات n). بناءً على تحليل السعة الحرارية والمغناطيسية تحت الشبكة، تم تحديد مخطط الطور للنموذج. بالنسبة للحالة الحديدية المغناطيسية (J<0)، يمكن أن يكون الانتقال من الطور المرتب انتقالاً من الدرجة الثانية (عند درجات حرارة عالية بما يكفي kBT/|J|≳0.2) أو انتقالاً من الدرجة الأولى (عند الملء النصفي -1<U/|J|<-0.65). بالنسبة لحالة J>0، يظهر الطور المرتب في النطاق -1/2<U/|J|<0 (عند n=1)، وبالنسبة لقيم U الأكبر، تم اكتشاف حالات ذات ترتيب قصير المدى أيضاً.

خلفية البحث والدافع

  1. المشكلة الأساسية: دراسة سلوك انتقالات الطور المغناطيسية للفرميونات المحلية في الأنظمة الهندسية المحبطة، خاصة نموذج هوبارد الممتد على الشبكة الثلاثية في الحد الذري.
  2. الأهمية:
    • تتمتع الشبكة الثلاثية بخصائص إحباط هندسي متأصلة، مما يؤدي إلى عدم إمكانية تلبية جميع التفاعلات المضادة للحديد في نفس الوقت
    • يمثل هذا النموذج امتداداً غير تافه للنموذج الكلاسيكي لإيزينج إلى الأنظمة الفرميونية
    • له أهمية كبيرة في فهم الأنظمة الإلكترونية المترابطة بقوة والمواد المغناطيسية
  3. قيود الطرق الموجودة:
    • تم دراسة الشبكات فوق المكعبة بشكل كافٍ، لكن النتائج التحليلية الدقيقة للشبكة الثلاثية ناقصة
    • تبالغ التقريبات متوسطة المجال في درجات حرارة الانتقال
    • يفتقد الأدب إلى دراسات عددية مفصلة لتأثيرات الإحباط الهندسي
  4. الدافع للبحث:
    • سد الفجوة في الدراسات العددية لهذا النموذج على الشبكة الثلاثية
    • استكشاف تأثير الإحباط الهندسي على سلوك الانتقالات
    • توفير أساس نظري للمواد التجريبية (مثل Li2MnTeO6 و KCeO2 وغيرها)

المساهمات الأساسية

  1. أول دراسة منهجية لمونت كارلو: محاكاة عددية شاملة لنموذج هوبارد الممتد مع تفاعلات إيزينج المغناطيسية على الشبكة الثلاثية
  2. بناء مخطط طور كامل: تحديد مخطط الطور ثلاثي الأبعاد (درجة الحرارة-الجهد الكيميائي-قوة التفاعل) للحالات الحديدية المغناطيسية (J<0) والمضادة للحديد (J>0)
  3. تحليل خصائص الانتقال: تحديد دقيق لحدود الانتقالات من الدرجة الأولى والثانية، مع اكتشاف أن جميع الانتقالات في الحالة المضادة للحديد من الدرجة الثانية
  4. كشف تأثيرات الإحباط: اكتشاف أن الإحباط الهندسي يؤثر بشكل كبير على بنية مخطط الطور، حيث تكون درجات حرارة استقرار الطور المرتب في الحالة المضادة للحديد أقل بمرتبة من الحالة الحديدية المغناطيسية
  5. تحديد حالات الترتيب قصيرة المدى: اكتشاف وتوصيف منطقة وجود حالات الترتيب قصيرة المدى في الحالة المضادة للحديد

شرح الطريقة

تعريف المهمة

دراسة سلوك الانتقالات في نظام فرميوني محلي على الشبكة الثلاثية، حيث المدخلات هي معاملات النموذج (J, U, μ, T) والمخرجات هي مخطط الطور وخصائص الانتقالات.

بنية النموذج

هاميلتونيان: H^=U∑in^i↑n^i↓+2Jz∑⟨i,j⟩s^is^j−(μ+U2)∑in^i\hat{H} = U\sum_i \hat{n}_{i\uparrow}\hat{n}_{i\downarrow} + \frac{2J}{z}\sum_{\langle i,j\rangle} \hat{s}_i\hat{s}_j - \left(\mu + \frac{U}{2}\right)\sum_i \hat{n}_i

حيث:

  • s^i=(n^i↑−n^i↓)/2\hat{s}_i = (\hat{n}_{i\uparrow} - \hat{n}_{i\downarrow})/2: عامل الدوران المحلي
  • n^i=n^i↑+n^i↓\hat{n}_i = \hat{n}_{i\uparrow} + \hat{n}_{i\downarrow}: عامل عدد الجزيئات
  • z=6z = 6: رقم التنسيق للشبكة الثلاثية

حالات الموقع: يمكن لكل موقع أن يكون في أربع حالات:

  • ∣0⟩|0\rangle: موقع فارغ
  • ∣↑⟩|\uparrow\rangle: احتلال فردي (دوران لأعلى)
  • ∣↓⟩|\downarrow\rangle: احتلال فردي (دوران لأسفل)
  • ∣↑↓⟩|\uparrow\downarrow\rangle: احتلال مزدوج

نقاط الابتكار التقني

  1. محاكاة المجموعة الكنسية: استخدام المجموعة الكنسية مع تقلبات عدد الجزيئات لوصف أكثر دقة للنظام المتحكم به بواسطة الجهد الكيميائي
  2. تحليل الحجم المحدود: دراسة منهجية لتأثيرات الحجم المحدود من L = 12×12 إلى 36×36
  3. طرق تحديد الانتقال:
    • استخدام طريقة دالة توزيع الطاقة P(E) من قبل Challa وآخرين للتمييز بين الانتقالات من الدرجة الأولى والثانية
    • الجمع بين تحليل السعة الحرارية والمغناطيسية تحت الشبكة مع سلوك التحجيم
  4. تحليل تحت الشبكة: تقسيم الشبكة الثلاثية إلى ثلاث شبكات فرعية متكافئة لتحليل دقيق لبنية الترتيب المغناطيسي

الإعدادات التجريبية

مجموعة البيانات

  • أحجام النظام: مجموعات مثلثة، L = 12×12, 24×24, 36×36 موقع
  • نطاق المعاملات:
    • درجة الحرارة: kBT/|J| = 0.005 - 1.0
    • التفاعل: U/|J| = -1.0 - 5.0
    • الجهد الكيميائي: μ/|J| = -0.6 - 0 (باستخدام تماثل الجزيء-الفراغ)

مؤشرات التقييم

  1. السعة الحرارية: C=⟨ε2⟩−⟨ε⟩2/T2C = \langle\varepsilon^2\rangle - \langle\varepsilon\rangle^2/T^2
  2. المغناطيسية تحت الشبكة: mα=(6/L)∑i∈α⟨s^i⟩m_\alpha = (6/L)\sum_{i\in\alpha}\langle\hat{s}_i\rangle
  3. كثافة الجزيئات: n=⟨n^i⟩n = \langle\hat{n}_i\rangle
  4. دالة توزيع الطاقة: P(E)P(E) لتحديد خصائص الانتقال

تفاصيل التنفيذ

  • الخوارزمية: خوارزمية مونت كارلو من نوع Metropolis
  • خطوات التوازن: 10610^6 خطوة
  • خطوات الإحصاء: 10610^6 خطوة
  • مسح درجة الحرارة: تبريد تدريجي من درجات حرارة عالية إلى منخفضة
  • شروط الحدود: شروط حدود دورية

النتائج التجريبية

النتائج الرئيسية

الحالة الحديدية المغناطيسية (J < 0):

  • أعلى درجة حرارة انتقال: kBT/|J| ≈ 0.605 (عند U → ∞، قريبة من 0.607 لنموذج إيزينج الكلاسيكي)
  • موقع النقطة الثلاثية الحرجة: kBT/|J| ≈ 0.18-0.20، مستقل تقريباً عن U و μ
  • خصائص الانتقال: انتقال من الدرجة الأولى عند درجات حرارة منخفضة، انتقال من الدرجة الثانية عند درجات حرارة عالية
  • الطور المرتب: طور حديدي مغناطيسي قياسي، جميع المغناطيسيات تحت الشبكة متساوية

الحالة المضادة للحديد (J > 0):

  • أعلى درجة حرارة انتقال: kBT/|J| ≈ 0.06 (أقل بكثير من الحالة الحديدية المغناطيسية)
  • نطاق وجود الطور المرتب: -1/2 < U/|J| < 0 (عند الملء النصفي)
  • خصائص الانتقال: جميع الانتقالات من الدرجة الثانية
  • بنية الطور المرتب: شبكتان فرعيتان تشكلان ترتيباً مضاداً للحديد، المغناطيسية تحت الشبكة الثالثة تساوي صفراً

التجارب الاستئصالية

التحقق من سلوك الحد الديناميكي الحراري للانتقال من خلال تغيير حجم النظام:

  • الانتقال من الدرجة الأولى: ظهور بنية ثنائية القمة في توزيع الطاقة
  • الانتقال من الدرجة الثانية: ظهور شذوذ λ في السعة الحرارية، مع نمو القمة مع حجم النظام
  • الترتيب قصير المدى: ظهور قمة عريضة في السعة الحرارية لا تتحجم مع حجم النظام

دراسات الحالات

سلوكيات الانتقال النموذجية:

  1. U = -0.7, μ = -0.10: انتقال F-NO من الدرجة الأولى، مع قفزات غير متصلة في المغناطيسية وعدد الجزيئات
  2. U = -0.6, μ = -0.12: انتقال F-NO من الدرجة الثانية، مع تغير متصل في معامل الترتيب
  3. U = -0.2, μ = -0.15: انتقال AF-NO من الدرجة الثانية، مع ترتيب مضاد للحديد لشبكتين فرعيتين

الاكتشافات التجريبية

  1. التأثير القوي للإحباط الهندسي: درجات حرارة الانتقال في الحالة المضادة للحديد أقل بمرتبة من الحالة الحديدية المغناطيسية
  2. عدم التماثل في خصائص الانتقال: الحالة الحديدية المغناطيسية تحتوي على انتقالات من الدرجة الأولى والثانية، بينما الحالة المضادة للحديد تحتوي على انتقالات من الدرجة الثانية فقط
  3. الانتشار الواسع للترتيب قصير المدى: ملاحظة واسعة لحالات الترتيب قصيرة المدى في الحالة المضادة للحديد

الأعمال ذات الصلة

  1. نموذج إيزينج الكلاسيكي: الحل الدقيق على الشبكة الثلاثية (Houtappel, 1950) يوفر معياراً لهذا البحث
  2. نموذج هوبارد الممتد: الدراسات التباينية والعددية على الشبكات فوق المكعبة توفر أساساً منهجياً
  3. أنظمة المغناطيسية المحبطة: الدراسات النظرية والتجريبية للمواد المضادة للحديد على الشبكة الثلاثية توفر خلفية لتفسير النتائج
  4. تجارب الذرات الباردة: تحقيق غازات الذرات الباردة في الشبكات البصرية يوفر إمكانية للتحقق التجريبي من النموذج

الاستنتاجات والمناقشة

الاستنتاجات الرئيسية

  1. يؤثر الإحباط الهندسي للشبكة الثلاثية بشكل كبير على سلوك الانتقالات في نظام الفرميونات المحلية
  2. التفاعلات الحديدية المغناطيسية والمضادة للحديد تؤدي إلى بنى مخطط طور وخصائص انتقال مختلفة تماماً
  3. الطور المرتب في الحالة المضادة للحديد يتمتع ببنية ترتيب مضادة للحديد فريدة لشبكتين فرعيتين
  4. جميع الانتقالات المضادة للحديد هي انتقالات متصلة، بينما توجد تحويلات من الدرجة الأولى إلى الثانية في الحالة الحديدية المغناطيسية

القيود

  1. قيود حجم النظام: بسبب قيود الموارد الحسابية، أكبر نظام يبلغ 36×36 فقط
  2. تقريب الحد الذري: تجاهل حدود قفز الجزيئات قد يؤثر على التنبؤات الكمية للمواد الفعلية
  3. تقريب الدوران الكلاسيكي: قد لا يكون دقيقاً كافياً في مناطق التقلبات الكمية القوية

الاتجاهات المستقبلية

  1. دراسات تتضمن تفاعلات ذات مدى أطول
  2. تطبيق طرق مونت كارلو الكمية
  3. مقارنات كمية مع المواد الفعلية
  4. التحقق التجريبي في أنظمة الذرات الباردة

التقييم المتعمق

المزايا

  1. صرامة منهجية: استخدام طرق مونت كارلو الناضجة وتقنيات تحديد الانتقالات
  2. شمولية النتائج: بناء مخطط طور ثلاثي الأبعاد كامل يغطي مناطق معاملات متعددة
  3. عمق الرؤى الفيزيائية: كشف التأثير الجوهري للإحباط الهندسي على سلوك الانتقالات
  4. قوة الارتباط التجريبي: ارتباط وثيق مع مواد فعلية متعددة وأنظمة ذرات باردة

أوجه القصور

  1. نقص التحليل النظري: افتقار إلى شرح نظري عميق لآليات الانتقال
  2. غياب خصائص الديناميكا: عدم تناول الخصائص الديناميكية والنقل في النظام
  3. تجاهل التأثيرات الكمية: تجاهل تأثيرات التقلبات الكمية المهمة في الحد الذري

القيمة التأثيرية

  1. القيمة الأكاديمية: توفير معايير عددية مهمة لنظرية الأنظمة المغناطيسية المحبطة
  2. القيمة العملية: توفير إرشادات لفهم وتصميم المواد ذات الصلة
  3. قابلية إعادة الإنتاج: وصف منهجي تفصيلي وتتمتع النتائج بقابلية إعادة إنتاج جيدة

السيناريوهات المطبقة

  1. النمذجة النظرية للمواد المغناطيسية على الشبكة الثلاثية
  2. الإرشادات النظرية لتجارب الشبكات البصرية للذرات الباردة
  3. دراسة الانتقالات في الأنظمة الإلكترونية المترابطة بقوة
  4. اختبارات المعايير لأنظمة الكم متعددة الأجسام المحبطة

المراجع

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