نظرية النحلة – الاشتقاق النظري – 2025 مايو 17 مع كلود

من ψ = ψ = exp(-αr) إلى F = -G/R²: اشتقاق نظرية النحلة الكاملة

لماذا تكون الدالة الموجية ψ(r) = N exp(-αr) صحيحة – ولكن يجب التعامل مع طريقة إسقاطها بالقرب من جسيم ثانٍ بعناية. وينتج عن الإسقاط المصحح قانون قوة يوكاوا-نيوتن الذي يختزل إلى قانون المربع العكسي لنيوتن داخل طول التماسك.

موقع BeeTheory.com – ملحق وتصحيح BeeTheory v2 (Dutertre 2023)

ψ(r) = N ^e-r/r

شكل الدالة الموجية الجسيمية الصحيحة

ψA _|B =_CA ^e-αr/R

خطوة الإسقاط المحلي الرئيسية

V(R) ∝ ^-e-αR/R

إمكانات نظرية النحل

و → – ك/ر²

قانون نيوتن داخل طول التماسك

0. الجواب – مذكور أولاً

الدالة الموجية لنظرية النحل ψ(r) = N exp(-αr) صحيحة ولا تحتاج إلى تغيير. التعديل الذي ينتج عنه F ∝ 1/R² ليس في شكل ψ، ولكن في كيفية تقييم _ψA بالقرب من الجسيم B.

عندما يتم إسقاط موجة A حول موقع B، باستخدام مفكوك تايلور لـ exp(-α|AP|) لـ r = |BP|، تكون النتيجة

\(\psi_A\big|_{\mathrm{near}\ B}=C_A(R)e^{-\alpha r\cos\theta}\)

يُعطي متوسط أحادي القطب الكروي:

\(\\ppsi_A\big|{\mathrm{near}\B}=C_A(R)\frac{\sinh(\alpha r)}{\alpha r}\)

عند الرتبة الرئيسية في r، تكون sinh(αr)/(αr)≈ 1، ومن ثم تبدو الموجة من A ثابتة تقريبًا بالقرب من B. يدخل الاعتماد على R من خلال السعة:

\(C_A(R)=Ne^^{- \ألفا R}\)

باستخدام الإسقاط المحلي لنظرية النحلة، يُكتب معدل التضاؤل المحلي الفعال على الصورة α/R. ينتج عن تطبيق لابلاسيان على هذه الموجة المسقطة المحلية حد مهيمن يتناسب مع 1/(Rr). يعمل هذا كجهد شبيه بجهد كولوم 1/ص بالقرب من B. بعد التكامل على الدالة الموجية B، يصبح جهد التفاعل

\(V(R)=- \frac{K e^{- \ألفا R}}{R}\)

القوة إذن:

\(F(R)=-\frac{dV}{dR}=-\frac{K(1+\alpha R)e^{-\alpha R}}{R^2}\)

داخل طول التماسك، حيث R ≪ ≪ ℓ = 1/α، يكون لدينا ^e-αR ≈ 1 و1 + αR ≈ 1، إذن:

\(\boxed{F(R)\xrightarrow{\alpha R\ll1}-\frac{K}{R^2}}\)

هذا هو قانون المربع العكسي لنيوتن. طول التماسك ℓ هو المدى الذي تتصرف فيه الجاذبية كقوة نيوتونية.

1. دالة موجة الجسيمات – الشكل الثلاثي الأبعاد الدقيق

تُصوِّر نظرية النحلة كل جسيم ذي كتلة على أنه دالة موجية متماثلة كرويًّا تتحلل أسيًّا من مركزها. بالنسبة لجسيم طول تماسكه ℓ = 1/α:

\(\psi(\mathbf{r})=Ne^{-\alpha|\mathbf{r}|}=\frac{\alpha^{3/2}}{\sqrt{\pi}}e^{-r/\ell}\)

شرط التطبيع هو:

\(\\int \\psi |^2d^3r=4\pi N^2\int_0^\^\infty e^{2\alpha r}r^2dr=1\) \(N=\frac{\alpha^{3/2}}{\sqrt{\pi}}\)

هذا الشكل له مركز مدمج على شكل جرس، ويصل إلى الحد الأقصى عند r = 0، ويظل محدودًا في كل مكان، ويتضاءل إلى الصفر مع اقتراب r من اللانهاية. وهو يمثل جسيمًا موضعيًا ذا طابع موجي يمتد إلى ما بعد مركزه.

