الإطار العددي — الورقة التفاعلية الشاملة

🔬 المحاكاة التفاعلية (Interactive Simulation)

استخدم شريط التمرير أدناه لتغيير قيمة القطع \(x_p\). تم توسيع نطاق الرسم البياني إلى \(x=50\) لتأكيد أن دالة التجميع \(A(x; x_p)\) مستمرة في الصعود ولا تمتلك قمة عند \(x=42\).

1. صفر ريمان والقطع المقترح

تعريف (القطع الصحيح المقترح). ليكن \(\gamma_2 \approx 21.022039\) هو الإحداثي الصادي لثاني صفر غير بديهي موجب لدالة زيتا الريمانية. نعرّف معامل النموذج:

تنبيه اصطلاحي. اختيار \(x_p = 21\) اختيار نمذجي حرّ، مستوحى من \(\lfloor \gamma_2 \rfloor\) لأغراض التقريب العددي فقط. ليس اشتقاقًا من دالة زيتا، ولا يدّعي النموذج وجود علاقة رياضية بنيوية بين أصفار ريمان والدالة \(Z\).

2. تحليل الإشارة (البرهان الرياضي)

من أجل \(x_p>0\)، تحقق الدالة الحقيقية الشروط التالية:

• \(Z(x;x_p) < 0\) من أجل \(1 < x < 2\).
• \(Z(1;x_p) = Z(2;x_p) = 0\) (جذرا الدالة في \([1,\infty)\)).
• \(Z(x;x_p) > 0\) من أجل \(x > 2\).

3. دالة التجميع

