پژوهش یو دنگ، زاهر هانی و شیائو ما تحت عنوان «استنتاج معادله بولتزمن از کره‌های سخت در زمان‌های طولانی» مستقیماً به حل مسئله حد بولتزمن-گراد برای بازه‌های زمانی کلان در فضای اقلیدسی و چنبره سه‌بعدی می‌پردازد ((یو دنگ، زاهر هانی و شیائو ما، مقاله Long-time derivation of the Boltzmann equation from hard spheres، آرشیو arXiv:2408.07818، سال ۲۰۲۴-۲۰۲۵ میلادی)). در حالی که قضایای پیشین به رهبری اسکار لنفورد (۱۹۷۵) تنها برای کسری از میانگین زمان بین دو برخورد متوالی معتبر بودند، این مقاله با معرفی روش «برش خوشه‌ها» و دسته‌بندی توپولوژیک گراف‌های تاریخچه برخورد، همگرایی دینامیک میکروسکوپیک ذرات به [[معادله بولتزمن]] را در تمام طول عمر جواب‌های کلاسیک این معادله اثبات می‌کند.
<قضیه اصلی (قضیه ۱.۱ مقاله): فرض کنید تابع توزیع اولیه ${f_0(x, v)}$ هموار، نامنفی و دارای افت نمایی در سرعت باشد به طوری که معادله بولتزمن در بازه زمانی ${[0, T_{\text{max}}]}$ دارای جواب قوی و پیوسته ${f(t, x, v)}$ باشد. آنگاه در حد بولتزمن-گراد (${N \to \infty}$ و ${\varepsilon \to 0}$ با قید ${N \varepsilon^{d-1} = 1}$)، چگالی حاشیه‌ای تک‌ذره‌ای سیستم کره‌های سخت به صورت یکنواخت در کل بازه ${[0, T_{\text{max}}]}$ به جواب معادله بولتزمن همگرا می‌شود.>
# دستگاه دینامیکی میکروسکوپیک و رژیم حدی بولتزمن-گراد #
یک سیستم متشکل از ${N}$ کره سخت همسان با قطر ${\varepsilon > 0}$ و جرم واحد را در دامنه فضایی ${\Lambda}$ (که می‌تواند ${\mathbb{R}^d}$ یا چنبره متناوب ${\mathbb{T}^d}$ با ${d \ge 2}$ باشد) در نظر می‌گیریم. وضعیت میکروسکوپیک سیستم در هر لحظه ${t}$ با بردار فاز زیر مشخص می‌شود:
$$
\boldsymbol{z}_N = (z_1, z_2, \dots, z_N) \in \mathcal{D}_N, \quad z_j = (x_j, v_j) \in \Lambda \times \mathbb{R}^d
$$
فضای فاز مجاز ${\mathcal{D}_N}$ شامل پیکربندی‌هایی است که در آن‌ها مراکز کره‌ها به اندازه حداقل ${\varepsilon}$ از یکدیگر فاصله دارند:
$$
\mathcal{D}_N = \left\{ \boldsymbol{z}_N \in (\Lambda \times \mathbb{R}^d)^N : \forall i \ne j, \, |x_i - x_j| \ge \varepsilon \right\}
$$
در درون فضای ${\mathcal{D}_N}$، ذرات به صورت آزاد و با سرعت ثابت حرکت می‌کنند (${\dot{x}_j = v_j, \, \dot{v}_j = 0}$). به محض رسیدن فاصله دو ذره ${i}$ و ${j}$ به حد تماس ${|x_i - x_j| = \varepsilon}$ با بردار یکه برخورد ${\omega = (x_i - x_j)/\varepsilon \in \mathbb{S}^{d-1}}$ (به طوری که ${(v_i - v_j) \cdot \omega < 0}$)، یک برخورد کشسان آنی رخ می‌دهد و سرعت‌های پس از برخورد (${v'_i, v'_j}$) طبق قوانین پایستگی تکانه و انرژی جنبشی تعیین می‌شوند:
$$
v'_i = v_i - \left[ (v_i - v_j) \cdot \omega \right] \omega, \quad v'_j = v_j + \left[ (v_i - v_j) \cdot \omega \right] \omega
$$
**حد بولتزمن-گراد «Boltzmann-Grad Limit»:** برای آنکه میانگین مسیر آزاد ذرات ${\ell_{\text{mfp}}}$ از مرتبه ${1}$ ثابت بماند (یعنی هر ذره در واحد زمان به طور میانگین تعداد متناهی برخورد را تجربه کند)، قطر ذرات ${\varepsilon \to 0}$ میل داده می‌شود در حالی که میانگین تعداد ذرات ${N}$ طبق رابطه مقیاس‌بندی زیر به بی‌نهایت می‌رود:
$$
N \varepsilon^{d-1} = 1
$$
# توزیع آماری گروه بزرگ و معادلات سلسله‌مراتبی #
به جای فیکس کردن تعداد دقیق ذرات ${N}$، مقاله از «آنسامبل کانونی بزرگ» «Grand Canonical Ensemble» استفاده می‌کند که در آن تعداد کل ذرات ${\mathcal{N}}$ خود یک متغیر تصادفی پواسون با میانگین متناسب با ${\varepsilon^{-(d-1)}}$ است. در لحظه ${t = 0}$، چگالی احتمال یافتن دقیقاً ${n}$ ذره در موقعیت فاز ${\boldsymbol{z}_n}$ به صورت حاصل‌ضرب توزیع‌های مستقل اولیه (فرض آشفتگی مولکولی اولیه) به همراه شرط عدم هم‌پوشانی کره‌ها تعریف می‌شود:
$$
W_0^{(n)}(\boldsymbol{z}_n) = \frac{1}{\mathcal{Z}} \frac{\varepsilon^{-n(d-1)}}{n!} \mathbb{I}_{\mathcal{D}_n}(\boldsymbol{z}_n) \prod_{j=1}^n f_0(x_j, v_j)
$$
که در آن ${\mathcal{Z}}$ تابع افراز بهنجارساز و ${\mathbb{I}_{\mathcal{D}_n}}$ تابع مشخصه عدم تداخل کره‌ها است. با تحول زمانی سیستم تحت جریان همیلتونی کره‌های سخت ${\Phi_t^{(n)}}$، توابع هم‌بستگی حاشیه‌ای ${s}$-ذره‌ای «Correlation Functions» در زمان ${t}$ عبارت‌اند از:
$$
f_s^\varepsilon(t, \boldsymbol{z}_s) = \varepsilon^{s(d-1)} \sum_{k=0}^\infty \frac{1}{k!} \int_{(\Lambda \times \mathbb{R}^d)^k} W_t^{(s+k)}(\boldsymbol{z}_s, \boldsymbol{z}'_k) \, d\boldsymbol{z}'_k
$$
هدف نهایی اثبات این است که به ازای هر زمان کلان ${t \in [0, T_{\text{max}}]}$، تابع هم‌بستگی یک‌ذره‌ای ${f_1^\varepsilon(t, z_1)}$ به جواب معادله غیرخطی بولتزمن همگرا شود:
$$
\left( \partial_t + v \cdot \nabla_x \right) f(t, x, v) = \int_{\mathbb{R}^d} \int_{\mathbb{S}^{d-1}} \left[ (v - v_*) \cdot \omega \right]_+ \left( f(t, x, v') f(t, x, v'_*) - f(t, x, v) f(t, x, v_*) \right) d\omega \, dv_*
$$
و برای ${s \ge 2}$، خاصیت فاکتورگیری یا **«انتشار آشفتگی مولکولی»** «Propagation of Chaos» خارج از مجموعه‌ای با اندازه صفر برقرار بماند:
$$
\lim_{\varepsilon \to 0} f_s^\varepsilon(t, z_1, \dots, z_s) = \prod_{j=1}^s f(t, z_j)
$$
# علت شکست روش لنفورد و سری‌های دایسون در زمان‌های طولانی #
در روش کلاسیک اسکار لنفورد، تابع هم‌بستگی ${f_1^\varepsilon(t)}$ با استفاده از انتگرال‌گیری مکرر از سلسله‌مراتب معادلات BBGKY به صورت یک سری توانی بسط داده می‌شود. جملات این سری متناظر با «درخت‌های برخورد» هستند که تاریخچه برخورد ذره شماره ۱ با ذرات دیگر را تا زمان ${t = 0}$ به عقب ردیابی می‌کنند.
مسئله بنیادین در زمان‌های طولانی از رشد سریع فاکتوریل تعداد درخت‌های ممکن ناشی می‌شود:
* در یک درخت با ${k}$ برخورد، تعداد توپولوژی‌های ممکن برای درخت از مرتبه ${k!}$ است.
* از سوی دیگر، انتگرال زمانی روی ${k}$ لحظه برخورد متوالی ${0 < t_k < \dots < t_1 < t}$ حجمی برابر با ${\frac{t^k}{k!}}$ تولید می‌کند.
* در نتیجه، جمله ${k}$-ام سری با ضریب حدود ${(C t)^k}$ رشد می‌کند و سری تنها برای زمان‌های بسیار کوتاه ${t < t_{\text{Lanford}} = C^{-1}}$ همگراست!
علت فیزیکی این واگرایی آن است که در بسط مستقیم BBGKY، جملات مثبت (ورود به حالت از طریق برخورد) و جملات منفی (خروج از حالت به دلیل برخورد یا هم‌پوشانی ممنوعه کره‌ها) از هم جدا می‌شوند و اثر حذف‌کنندگی متقابل آن‌ها که ضامن پایستگی احتمال کل (${\|f_1^\varepsilon\|_{L^1} = 1}$) است، از دست می‌رود.
# نوآوری اول: خوشه‌های برخورد و الگوریتم برش لایه‌بندی‌شده #
دنگ، هانی و ما به جای بسط سلسله‌مراتب BBGKY، مستقیماً دینامیک سیستم را در فضای پیکربندی کامل تحلیل کرده و زمان کل ${[0, T_{\text{max}}]}$ را به تعداد زیادی گام زمانی کوچک به طول ${\tau \ll 1}$ تقسیم می‌کنند:
$$
0 = \tau_0 < \tau_1 < \dots < \tau_L = T_{\text{max}}, \quad \tau_{\ell} - \tau_{\ell-1} = \tau
$$
در یک بازه زمانی، دو ذره را در صورتی که با یکدیگر برخورد کنند (یا در فاصله کمتر از ${\varepsilon}$ از هم رد شوند به عنوان برخورد مجازی) با یک یال به هم وصل می‌کنیم. مؤلفه‌های هم‌بند این گراف در فضا-زمان را **«خوشه‌های دینامیکی»** «Dynamical Clusters» می‌نامیم.
اگر اندازه یک خوشه (تعداد ذرات درگیر در زنجیره برخوردهای متوالی با ذره ریشه) از یک حد آستانه بزرگ مانند ${K_{\text{max}} = |\log \varepsilon|^C}$ فراتر رود، محاسبه مستقیم آن ناممکن می‌شود. نویسندگان برای حل این مشکل، **عملیات برش خوشه** «Cluster Cutting Operation» را معرفی می‌کنند:
* در هر لایه زمانی ${[\tau_{\ell-1}, \tau_\ell]}$، اگر اندازه درخت برخوردهای متصل به یک ذره به عدد ثابت و بزرگی مانند ${\Lambda_0}$ برسد، آخرین برخورد «قطع» می‌شود؛ یعنی به جای اعمال برخورد، فرض می‌شود دو ذره بدون دیدن یکدیگر از هم عبور کرده‌اند!
* خطای حاصل از این برش در یک گام زمانی ${\tau}$ از مرتبه ${\tau^{\Lambda_0}}$ است که با انتخاب ${\Lambda_0}$ به اندازه کافی بزرگ، در طول کل زمان ${T_{\text{max}}}$ کاملاً ناچیز می‌شود.
# نوآوری دوم: فرمول شمول-عدم شمول و گراف‌های با حلقه #
با اعمال قانون برش، در هر گام زمانی تنها با خوشه‌هایی سر و کار داریم که حداکثر شامل تعداد پلی‌لگاریتمی ذره هستند. اما برای چسباندن گام‌های زمانی ${\tau_\ell}$ به یکدیگر، باید شرط «عدم هم‌پوشانی و عدم برخورد با ذرات خارج از خوشه» را به زبان ریاضی بازنویسی کنیم.
فرض کنید ${\mathbb{I}_{\text{no-ext}}(H)}$ تابع مشخصه‌ای باشد که تضمین می‌کند ذرات داخل خوشه ${H}$ با هیچ ذره خارجی دیگری برخورد نکرده‌اند. به جای نگه‌داشتن این شرط غیرخطی، از **بسط ترکیبیاتی شمول و عدم شمول (بسط موبیوس)** استفاده می‌شود:
$$
\prod_{e \in \mathcal{E}_{\text{ext}}} (1 - \chi_e) = \sum_{E \subset \mathcal{E}_{\text{ext}}} (-1)^{\# E} \prod_{e \in E} \chi_e
$$
که در آن ${\chi_e}$ نشان‌دهنده وقوع یک «هم‌پوشانی یا برخورد مجازی (سایه‌ای)» روی یال ${e}$ است. این بسط به ما اجازه می‌دهد احتمال دقیق یک خوشه را به صورت مجموع علامت‌داری روی **«گراف‌های تعمیم‌یافته برخورد»** ${\mathcal{G}}$ بنویسیم که شامل دو نوع یال هستند:
1. **یال‌های برخورد واقعی (فیزیکی):** که مسیر و سرعت ذرات را تغییر می‌دهند.
2. **یال‌های هم‌پوشانی (سایه‌ای / منفی):** که وزن ${-1}$ دارند و ذرات در آن‌ها بدون تغییر سرعت از درون هم عبور می‌کنند.
<تفکیک گراف‌های درختی و گراف‌های دارای دور (حلقه): اگر گراف تاریخچه برخورد میان ذرات کاملاً یک «درخت بدون حلقه» باشد، به این معناست که هر دو ذره حداکثر یک بار با هم روبرو شده‌اند و هیچ هم‌بستگی پیشینی میان آن‌ها وجود ندارد. این گراف‌های درختی دقیقاً جملات سری معادله بولتزمن را تولید می‌کنند! در مقابل، گراف‌هایی که دارای حلقه «Cycle» هستند، نماینده پدیده مخرب **برخورد مجدد «Recollision»** می‌باشند که در آن دو ذره که در گذشته مستقیم یا غیرمستقیم با هم برخورد کرده‌اند، دوباره در آینده به هم می‌رسند و استقلال آماری را نقض می‌کنند.>
# مهار برخوردهای مجدد و اثبات ناچیز بودن گراف‌های حلقه‌دار #
بخش عمده مقاله (بیش از ۱۰۰ صفحه از اثبات) به اثبات یک نامساوی ترکیبیاتی-هندسی اختصاص دارد که نشان می‌دهد مجموع سهم تمام گراف‌های دارای حلقه (یعنی تمام تاریخچه‌هایی که در آن‌ها برخورد مجدد رخ می‌دهد) در حد ${\varepsilon \to 0}$ با توانی از ${\varepsilon}$ به صفر میل می‌کند:
$$
\sum_{\mathcal{G} \text{ has cycles}} \left| \mathcal{I}(\mathcal{G}) \right| \lesssim \varepsilon^{\vartheta}, \quad \vartheta > 0
$$
برای اثبات این کران، نویسندگان تکنیک‌های زیر را ابداع کرده‌اند:
## ۱. انتگرال‌گیری روی پارامترهای برخورد در حضور یک حلقه ##
هر حلقه بنیادی در گراف تاریخچه به معنای وجود یک شرط هندسی اضافی است: یک ذره پس از طی زنجیره‌ای از پراکندگی‌ها، باید دقیقاً در لحظه مناسب از استوانه‌ای به شعاع بسیار کوچک ${\varepsilon}$ در اطراف ذره‌ای دیگر عبور کند. در فضای سه‌بعدی (${d = 3}$)، احتمال اینکه دو خط مسیر تصادفی از فاصله ${\varepsilon}$ یکدیگر عبور کنند از مرتبه ${\varepsilon^{d-1}}$ است، که در نگاه اول ضریب بزرگِ تعداد ذرات ${\varepsilon^{-(d-1)}}$ را خنثی نمی‌کند؛ بلکه به دلیل محدود بودن زاویه فضایی پراکندگی، هر حلقه مستقل یک ضریب کوچک اضافی از مرتبه زیر وارد انتگرال می‌کند:
$$
\text{Gain per independent cycle} \lesssim \varepsilon^{d-1} \cdot \left( \frac{\varepsilon}{|t_1 - t_2|} \right) \implies \varepsilon^{1 - \delta} \text{ after time integration}
$$
## ۲. الگوریتم تجمیع هسته‌ها برای حلقه‌های متداخل ##
در زمان‌های طولانی، ممکن است یک گراف دارای ده‌ها حلقه تو در تو و متقاطع باشد (گراف‌های با عدد بتی اول ${\beta_1(\mathcal{G}) = \# E - \# V + 1}$ بزرگ). اگر حلقه‌ها به هم گره خورده باشند، شرایط هندسی آن‌ها ممکن است به هم وابسته شوند (تکینگی ژاکوبین تبدیل مختصات).
دنگ، هانی و ما یک **الگوریتم کاهش گام‌به‌گام گراف** طراحی کردند که در آن رئوس و یال‌های گراف بر اساس ترتیب زمانی برخوردها یکی‌یکی حذف و انتگرال‌گیری می‌شوند. در هر مرحله، ذرات به «مولکول‌های مجازی» تبدیل می‌شوند و اثبات می‌شود که به ازای هر حلقه توپولوژیک، دست‌کم یک درجه آزادی مکانی یا سرعتی در فضای فاز به عرض ${\varepsilon}$ مقید می‌شود، بدون آنکه تکینگی‌های انتگرال‌پذیر با هم تداخل مخرب داشته باشند.
# معادله کومولانت، جفت‌سازی با معادله بولتزمن و بستن حلقه اثبات #
پس از حذف گراف‌های دارای حلقه و برش خوشه‌های غول‌پیکر، تابع توزیع سیستم در هر گام زمانی ${\tau_{\ell}}$ به صورت زیر تجزیه می‌شود:
$$
f_1^\varepsilon(\tau_\ell) = \mathcal{T}_{\text{tree}}\left[ f_1^\varepsilon(\tau_{\ell-1}) \right] + \mathcal{R}_{\text{cycle}}(\tau_\ell) + \mathcal{R}_{\text{cut}}(\tau_\ell)
$$
که در آن عملگر درختی ${\mathcal{T}_{\text{tree}}}$ دقیقاً عملگر تحول معادله بولتزمن در بازه زمانی ${[\tau_{\ell-1}, \tau_\ell]}$ است و جملات باقی‌مانده ${\mathcal{R}_{\text{cycle}}}$ و ${\mathcal{R}_{\text{cut}}}$ از مرتبه ${\mathcal{O}(\varepsilon^\vartheta)}$ هستند.
در گام نهایی، با خطی‌سازی معادله بولتزمن حول جواب هموار ${f(t, x, v)}$ و استفاده از پایداری عملگر خطی‌شده بولتزمن در فضاهای وزن‌دار سرعت، نشان داده می‌شود که خطای انباشته‌شده در طول ${L = T_{\text{max}}/\tau}$ گام زمانی، از مرتبه ${\varepsilon^{\vartheta'}}$ باقی می‌ماند:
$$
\sup_{t \in [0, T_{\text{max}}]} \left\| f_1^\varepsilon(t) - f(t) \right\|_{L^\infty_{x, v}} \le C(T_{\text{max}}) \, \varepsilon^{\vartheta'} \xrightarrow{\varepsilon \to 0} 0
$$
این نتیجه، بخش جنبشی **مسئله ششم هیلبرت** را پس از ۱۲۵ سال به طور کامل حل می‌کند: قوانین میکروسکوپیک نیوتن که در زمان کاملاً برگشت‌پذیر هستند (${t \mapsto -t}$)، در حد آماری بولتزمن-گراد و بدون هیچ محدودیت زمانی مصنوعی، به معادله برگشت‌ناپذیر بولتزمن (دارای قضیه افزایش آنتروپی ${H}$) منجر می‌شوند؛ و منشأ این برگشت‌ناپذیری در انتخاب جهت زمانیِ فرض استقلال آماری اولیه پیش از برخوردها نهفته است.