دراسة التشخيص الدقيق لإحصائيات مستويات الطاقة في البوزونات المتفاعلة المحصورة بشكل توافقي Exact diagonalization study of energy-level statistics in harmonically confined interacting bosons

شارك:
المجلة: Physical review. E، المجلد: 113، العدد: 1
DOI: https://doi.org/10.1103/7243-bpjp
PMID: https://pubmed.ncbi.nlm.nih.gov/41715869
تاريخ النشر: 2026-01-08
المؤلف: Zhenyun Du وآخرون
الموضوع الرئيسي: الفوضى الكمومية والأنظمة الديناميكية

نظرة عامة

تبحث هذه الورقة البحثية في الخصائص الطيفية للبوزونات المحصورة في مستوى شبه ثنائي الأبعاد وتتفاعل من خلال جهد غاوسي طارد، باستخدام تقنيات التشخيص الدقيق. تركز الدراسة على أنظمة تحتوي على 12 و16 و20 بوزون، مع فحص نظامين من التفاعل: معتدل (حيث تكون طاقة التفاعل صغيرة مقارنة بطاقة الفخ) وقوي (حيث تكون متقاربة). تشمل التحليلات كل من الحالات غير الدوارة (Lz = 0) وحالات الدوامة الفردية الدوارة (Lz = N)، مع اعتبار الزخم الزاوي الأعلى (Lz = 2N وLz = 3N) فقط في نظام التفاعل القوي. تشمل الأدوات الإحصائية الرئيسية المستخدمة توزيع تباعد الجيران الأقرب (NNSD)، ونسب تباعد المستويات المتتالية (P(r))، وإحصائية دايسون-ميتا ∆3، وتباين عدد المستويات (Σ²(L)).

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

مقدمة

تناقش مقدمة هذه الورقة البحثية تطبيق نظرية المصفوفات العشوائية (RMT) لتحليل أطياف مستويات الطاقة للأنظمة الكمومية المعقدة، وخاصة في النوى الثقيلة. تم تقديمها في البداية من قبل E. P. Wigner وتطويرها لاحقًا بواسطة ميهتا ودايسون، توفر RMT إطارًا إحصائيًا لفهم أطياف الطاقة. جانب حاسم من هذا التحليل هو عملية “الفك”، التي تقوم بتطبيع كثافة المستوى المتوسطة إلى الواحد، مما يسمح باستخراج ميزات محددة للنظام مع تقليل الأخطاء من خلال التناسبات متعددة الحدود أو المتوسطات المتحركة. تبرز الورقة مقاييس إحصائية متنوعة، مثل توزيع تباعد الجيران الأقرب (NNSD) $P(s)$ وتوزيع نسبة تباعد المستويات المتتالية $P(r)$، والتي تعتبر ضرورية لاستكشاف الارتباطات القصيرة والطويلة المدى في مستويات الطاقة.

يشير المؤلفون إلى أنه بينما تتطلب NNSD وإحصائيات أخرى الفك وقد تقدم أخطاء، فإن توزيع نسبة تباعد المستويات المتتالية لا يتطلب ذلك، مما يجعله أداة أكثر موثوقية لتحليل الارتباطات القصيرة المدى. تشير المقدمة أيضًا إلى تخمينات مهمة في الفوضى الكمومية، وهي تخمينات بوهغاس-جيانوني-شميت (BGS) وتخمين بيري-تابور، التي تربط الخصائص الإحصائية لمستويات الطاقة بالسلوك الكلاسيكي للأنظمة. تهدف الورقة إلى التحقيق في الخصائص الإحصائية لمستويات الطاقة في البوزونات المتفاعلة المحصورة تحت قوى تفاعل متغيرة وتأثيرات دورانية، وهو مجال لم يتم استكشافه بشكل شامل. تعد الدراسة بتقديم رؤى قيمة حول التفاعل بين طاقة التفاعل وطاقة الفخ، مع آثار لفهم الأنظمة الكمومية المعقدة وسلوكها الفوضوي.

النتائج

في هذا القسم، يبحث المؤلفون في الخصائص الطيفية لنظام من $N = 12، 16، 20$ ذرات بوزونية من $^{87}\text{Rb}$ محصورة في فخ هارموني شبه ثنائي الأبعاد، يتميز بمعامل أنيسوتروبي $\lambda_z = 4$ وتردد محوري $\omega_z = 2\pi \times 220 \text{ Hz}$. يتم نمذجة التفاعل بين الذرات باستخدام جهد غاوسي طارد مع مدى محدد بواسطة $\sigma = 0.1 a_\perp$، حيث $a_\perp = 1.446 \, \mu\text{m}$ هو طول الفخ. تستخدم الدراسة التشخيص الدقيق لتحليل مستويات الطاقة، مع التركيز على نظامين من التفاعل: معتدل ($g_2 = 0.3669$) وقوي ($g_2 = 3.669$)، تم تحقيقه من خلال تغيير طول تشتت الموجة s $a_{sc}$.

