DOI: https://doi.org/10.1038/s41535-026-00908-1
تاريخ النشر: 2026-06-06
المؤلف: Yang Yang وآخرون
الموضوع الرئيسي: فيزياء المادة المكثفة المتقدمة
نظرة عامة
تبحث هذه الدراسة في ديناميات الاسترخاء غير المتوازنة للسوائل الكمومية المغزلية، مع التركيز بشكل خاص على نموذج كيتاف بعد التبريد الحراري. تؤكد الدراسة على دور الفيزونات، وهي إثارات تدفق $Z_2$ المقطوعة في مجال القياس الناشئ في سائل كيتاف المغزلي. يتميز الديناميات على المدى الطويل بالانتشار والانقراض الثنائي للفيزونات، المتأثرة بحاجز التنشيط لحركة الفيزون وتفاعلاتها التي تتوسطها فيرميونات مايورانا المجزأة.
باستخدام نهج مونت كارلو الحركي الفعال، يكشف المؤلفون عن مجموعة معقدة من الأنظمة الديناميكية المعتمدة على درجة الحرارة. مع انخفاض درجة الحرارة، ينتقل النظام من الانقراض المحدود بالانتشار إلى الديناميات المحدودة بسرعة نهائية، مما يؤدي في النهاية إلى اعتقال ديناميكي كبير وتجميد. يرتبط هذا التجميد بظهور بلورات فيزون $\sqrt{3} \times \sqrt{3}$ طويلة العمر، والتي تمنع المزيد من الانقراض. ومن الجدير بالذكر أنه حتى بدون ترتيب بلوري حقيقي طويل المدى، يظهر النظام ديناميات تكتل خفية للفيزونات المرتبطة بتماثل $Z_3$ الناشئ.
مقدمة
تناقش مقدمة الورقة نموذج كيتاف، الذي يعد إطارًا مهمًا لدراسة السوائل الكمومية المغزلية التي تتميز بإثارات مجزأة وهياكل قياس ناشئة. أدت القابلية الدقيقة لحل النموذج وزيادة عدد تجسيدات المواد إلى فهم شامل لخصائصه المتوازنة، بما في ذلك مراحل الحالة الأساسية والسلوك الديناميكي الحراري. ومع ذلك، لا تزال الديناميات غير المتوازنة، خاصة بعد التبريد الحراري، غير مستكشفة بشكل كاف. في سائل كيتاف المغزلي، تتجزأ المغازل إلى فيرميونات مايورانا بدون فجوة مرتبطة بمجال قياس \( \mathbb{Z}_2 \) الناشئ، حيث تمثل الفيزونات الإثارات الأساسية للقياس.
تسلط الورقة الضوء على تعقيد ديناميات الاسترخاء في الأنظمة التي تفتقر إلى كسر التماثل التقليدي، مثل السوائل الكمومية المغزلية، حيث يتم التحكم في الاسترخاء بواسطة ديناميات درجات الحرية الناشئة بدلاً من نمو المجالات. يقترح المؤلفون إطارًا نظريًا لدراسة ديناميات الاسترخاء في سائل كيتاف المغزلي، مستفيدين من القابلية الدقيقة لحل النموذج. يقومون بتطوير وصف مونت كارلو الحركي لديناميات الفيزون، كاشفين عن ثلاث مسارات استرخاء متميزة تتأثر بالانتشار المعتمد على الحرارة والتفاعلات التي تتوسطها الفيرميونات. تهدف هذه الدراسة إلى توضيح آليات الاسترخاء غير المتوازن وتقديم رؤى ذات صلة تجريبياً حول المواد المرشحة لنموذج كيتاف.
النتائج
تناقش قسم النتائج من ورقة البحث ديناميات نموذج كيتاف المتساوي، المميز بالهاميلتونيان \(\hat{H} = -J \sum_{\langle ij \rangle \parallel \alpha} \sigma_i^\alpha \sigma_j^\alpha\). يسمح النموذج بحل دقيق من خلال تجزئة درجات حرية المغزل إلى فيرميونات مايورانا، التي تتفاعل مع مجال قياس ناشئ \(Z_2\). تكشف الدراسة أن الفيزونات، الإثارات في هذا النموذج، تتصرف بشكل مشابه لمتغيرات إيسين الكلاسيكية في غياب الاضطرابات الخارجية. ومع ذلك، عندما يرتبط النظام ببيئة حرارية، تظهر ديناميات الفيزون، مما يؤدي إلى مشهد غني من عمليات الاسترخاء المميزة بثلاث مقاييس زمنية متميزة: زمن استرخاء فيرميونات مايورانا \(\tau_M\)، زمن تدهور الفيزون \(\tau_2\)، وزمن قفز الفيزون \(\tau_1\).
تشير النتائج إلى أنه عند درجات حرارة عالية، يتدهور كثافة الفيزون كـ \(\rho(t) \sim t^{-1} \ln t\)، متماشياً مع الانقراض المحدود بالانتشار. عند درجات حرارة متوسطة، تسرع التفاعلات الجذابة بين الفيزونات الانقراض، مما يؤدي إلى تدهور أسرع لـ \(\rho(t) \sim t^{-2}\). عند درجات حرارة منخفضة، تستقر كثافة الفيزون عند قيمة غير صفرية بسبب تكوين تكوينات بلورية فيزون غير مستقرة، مما يثبط المزيد من الانقراض. تبرز الدراسة أيضًا ظهور بلورة فيزون $\sqrt{3} \times \sqrt{3}$، والتي هي حالة غير مستقرة بدلاً من مرحلة توازن، وتناقش التكتل البطيء للفيزونات العملاقة، مما يشير إلى تفاعلات معقدة تنحرف عن النماذج التقليدية. بشكل عام، توفر الدراسة رؤى حول الديناميات غير المتوازنة للفيزونات في نموذج كيتاف، مع التأكيد على التفاعل بين التقلبات الحرارية وتفاعلات الفيزون.
مناقشة
في هذا القسم، يعترف المؤلفون بمساهمات من شانغ-شون زانغ وناتاليا ب. بيركنز، ويستعرضون النتائج من محاكياتهم حول ديناميات كثافة الفيزون بعد التبريد الحراري. تشير النتائج إلى تدهور كثافة الفيزون المعتمد على درجة الحرارة، $\rho$، المميز بسلوكيات متميزة عند درجات حرارة مختلفة بعد التبريد ($T/|J| = 10^{-1}, 10^{-2}, 10^{-3}$). يتم تمثيل ديناميات التدهور من خلال نتائج المحاكاة المتوسطة، مع أشكال مقياس مرجعية تشير إلى تدهور أسي، تدهور بقوة، وسلوك لوغاريتمي. ومن الجدير بالذكر أن التدهور يظهر سلوكاً $\sim 1/t^2$ عند درجات حرارة منخفضة، كما هو موضح بواسطة الخط الأخضر المتقطع.
علاوة على ذلك، يستكشف المؤلفون تكوين بلورات الفيزون والعملاق عبر الشبكات الفرعية، محددين ثلاث مجالات متميزة بناءً على مواقع الفيزونات الأقرب. توضح الدراسة تسلسلًا هرميًا لمقاييس الطاقة المنخفضة التي تحكم ديناميات الفيزون، تمتد من مرحلة بارامغنت كيتاف عند درجات حرارة عالية إلى عملاق بلوري مجمد عند درجات حرارة منخفضة. درجة الحرارة الانتقالية للديناميات المجمدة تقريبًا $T_{\text{freezing}}/|J| \sim 0.006$، مع مقاييس طاقة التفاعل وخلق الأزواج المشار إليها كـ $T_{\text{int}} \sim 0.04$ و $T_{\text{pair}} \sim 0.26$، على التوالي. توضح هذه النتائج الأنظمة الديناميكية المعقدة للفيزونات، بما في ذلك ديناميات الانتشار والانقراض والديناميات الموقوفة-الهاربة، مما يساهم في فهم أعمق للأنظمة الكمومية المغزلية.
DOI: https://doi.org/10.1038/s41535-026-00908-1
Publication Date: 2026-06-06
Author(s): Yang Yang et al.
Primary Topic: Advanced Condensed Matter Physics
Overview
This research investigates the nonequilibrium relaxation dynamics of quantum spin liquids, specifically focusing on the Kitaev honeycomb model after a thermal quench. The study emphasizes the role of visons, which are gapped $Z_2$ flux excitations of the emergent gauge field in the Kitaev spin liquid. The long-time dynamics are characterized by the diffusion and pair annihilation of visons, influenced by the activation barrier for vison motion and their interactions mediated by fractionalized Majorana fermions.
Utilizing an effective kinetic Monte Carlo approach, the authors reveal a complex set of temperature-dependent dynamical regimes. As the temperature decreases, the system transitions from diffusion-limited annihilation to terminal-velocity-limited kinetics, eventually leading to significant dynamical arrest and freezing. This freezing is associated with the emergence of long-lived, metastable $\sqrt{3} \times \sqrt{3}$ vison crystals, which inhibit further annihilation. Notably, even without true long-range crystalline order, the system displays hidden coarsening dynamics of vison superclusters linked to an emergent $Z_3$ symmetry.
Introduction
The introduction of the paper discusses the Kitaev honeycomb model, which serves as a significant framework for investigating quantum spin liquids characterized by fractionalized excitations and emergent gauge structures. The model’s exact solvability and the increasing number of material realizations have led to a comprehensive understanding of its equilibrium properties, including ground-state phases and thermodynamic behavior. However, the nonequilibrium dynamics, particularly following thermal quenches, remain underexplored. In the Kitaev spin liquid, spins fractionalize into gapless Majorana fermions coupled with an emergent \( \mathbb{Z}_2 \) gauge field, with visons representing the elementary gauge excitations.
The paper highlights the complexity of relaxation dynamics in systems lacking conventional symmetry breaking, such as quantum spin liquids, where relaxation is governed by the dynamics of emergent degrees of freedom rather than domain growth. The authors propose a theoretical framework to study the relaxation dynamics of the Kitaev spin liquid, leveraging the model’s exact solvability. They develop a kinetic Monte Carlo description of vison dynamics, revealing three distinct relaxation pathways influenced by thermally activated diffusion and fermion-mediated interactions. This work aims to elucidate the mechanisms of nonequilibrium relaxation and provide experimentally relevant insights into candidate Kitaev materials.
Results
The results section of the research paper discusses the dynamics of the isotropic Kitaev honeycomb model, characterized by the Hamiltonian \(\hat{H} = -J \sum_{\langle ij \rangle \parallel \alpha} \sigma_i^\alpha \sigma_j^\alpha\). The model allows for an exact solution through the fractionalization of spin degrees of freedom into Majorana fermions, which interact with an emergent \(Z_2\) gauge field. The study reveals that visons, the excitations in this model, behave similarly to classical Ising variables in the absence of external perturbations. However, when the system is coupled to a thermal environment, vison dynamics emerge, leading to a rich landscape of relaxation processes characterized by three distinct time scales: the relaxation time of Majorana fermions \(\tau_M\), the vison decoherence time \(\tau_2\), and the vison hopping time \(\tau_1\).
The findings indicate that at high temperatures, the vison density decays as \(\rho(t) \sim t^{-1} \ln t\), consistent with diffusion-limited annihilation. At intermediate temperatures, attractive interactions between visons accelerate annihilation, resulting in a faster decay of \(\rho(t) \sim t^{-2}\). At low temperatures, the vison density stabilizes at a nonzero value due to the formation of metastable vison crystal configurations, which suppress further annihilation. The study also highlights the emergence of a \(\sqrt{3} \times \sqrt{3}\) vison crystal, which is a metastable state rather than an equilibrium phase, and discusses the slow coarsening of vison superclusters, suggesting complex interactions that deviate from conventional models. Overall, the research provides insights into the non-equilibrium dynamics of visons in the Kitaev model, emphasizing the interplay between thermal fluctuations and vison interactions.
Discussion
In this section, the authors acknowledge contributions from Shang-Shun Zhang and Natalia B. Perkins, and detail the findings from their simulations on vison density dynamics following a thermal quench. The results indicate a temperature-dependent decay of vison density, $\rho$, characterized by distinct behaviors at various post-quench temperatures ($T/|J| = 10^{-1}, 10^{-2}, 10^{-3}$). The decay dynamics are represented through averaged simulation results, with reference scaling forms indicating exponential decay, power-law decay, and logarithmic behavior. Notably, the decay exhibits a $\sim 1/t^2$ behavior at lower temperatures, as highlighted by the green dashed line.
Furthermore, the authors explore the formation of vison crystals and super-clusters across sublattices, identifying three distinct domains based on the nearest vison locations. The study delineates a hierarchy of low-energy scales that govern vison dynamics, spanning from the Kitaev paramagnet phase at high temperatures to frozen-crystal super-clusters at low temperatures. The crossover temperature for the transition to frozen-crystal dynamics is approximately $T_{\text{freezing}}/|J| \sim 0.006$, with interaction and pair-creation energy scales indicated as $T_{\text{int}} \sim 0.04$ and $T_{\text{pair}} \sim 0.26$, respectively. These findings elucidate the complex dynamical regimes of visons, including diffusion-annihilation and arresting-escaping dynamics, contributing to a deeper understanding of quantum spin systems.