نعرّف دالة التجميع رياضياً على المجال \(x \ge 1\) بالعلاقة: \(A(x) := \int_1^x Z(t;21)\,dt\). وفقاً للنظرية الأساسية للتفاضل والتكامل، \(A'(x) = Z(x;21)\).

بما أن \(Z(x)\) سالبة على \((1, 2)\) وموجبة على \((2, \infty)\)، فإن \(A(x)\) تتناقص على \((1, 2)\)، وتصل إلى قيمة صغرى مطلقة (Global Minimum) عند \(x=2\)، ثم تتزايد على \((2, \infty)\).

سلوك \(A(x)\) عند اللانهاية

التكامل \(A_\infty := \int_1^\infty Z(t;21)\,dt\) متقارب تقاربًا مطلقًا بسبب العامل \(e^{-t/21}\). بتحويل \(t = 2\pi/x\) نحصل على:

بتقسيم التكامل على \((0,\pi)\) و\((\pi, 2\pi)\) واستخدام التعويض \(t \mapsto 2\pi - t\) في الشطر الثاني:

النتيجة. بما أن \(A(2)<0\)، و\(A'(x) > 0\) على \((2,\infty)\)، و\(A_\infty < 0\)، فإن \(A(x) < 0\) لكل \(x > 1\)، ولا تعود \(A(x)\) إلى الصفر بعد \(x=2\).

4. الامتداد المركب والغلاف التقاربي

لفحص السلوك على الخط الحرج \(s=\tfrac12+it\)، نعرّف امتداداً مركباً صريحاً باستخدام الفرع الرئيسي للوغاريتم \(\operatorname{Log}(s)\):

عندما \(t \to \infty\)، يُعطى الغلاف التقاربي للامتداد المركب بواسطة:

5. تحليل الحساسية وحالة الأدلة

مصدر القيمة المستهدفة. القيمة \(K \approx 6.1353\) هي \(K(21) = 2\pi e^{-1/42}\) نفسها، محسوبة من الصيغة البنائية للنموذج. لذلك فإن «التطابق» بين \(K(21)\) و\(6.1353\) تحقق دائري بالبناء، وليس اختبارًا مستقلًا.

يُقرأ الجدول أدناه بوصفه تحليل حساسية لاختيار \(x_p\) ضمن نطاق محدود (19–22)، لا بوصفه دليلًا على تفوّق \(21\) على أي قيمة أخرى. المقارنة مفيدة لمعرفة درجة انحناء \(K(x_p)\) حول القيمة المختارة، لا لإثبات تفرّدها.

6. امتدادات كونية تخمينية (RiemannOS)

🛠️ 7. الأدوات البرمجية والتحقق العددي (Computational Tools)

لضمان إمكانية إعادة إنتاج النتائج والتحقق من صحة المشتقات، تم تزويد هذا الإطار بالأدوات التفاعلية والبرمجية التالية:

8. الملحق D — التحقق العددي التفصيلي

يقدّم هذا الملحق تفاصيل التحقق العددي المنفصل عن البرهان التحليلي الذي يدعم النتائج، دون أن يحل محل البرهان:

1. إشارة \(A_\infty\): تقدير صريح لباقي الذيل، مع حساب \(A(100)\) بدقة عالية.
2. تقارب \(C(t)\): شيفرة Python لتقييم \(|\mathcal{Z}(1/2+it;21)|\) ومقارنته بالتوقع التقاربي (النتائج في جدول §7).

D.1 تقدير باقي الذيل لـ \(A_\infty\)

لكل \(t \ge 2\)، المتباينة \(|\sin(2\pi/t)| \le 2\pi/t\) صحيحة. بما أن \(\ln t > 0\) على هذا المجال، فإن:

لإيجاد حد للتكامل، نلاحظ أن الدالة \(g(t)=t^{-6}\ln t\) تتناقص لكل \(t \ge 2\)، لأن:

لذلك، لكل \(X \ge 2\):

عند \(X=100\):

D.2 الشيفرة (Python — NumPy / SciPy)

import numpy as np
from scipy.integrate import quad

# D.1 — التكامل A(100) وحدّ الذيل
def Z_real(t, xp=21.0):
    if t <= 1: return 0.0
    return (t**-5) * np.log(t) * np.sin(2*np.pi/t) * np.exp(-t/xp)

A_100, err_A100 = quad(Z_real, 1, 100, limit=500,
                       epsabs=1e-15, epsrel=1e-12)
X = 100.0
tail = 42*np.pi * X**-6 * np.log(X) * np.exp(-X/21.0)
print(f"A(100)      = {A_100:.12f}   (quad err: {err_A100:.1e})")
print(f"Tail bound  = {tail:.3e}")
print(f"A_inf < 0 (numerically supported): {A_100 + err_A100 + tail < 0}")
print()

# D.2 — تقييم |Z(1/2+it;21)|
def Z_cplx(s, xp=21.0):
    return (s**-5) * np.log(s) * np.sin(2*np.pi/s) * np.exp(-s/xp)

K21 = 2*np.pi*np.exp(-1/42)
print(f"K(21) = 2*pi*e^(-1/42) = {K21:.6f}")
print(f"{'t':>8} {'K ln(t)/t^6':>16} {'|Z|':>16} "
      f"{'C(t)':>10} {'|C-K|':>10} {'rel':>10}")
print("-" * 80)
for t in [1e2, 1e3, 1e4, 1e5]:
    z = abs(Z_cplx(0.5 + 1j*t))
    C = z * t**6 / np.log(t)
    err = abs(C - K21); rel = err / K21
    asymp = K21 * np.log(t) / t**6
    print(f"{t:>8.0e} {asymp:>16.4e} {z:>16.4e} "
          f"{C:>10.6f} {err:>10.6f} {rel:>10.2e}")

D.3 نتيجة \(A_\infty\)

التقييم العددي وحدّ الذيل. أعطى التكامل العددي \(A(100) \approx -0.0246463895708\)، وأبلغت scipy.integrate.quad عن تقدير خطأ يقارب \(4.3 \times 10^{-16}\) في التنفيذ المرجعي. أما حدّ الذيل التحليلي عند \(X=100\) فهو أقل من \(5.2 \times 10^{-12}\).

تقدير الخطأ الذي تُرجعه quad ليس حدًا صارمًا مضمونًا؛ لذلك تُعرض هذه النتائج بوصفها تحققًا عدديًا قويًا، لا شهادة احتواء حسابية مضمونة. وبما أن إشارة \(A_\infty < 0\) قد ثبتت بالفعل بالبرهان التحليلي في §3، فإن النتيجة النهائية لا تعتمد على افتراض أن تقدير خطأ التكامل حدّ صارم.

D.4 ملاحظات إعادة الإنتاج

• البيئة: Python 3.9+ مع NumPy وSciPy. التكامل يستخدم scipy.integrate.quad بتفاوتات epsabs=1e-15 وepsrel=1e-12 وlimit=500.
• الدقة: القيم المعروضة تقريبية، وقد تختلف في الرقم الأخير بحسب إصدار المكتبة. الاستنتاجات النوعية (إشارة \(A_\infty\)، اتجاه \(C(t)\)) مستقرة ولا تعتمد على هذا الاختلاف.
• الفرع الرئيسي للوغاريتم: يُستخدم numpy.log للفرع الرئيسي، بجزء تخيلي في \((-\pi, \pi]\)، وهو التعريف المعتمد في §4.
• قابلية التصحيح: أي تصحيح في الصيغ الرياضية يجب أن يظهر في §3 و§4 أولًا، ثم يُعاد اشتقاق الشيفرة؛ لا يُعدَّل الحد العددي استقلالًا عن البرهان.