تكشف النتائج عن اعتماد إحصائيات مستويات الطاقة على عدد البوزونات وقوة التفاعل. تشير معلمات برودي الملائمة $b$ للحالة غير الدوارة إلى انتقال من إحصائيات بواسون (لـ $b = 0$) إلى إحصائيات مجموعة غاوسي متعامدة (GOE) (لـ $b = 1$). على وجه التحديد، بالنسبة للتفاعلات المعتدلة، تكون معلمات برودي $b = 0.33$ لـ $N = 12$، و$b = 0.16$ لـ $N = 16$، و$b = 0.04$ لـ $N = 20$. في نظام التفاعل القوي، تكون المعلمات $b = 0.53$، و$b = 0.85$، و$b = 0.51$، على التوالي. تشير هذه النتائج إلى أنه مع زيادة عدد البوزونات، تتغير إحصائيات المستويات، مما يعكس تأثير التفاعلات الثنائية وأبعاد النظام على الارتباطات الطيفية.

المناقشة

في هذا القسم، يقدم المؤلفون تحليلًا مفصلًا لنظام من البوزونات المتفاعلة المحصورة في جهد هارموني شبه ثنائي الأبعاد وتتعرض للدوران. يتم صياغة هاملتونيان العديد من الأجسام في إطار دوار مشترك، مع تضمين الطاقة الحركية، وجهد خارجي محصور، وتفاعل ثنائي يتم نمذجته بواسطة جهد غاوسي. يتم اشتقاق قوة التفاعل الفعالة في بعدين، مما يكشف أن النظام يعمل في نظام طارد مع طول تشتت $a_{sc} > 0$. يتم الحصول على الطيف الذاتي بطريقة متغيرة، مما يؤدي إلى مشكلة قيمة ذاتية مصفوفة يتم حلها لاستخراج مستويات الطاقة.

تستخدم التحليلات الإحصائية لمستويات الطاقة أدوات متنوعة للتحقيق في كل من الارتباطات القصيرة والطويلة المدى. يتم استخدام توزيع تباعد الجيران الأقرب (NNSD) وتوزيع نسبة تباعد المستويات المتتالية لتقييم الارتباطات الطيفية. تشير النتائج إلى أنه في نظام التفاعل المعتدل، مع زيادة عدد البوزونات، يقترب NNSD من توزيع بواسون، مما يشير إلى عدم وجود طرد للمستويات ويشير إلى سلوك منتظم. على العكس، في نظام التفاعل القوي، يتماشى NNSD بشكل أقرب مع توزيع مجموعة غاوسي متعامدة (GOE)، مما يشير إلى سلوك فوضوي يتميز بطرد قوي للمستويات. تُستخدم إحصائية دايسون-ميتا $\Delta_3(L)$ وتباين عدد المستويات $\Sigma^2(L)$ أيضًا لتحليل الارتباطات الطويلة المدى، مما يكشف عن انتقال من سلوك قابل للتكامل إلى غير قابل للتكامل مع تغير قوة التفاعل وعدد البوزونات. بشكل عام، تسلط النتائج الضوء على التفاعل المعقد بين قوة التفاعل، وعدد البوزونات، والارتباطات الطيفية في هذا النظام الكمومي متعدد الأجسام.

Journal: Physical review. E, Volume: 113, Issue: 1
DOI: https://doi.org/10.1103/7243-bpjp
PMID: https://pubmed.ncbi.nlm.nih.gov/41715869
Publication Date: 2026-01-08
Author(s): Zhenyun Du et al.
Primary Topic: Quantum chaos and dynamical systems

Overview

This research paper investigates the spectral properties of bosons confined in a quasi-2D plane and interacting through a repulsive Gaussian potential, utilizing exact diagonalization techniques. The study focuses on systems with 12, 16, and 20 bosons, examining two interaction regimes: moderate (where interaction energy is small compared to trap energy) and strong (where they are comparable). The analysis includes both non-rotating states (Lz = 0) and rotating single-vortex states (Lz = N), with higher angular momenta (Lz = 2N and Lz = 3N) considered only in the strong interaction regime. Key statistical tools employed include the nearest-neighbor spacing distribution (NNSD), ratios of consecutive level spacings (P(r)), the Dyson-Mehta ∆3 statistic, and level number variance (Σ²(L)).