في ميكانيكا الكم، بالنسبة لذرة الهيدروجين، هذه هي بالضبط الدالة الموجية للحالة الأرضية 1s مع α = 1/a0، حيث a0 هو نصف قطر بوهر. وهذا يعطي طاقة الحالة الأرضية المعروفة _E1s = -13.6 eV.

اللابلاسيان الدقيق في الإحداثيات الكروية ثلاثية الأبعاد

\(\nabla^2\left(e^{-\alpha r}\right)=\frac{d^2}{dr^2}e^{-\alpha r}+\frac{2}{r}\frac{d}{dr}e^{-\alpha r}\) \(=\alpha^2e^{-\alpha r}-\frac{2\alpha}{r}e^{-\alpha r}=e^{-\alpha r}\left(\alpha^2-\frac{2\alpha}{r}\right)\)

ما الذي تفعله الورقة الأصلية ولماذا تفقد الاعتماد على R

يكتب البحث الأصلي _ψA(r بالقرب من B) = C exp(_-αr/RAB) ويحسب لابلاسيان على أنه _-3α/RAB تقريبًا، وهو ثابت. يعطي هذا التقريب أحادي القطب طاقة ثابتة بدلًا من إمكانية، لأنه يفقد اعتماد القوة على R.

يوضح الاشتقاق المصحح أن نتيجة -3α/R يمكن تفسيرها على أنها معامل محلي، ولكن لا يجب إيجاد قيمة لابلاسيان عند r = 0 فقط. يجب أن يتم تكاملها على حجم الدالة الموجية B. وهذا ما يعيد قانون القوة الصحيح.

2. إسقاط _ψA بالقرب من B – الخطوة الرئيسية

ضع الجسيم A عند نقطة الأصل والجسيم B عند الموضع R على طول المحور z. افترض أن نقطة المجال P عند الموضع r عند الموضع r مقيسة من B، بزاوية قطبية θ من المحور AB، مع r ≪ R.

2.1 المسافة الدقيقة من A إلى P

\(|AP|^2=R^2+r^2+2Rr\cos\theta\) \(|AP|=R\sqrt{1+\frac{2r\cos\theta}{R}+\frac{r^2}{R^2}}\) \(|AP]|AP \r\r\r\cos\theta\qquad\text{لـ}r\ll R\)

لذلك:

\(\psi_A(P)=Ne^^{-\alpha \alpha |AP|}=Ne^{-\alpha R}e^{-\alpha r\cos\theta}=C_A(R)e^{-\alpha r\cos\theta}\) [اللاتكس]C_A(R)=Ne^^{- \ألفا R}[/latex]

2.2 متوسط القطب الأحادي الكروي

بحساب المتوسط على جميع الاتجاهات θ، وهو أمر مناسب عندما تكون الدالة الموجية B متماثلة كرويًّا، نحصل على

\(\langle\psi_A\rangle_\theta=C_A(R)\frac{1}{2}\int_0^\pi e^{-\alpha r\cos\theta}\sin\theta\,d\theta\) \(=C_A(R)\frac{\sinh(\alpha r)}{\alpha r}\)

داخل طول التماسك، عندما يكون r ≪ ℓ = 1/α:

\(\frac{\sinh(\alpha r)}{\alpha r}\approx1+\frac{(\alpha r)^2}{6}+\cdots\approx1\)

تظهر الموجة A ثابتة محليًّا بالقرب من B. ويهيمن على التفاعل السعة_CA(R).

تستخدم ورقة “نظرية النحل” التقريب المحلي:

\(\\ppsi_A\big|{\mathrm{near}\B}(r)=C_A(R)e^{-(\ألفا/ر)r}\)

وهذا يعامل التضاؤل المحلي الفعال على أنه _βeff = α/R. هذه هي الخطوة التي تُدخِل 1/R في المشغِّل المحلي وتولِّد في النهاية القوة العكسية المربعة.

\(\beta_{\mathrm{eff}}=\frac{\alpha}{R}\)

مع نمو R، تبدو الموجة من A مسطحة بشكل متزايد في المنطقة المجاورة لـ B. هذه هي آلية نظرية النحل للقوة بعيدة المدى.

3. لابلاسيان الموجة المسقطة – من أين يأتي 1/R²

3.1 اللابلاسيان الدقيق لـ ^e-βr مع β = α/R