The findings reveal that in the moderate interaction regime, the non-rotating system exhibits Poisson statistics, indicative of regular energy spectra, while in the strong interaction regime, the system transitions to chaotic behavior characterized by Gaussian Orthogonal Ensemble (GOE) statistics. The rotating single-vortex state in the moderate regime shows weak chaos, whereas in the strong regime, it displays strong chaotic behavior, with rotation enhancing the chaotic characteristics across both interaction regimes. The study underscores the influence of interaction and trap energies on the transition from regular to chaotic behavior in quantum systems, suggesting the universal applicability of random matrix theory (RMT) in understanding spectral correlations in ultracold Bose systems. Future work is proposed to extend these findings to explore spectral form factors and the dynamics of energy-level statistics in more complex models.

Introduction

The introduction of this research paper discusses the application of Random Matrix Theory (RMT) to analyze the energy-level spectra of complex quantum many-body systems, particularly in heavy nuclei. Initially introduced by E. P. Wigner and further developed by Mehta and Dyson, RMT provides a statistical framework for understanding energy spectra. A critical aspect of this analysis is the “unfolding” process, which normalizes the mean level density to unity, allowing for the extraction of system-specific features while minimizing errors through polynomial fitting or moving averages. The paper highlights various statistical measures, such as the nearest-neighbor spacing distribution (NNSD) $P(s)$ and the distribution of the ratio of consecutive level spacings $P(r)$, which are essential for probing short- and long-range correlations in energy levels.

The authors note that while the NNSD and other statistics require unfolding and may introduce errors, the distribution of the ratio of consecutive level spacings does not, making it a more reliable tool for analyzing short-range correlations. The introduction also references significant conjectures in quantum chaos, namely the Bohigas-Giannoni-Schmit (BGS) conjecture and the Berry-Tabor conjecture, which relate the statistical properties of energy levels to the classical behavior of the systems. The paper aims to investigate the statistical properties of energy levels in trapped interacting bosons under varying interaction strengths and rotational effects, an area that has not been thoroughly explored. The study promises to contribute valuable insights into the interplay between interaction energy and trap energy, with implications for understanding quantum many-body systems and their chaotic behavior.

Results

In this section, the authors investigate the spectral properties of a system of $N = 12, 16, 20$ bosonic atoms of $^{87}\text{Rb}$ confined in a quasi-2D harmonic trap, characterized by an anisotropy parameter $\lambda_z = 4$ and an axial frequency $\omega_z = 2\pi \times 220 \text{ Hz}$. The interaction between the atoms is modeled using a repulsive Gaussian potential with a range defined by $\sigma = 0.1 a_\perp$, where $a_\perp = 1.446 \, \mu\text{m}$ is the trap length. The study employs exact diagonalization to analyze the energy levels, focusing on two interaction regimes: moderate ($g_2 = 0.3669$) and strong ($g_2 = 3.669$), achieved by varying the s-wave scattering length $a_{sc}$.

The results reveal the dependence of energy level statistics on the number of bosons and interaction strength. The fitted Brody parameters $b$ for the non-rotating case indicate a transition from Poisson statistics (for $b = 0$) to Gaussian Orthogonal Ensemble (GOE) statistics (for $b = 1$). Specifically, for moderate interactions, the Brody parameters are $b = 0.33$ for $N = 12$, $b = 0.16$ for $N = 16$, and $b = 0.04$ for $N = 20$. In the strong interaction regime, the parameters are $b = 0.53$, $b = 0.85$, and $b = 0.51$, respectively. These findings suggest that as the number of bosons increases, the level statistics shift, reflecting the influence of two-body interactions and the system’s dimensionality on spectral correlations.

Discussion

In this section, the authors present a detailed analysis of a system of interacting bosons trapped in a quasi-2D harmonic potential and subjected to rotation. The many-body Hamiltonian is formulated in a co-rotating frame, incorporating kinetic energy, an external confining potential, and a two-body interaction modeled by a Gaussian potential. The effective interaction strength in 2D is derived, revealing that the system operates in a repulsive regime with a scattering length $a_{sc} > 0$. The eigenspectrum is obtained variationally, leading to a matrix eigenvalue problem that is solved to extract the energy levels.

The statistical analysis of the energy levels employs various tools to investigate both short-range and long-range correlations. The nearest-neighbor spacing distribution (NNSD) and the distribution of the ratio of consecutive level spacings are utilized to assess spectral correlations. The results indicate that in the moderate interaction regime, as the number of bosons increases, the NNSD approaches a Poisson distribution, suggesting a lack of level repulsion and indicating regular behavior. Conversely, in the strong interaction regime, the NNSD aligns more closely with the Gaussian Orthogonal Ensemble (GOE) distribution, indicating chaotic behavior characterized by strong level repulsion. The Dyson-Mehta $\Delta_3(L)$ statistic and level number variance $\Sigma^2(L)$ are also employed to analyze long-range correlations, revealing a transition from integrable to non-integrable behavior as the interaction strength and number of bosons vary. Overall, the findings highlight the intricate interplay between interaction strength, boson number, and spectral correlations in this quantum many-body system.