\(\nabla^2\left(e^{-\beta r}\right)=e^{-\beta r}\left(\beta^2-\frac{2\beta}{r}\right)\) \(=e^{-\alpha r/R}\left(\frac{\alpha^2}{R^2}-\frac{2\alpha}{Rr}\right)\)

هذا له مصطلحان مختلفان من الناحية الهيكلية:

المدة	التعبير	السلوك	الدور البدني
ثابت الحركة	α²هـ ^{²ص/ص/ر/ر²}	منتهية عندما تكون r → 0	يساهم في تحول مستمر في الطاقة.
مولد كولوم	^{-2αe-αe-αr/R/}(Rr)	يتباعد بمقدار 1/ص	توليد جهد محلي شبيه بجهد كولوم بمعامل يتناسب مع 1/R.

تطبيق عامل الحركة على الموجة المحلية A بالقرب من B:

\(\hat{T}\left[C_A(R)e^{-\alpha r/R}\right]=C_A(R)e^{-\alpha r/R}\left[\frac{\hbar^2\alpha}{mR}\frac{1}{r}-\frac{\hbar^2\alpha^2}{2mR^2}\right]\)

3.2 طاقة التفاعل – التكامل على حجم B

طاقة تفاعل نظرية النحلة هي عنصر مصفوفة هذا العامل الحركي مع الدالة الموجية B:

\(V_{\mathrm{BT}}(R)=C_A(R)\left[\frac{\hbar^2\alpha}{mR}I_1(R)-\frac{\hbar^2\alpha^2}{2mR^2}I_2(R)\right]\)

حيث:

\(I_1(R)= \left\langle\psi_B\middle\\\frac{e^^^- \ألفا ص/ر}{r}{r}\middle\\psi_B\right\rangle\) \(I_2(R)=\left\langle\psi\B\B\middle\e^^{-\alpha r\R}\middle\\psi_B\right\rangle\)

بالوحدات الذرية، هذه التكاملات هي:

\(I_1(R)=\frac{4}{\pi}\int_0^\infty e^{-(2+1/R)r}r\,dr=\frac{4}{\pi(2+1/R)^2}\) \(I_2(R)=\frac{4}{\pi}\int_0^\infty e^{-(2+1/R)r}r^2\,dr=\frac{8}{\pi(2+1/R)^3}\)

عند R كبيرة، تقترب هذه الثوابت من الثوابت:

\(I_1(R)\xrightarrow{R\gg\ell}\frac{1}{\pi},\qquad I_2(R)\xrightarrow{R\gg\ell}\frac{1}{\pi}\)

تصبح الإمكانية:

\(V_{\mathrm{BT}}(R)\xrightarrow{R\gg\ell}-C_A(R)\frac{\hbar^2\alpha}{\pi mR}\left(1-\frac{\alpha}{2R}\right)\) \(V_{\mathrm{BT}}(R)\approx-\frac{\hbar^2\alpha}{\pi m}\frac{Ne^{-\alpha R}}{R}\) \(\boxed{V_{{\mathrm{BT}}(R)=-\frac{\Ke^{{-\ألفا R}}{R}، \qquad K=\frac{\hbar^2\ألفا N}{\pi m}}\)

4. القوة – بروز قانون نيوتن – قانون نيوتن

بدءاً من إمكانات نظرية النحلة:

\(V(R)=-K\frac{e^{-\alpha R}}{R}\)

القوة

\(F(R)=-\frac{dV}{dR}=-\frac{d}{dR}\left[-K\frac{e^{-\alpha R}}{R}\right]\) \(\boxed{F(R)=-\frac{K(1+\alpha R)}{R^2}e^{-\alpha R}}\)

تحتوي هذه الصيغة الواحدة على ثلاثة أنظمة.

I. نظام الجاذبية: R ≪ ≪ ℓ ℓ

\(e^^{- \ألفا R} \approx1,\qquad 1+ \ألفا R\approx1\) \(F(R)\approx-\frac{K}{R^2}\)

هذا هو قانون المربع العكسي لنيوتن. تظهر الجاذبية على شكل 1/R² عند مقاييس أصغر من طول التماسك.

ثانياً: النظام الانتقالي النظام الانتقالي: R∼ ℓ ℓ

\(F(R)=- \\frac{K(1+\ألفا R)}{R^2}e^{-\ألفا R}\)

يبدأ العامل الأسي في كبح القوة. هذا هو النظام الذي تصبح فيه الانحرافات عن القياس النيوتوني قابلة للقياس.

ثالثاً. نظام يوكاوا R ≫ ≫ ℓ ℓ

\(F(R)\approx- \frac{K\alpha}{R}e^^{-\ألفا R}\)

تصبح القوة مكبوتة أضعافًا مضاعفة. هذا هو نظام يوكاوا قصير المدى.

4.1 التحقق العددي: F(R) – R²

بالنسبة لقانون المربع العكسي النيوتوني المثالي، يجب أن يكون حاصل الضرب F(R) – R² ثابتًا. عامل تصحيح نظرية النحل هو

\(\frac{F(R)R^2}{K}=(1+\alpha R)e^{-\alpha R}=\left(1+\frac{R}{\ell}\right)e^{-R/\ell}\)

عندما يكون R/ℓ صغيرًا، يظل هذا العامل قريبًا من 1.

ص/ب	^هـ-ر/ب	(1 + ص/ص/م) ^{هـ-ص/ر/م}	الخطأ مقابل الخطأ النقي 1/ص²	النظام
0.01	0.9900	0.9999	<0.01%	نيوتن
0.05	0.9512	0.9988	0.12%	نيوتن
0.10	0.9048	0.9953	0.47%	نيوتن
0.30	0.7408	0.9631	3.7%	بدء المرحلة الانتقالية
0.50	0.6065	0.9098	9.0%	المرحلة الانتقالية
1.00	0.3679	0.7358	26.4%	النظام المختلط
2.00	0.1353	0.4060	59.4%	يوكاوا المهيمن
5.00	0.0067	0.0403	96%	الاضمحلال الأسي

الرسم البياني المقترح: ارسم F(R) – R² / K مقابل R لأطوال تماسك مختلفة ℓ. بالنسبة ل ℓ الكبير جدًا، يظل المنحنى مسطحًا تقريبًا عند 1، مما يُظهر السلوك النيوتوني. بالنسبة إلى ℓ الأصغر، ينخفض المنحنى أسيًا.

5. المعادلات الكاملة – جميع الخطوات

الخطوة 1 – دالة موجة الجسيمات

\(\psi(r)=\frac{\alpha^{3/2}}{\sqrt{\pi}}e^{-\alpha r},\qquad \alpha=\frac{1}{\ell},\qquad N=\frac{\alpha^{3/2}}{\sqrt{\pi}}\) \(\psi(0)=N,\qqquad \psi(\infty)=0,\qquad \int \\psi \^2d^3r=1\)

الخطوة 2 – إسقاط موجة A بالقرب من B

\(\psi_A\big|_{\mathrm{near}\ B}(r)=\underbrace{Ne^{-\alpha R}}_{C_A(R)}\underbrace{e^{-(\alpha/R)r}}_{\text{local shape}}\) \(\beta_{\mathrm{eff}}=\frac{\alpha}{R}\)

الخطوة 3 – اللبلاسيان الدقيق للموجة المحلية

\(\nabla^2\left(e^{-\alpha r/R}\right)=e^{-\alpha r/R}\left(\frac{\alpha^2}{R^2}-\frac{2\alpha}{Rr}\right)\)

الحد -2α/(Rr) هو أصل الجهد المحلي الشبيه بجهد كولوم وبالتالي القوة العكسية المربعة.

الخطوة 4 – عناصر المصفوفة على دالة الموجة B

\(\left\langle\psi_B\middle|\frac{e^{-\alpha r/R}}{r}\middle|\psi_B\right\rangle=\frac{4}{\pi(2+\alpha/R)^2}\) [لاتكس] \ latex]\ يسار \langle\psi_B\middle \middle \e^{-\alpha r\rangle=\frac{8}{\pi(2+\alpha/R)^3}[/latex] \(\xrightarrow{R\gg\ell}\frac{1}{\pi},\qquad \frac{1}{\pi}\)

الخطوة 5 – إمكانات التفاعل والقوة

\(V_{\mathrm{BT}}(R)=-\frac{Ke^{-\alpha R}}{R}\) \(K=\frac{\hbar^2\alpha N}{\pi m}=\frac{\hbar^2}{\pi m\ell}\frac{1}{\sqrt{\pi}\ell^{3/2}}\) \(\boxed{F(R)=-\frac{K(1+\alpha R)e^{-\alpha R}}{R^2}=-\frac{K}{R^2}\left(1+\frac{R}{\ell}\right)e^{-R/\ell}}\) [لاتكس]\مربع{R\R\ll\ll\ll\l\l\l\l\l\l\longright\longrightarrow\lquad F(R)=-\frac{K}{R\l^2}}[/latex]

مع:

\(K=\frac{\hbar^2\alpha^{5/2}}{\pi^{3/2}m},\qquad \alpha=\frac{1}{\ell}=\frac{mv_{\mathrm{wave}}}{\hbar}\)

5.1 تحديد ثابت نيوتن G

بالنسبة إلى كتلتين _m1 و _m2، تتطلب النهاية النيوتونية:

\(F(R)=-\frac{Gm_1m_2}{R^2}\)

بالمقارنة مع حد نظرية النحلة F = -K/R²، نحصل على:

\(G=\frac{K}{m_1m_2}=\frac{\hbar^2\alpha^{5/2}}{\pi^{3/2}m\,m_1m_2}\)

بالنسبة إلى _m1 = _m2 = م:

\(G=\frac{\hbar^2}{\pi^{3/2}m^2\ell^{5/2}}\)

إيجاد قيمة ℓ:

\(\ell=\left(\frac{\hbar^2}{\pi^{3/2}Gm^2}\right)^{2/5}\)

بالنسبة لكتلة البروتون_mp = 1.67 × ^10-27 كجم:

\(\ell_p=\left(\frac{(1.055\times10^{-34})^2}{\pi^{3/2}\times6.674\times10^{-11}\times(1.67\times10^{-27})^2}\right)^{2/5}\approx1.2\times10^{13}\,\mathrm{m}\approx80\,\mathrm{AU}\)

هذا هو طول ترابط الجاذبية للبروتون في هذا المقياس المبسط. بالنسبة للأجسام العيانية الكبيرة، فإن طول التماسك الفعال يتدرج مع المجال الموجي الكلي لجميع الجسيمات المكوِّنة.

6. ملخص: الورقة الأصلية مقابل الاشتقاق المصحح

ورقة نظرية النحلة v2

ψ = N exp(-αr): الصورة الصحيحة.

بالقرب من ب:_CA(R) exp(-αr/R): فكرة الإسقاط الصحيحة.

تقريب لابلاسيان: ≈ -3αr/R] ≈ -3α/R. هذا يقيِّم المعامل المحلي فقط ويتجاهل الحد 1/ص.

الاستنتاج F ∝ 1/R² صحيح فيزيائيًا، لكن الاستنتاج غير مكتمل.

الاشتقاق المصحح

ψ = N exp(-αr): لم يتغير.

بالقرب من B:_CA(R) exp(-αr/R): يتم الاحتفاظ به باعتباره الإسقاط المحلي الفعال.

\(\nabla^2(e^{-\alpha r/R})=e^{-\alpha r/R}\left(\frac{\alpha^2}{R^2}-\frac{2\alpha}{Rr}\right)\)

يؤدي الاشتقاق الكامل إلى تكامل المُشغِّل على الدالة الموجية B ويحصل على

\(V(R)=-\frac{Ke^{-\alpha R}}{R}\) \(F(R)=- \frac{K(1+\ألفا R)e^{-\ألفا R}}{R^2}\)

خاتمة الورقة البحثية صحيحة – كان الاستنباط بحاجة إلى استكمال.

يصل BeeTheory v2 إلى الإجابة الفيزيائية الصحيحة من خلال الحدس الصحيح، ولكن يجب إكمال تقريب أحادي القطب من خلال الاحتفاظ بالحد 1/r في لابلاسيان والتكامل على الدالة الموجية للجسيم الثاني.

المراجع

دوتيرتر، إكس. – نظرية النحلة™: نمذجة الجاذبية على أساس الموجة، BeeTheory.com v2، 2023.
يوكاوا، هـ. – حول تفاعل الجسيمات الأولية، Proc. Phys.-Math. Soc. Japan 17, 48, 1935.
أبراموفيتز، م.، ستيغون، أ. أ. – كتيب الدوال الرياضية، دوفر، 1972.
جاكسون، ج. د. – الديناميكا الكهربائية الكلاسيكية، الطبعة الثالثة، وايلي، 1999.
جريفيث، د. ج. – مقدمة في ميكانيكا الكم، الطبعة الثانية، بيرسون، 2005.

موقع BeeTheory.com – استكشاف الجاذبية من خلال فيزياء الكم القائمة على الموجات