تفرد BKL

جسم كروي يخضع لديناميكيات فوضوية من نوع BKL (Mixmaster) قريبة من حالة التفرد وفقًا للقواعد المذكورة في المعادلة 35. تمت المحاكاة باستخدام برنامج Mathematica مع قيم ابتدائيةu=13{\displaystyle u={\sqrt {13}}}[ ملاحظة 1 ]

تُعدّ نقطة بيلينسكي -خالاتنيكوف-ليفشيتز ( BKL ) نموذجًا للتطور الديناميكي للكون بالقرب من نقطة التفرد الجاذبي الأولية ، والموصوفة بحل فوضوي غير متناحٍ لمعادلة أينشتاين للمجال الجاذبي . [ 2 ] وفقًا لهذا النموذج، يتذبذب الكون بشكل فوضوي حول نقطة تفرد جاذبي يصبح فيها الزمان والمكان مساويين للصفر، أو بعبارة أخرى، يصبح انحناء الزمكان لانهائيًا. هذه النقطة حقيقية فيزيائيًا بمعنى أنها خاصية ضرورية للحل، وستظهر أيضًا في الحل الدقيق لتلك المعادلات. لا يتم إنشاء هذه النقطة بشكل مصطنع من خلال الافتراضات والتبسيطات التي تقوم عليها الحلول الخاصة الأخرى مثل حلول فريدمان-لوميتر-روبرتسون-ووكر ، وشبه المتناحية، وكاسنر .

سُمي النموذج على اسم مؤلفيه فلاديمير بيلينسكي وإسحاق خالاتنيكوف وإيفجيني ليفشيتز ، الذين كانوا يعملون آنذاك في معهد لانداو للفيزياء النظرية .

تتضمن الصورة التي طورتها شركة BKL عدة عناصر مهمة، وهي:

  • بالقرب من نقطة التفرد، ينفصل تطور الهندسة عند نقاط مكانية مختلفة، بحيث يمكن تقريب حلول المعادلات التفاضلية الجزئية بحلول المعادلات التفاضلية العادية بالنسبة للزمن، وذلك لمعاملات قياس مكانية محددة بشكل مناسب . يُعرف هذا باسم حدسية BKL .
  • بالنسبة لمعظم أنواع المادة، يصبح تأثير حقول المادة على ديناميكيات الهندسة ضئيلاً بالقرب من نقطة التفرد. أو كما قال جون ويلر ، "المادة لا تُهم" بالقرب من نقطة التفرد. افترضت دراسة BKL الأصلية تأثيرًا ضئيلاً لجميع أنواع المادة، لكنهم افترضوا لاحقًا أن "المادة الصلبة" (معادلة الحالة p = ε) المكافئة لحقل قياسي عديم الكتلة يمكن أن يكون لها تأثير مُعدِّل على الديناميكيات بالقرب من نقطة التفرد.
  • تأتي المعادلات التفاضلية العادية التي تصف السلوك التقاربي من فئة من الحلول المتجانسة مكانيًا والتي تشكل ديناميكيات Mixmaster : نموذج تذبذبي وفوضوي معقد يظهر خصائص مماثلة لتلك التي ناقشها BKL.

أصبحت دراسة ديناميكيات الكون في جوار التفرد الكوني مجالًا سريع التطور في الفيزياء النظرية والرياضية الحديثة . يتسم تعميم نموذج BKL للتفرد الكوني في النماذج الكونية متعددة الأبعاد ( من نوع كالوزا-كلاين ) بطابع فوضوي في الزمكانات التي لا يتجاوز بُعدها عشرة أبعاد، بينما في الزمكانات ذات الأبعاد الأعلى، يدخل الكون، بعد عدد محدود من التذبذبات، في نظام انكماش رتيب من نوع كاسنر. [ 3 ] [ 4 ] [ 5 ]

كشف تطور الدراسات الكونية القائمة على نماذج الأوتار الفائقة عن جوانب جديدة للديناميكيات في جوار التفرد. [ 6 ] [ 7 ] [ 8 ] في هذه النماذج، لا تُستحث آليات تغير حقب كاسنر بفعل التفاعلات الجاذبية، بل بفعل تأثير الحقول الأخرى الموجودة. وقد ثبت أن النماذج الكونية القائمة على ستة نماذج رئيسية للأوتار الفائقة، بالإضافة إلى نموذج الجاذبية الفائقة ذي أحد عشر بُعدًا، تُظهر ديناميكيات BKL الفوضوية باتجاه التفرد. كما تم اكتشاف صلة بين النماذج الكونية التذبذبية الشبيهة بـ BKL وفئة فرعية خاصة من جبر لي اللانهائي الأبعاد ، وهي ما يُسمى بجبر كاك-مودي الزائدي . [ 9 ] [ 10 ] [ 11 ]

مقدمة

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

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

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

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

يستحيل إيجاد تكامل عام لجميع الأماكن والأزمنة. مع ذلك، ليس هذا ضروريًا لحل المسألة، إذ يكفي دراسة الحل بالقرب من نقطة التفرد. وهذا من شأنه أن يحل جانبًا آخر من المسألة: خصائص تطور متريّة الزمكان في الحل العام عندما يصل إلى نقطة التفرد الفيزيائي، والتي تُفهم على أنها النقطة التي تصبح عندها كثافة المادة وثوابت موتر انحناء ريمان لانهائية .

وجود تفرد زمني فيزيائي

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

يُعدّ عدد دوال إحداثيات الفضاء المستقلة التي تحتويها الحلول معيارًا لعموميتها. ويشمل ذلك فقط الدوال "المستقلة فيزيائيًا" التي لا يمكن تقليل عددها باختيار أي إطار مرجعي . في الحل العام، يجب أن يكون عدد هذه الدوال كافيًا لتحديد الشروط الابتدائية (توزيع المادة وحركتها، وتوزيع مجال الجاذبية ) تحديدًا كاملًا عند لحظة زمنية مختارة كلحظة ابتدائية. هذا العدد هو أربعة للفضاء الفارغ، وثمانية للفضاء المملوء بالمادة و/أو الإشعاع. [ 12 ] [ 13 ]

أدت أعمال سابقة لمجموعة لانداو ( [ 14 ] [ 15 ] [ 16 ] ، والتي تمت مراجعتها في [ 12 ] ) إلى استنتاج مفاده أن الحل العام لا يحتوي على نقطة تفرد فيزيائية. وقد تم البحث عن فئة أوسع من الحلول التي تحتوي على نقطة تفرد، بشكل أساسي، عن طريق التجربة والخطأ، نظرًا لافتقارنا إلى منهجية منظمة لدراسة معادلات أينشتاين. إن النتيجة السلبية التي تم التوصل إليها بهذه الطريقة ليست مقنعة بحد ذاتها؛ إذ أن وجود حل يتمتع بدرجة كافية من العمومية من شأنه أن يُبطلها، وفي الوقت نفسه يؤكد أي نتائج إيجابية متعلقة بالحل المحدد.

في ذلك الوقت، كان المؤشر الوحيد المعروف لوجود تفرد فيزيائي في الحل العام مرتبطًا بشكل معادلات أينشتاين المكتوبة في إطار متزامن ، أي في إطار يكون فيه الزمن الخاص x₀ = t متزامنًا في جميع أنحاء الفضاء؛ في هذا الإطار، يكون عنصر المسافة المكانية dl منفصلًا عن الفترة الزمنية dt . [ ملاحظة 2 ] معادلة أينشتاين

إن كتابة المعادلات في إطار متزامن تعطي نتيجة يصبح فيها المحدد المتري g صفراً حتماً في زمن محدود بغض النظر عن أي افتراضات حول توزيع المادة. [ 12 ] [ 13 ]

مع ذلك، تم التخلي عن الجهود المبذولة لإيجاد تفرد فيزيائي عام بعد أن اتضح أن التفرد المذكور أعلاه مرتبط بخاصية هندسية محددة للإطار المتزامن: تقاطع إحداثيات خط الزمن . يحدث هذا التقاطع على بعض الأسطح الفائقة المحيطة ، وهي نظائر رباعية الأبعاد للأسطح الكاوية في البصريات الهندسية ؛ وتصبح قيمة g صفرًا عند هذا التقاطع تحديدًا. [ 16 ] لذلك، على الرغم من أن هذا التفرد عام، إلا أنه وهمي وليس فيزيائيًا؛ إذ يختفي عند تغيير الإطار المرجعي. ويبدو أن هذا الأمر قد ثبط عزيمة الباحثين عن مواصلة البحث في هذا الاتجاه.

مرت عدة سنوات قبل أن يعود الاهتمام بهذه المشكلة للظهور مجدداً عندما نشر بنروز ( 1965 ) نظرياته التي ربطت وجود حالة شاذة ذات طبيعة غير معروفة ببعض الافتراضات العامة جداً التي لا علاقة لها باختيار إطار مرجعي. وقد توصل هوكينغ [ 17 ] [ 18 ] وجيروش [ 19 ] لاحقاً إلى نظريات مماثلة (انظر نظريات بنروز-هوكينغ حول الحالات الشاذة ). وقد أدى ذلك إلى إحياء الاهتمام بالبحث عن حلول شاذة. 

حل متجانس معمّم

في فضاء متجانس ومتناحٍ، تُحدد المترية بشكل كامل، ولا يبقى سوى إشارة الانحناء. بافتراض تجانس الفضاء فقط دون أي تناظر إضافي كالتناحٍ، تتوفر حرية أكبر بكثير في اختيار المترية. يتعلق ما يلي بالجزء المكاني من المترية عند لحظة زمنية معينة بافتراض إطار مرجعي متزامن بحيث يكون t هو نفس الزمن المتزامن للفضاء بأكمله.

تخمين BKL

في دراستهما عام 1970، [ 2 ] ذكر BKL أنه مع اقترابنا من نقطة تفرد، تطغى الحدود التي تحتوي على مشتقات زمنية في معادلات أينشتاين على تلك التي تحتوي على مشتقات مكانية . عُرف هذا منذ ذلك الحين باسم حدسية BKL ، ويشير إلى أن المعادلات التفاضلية الجزئية لأينشتاين تُقارب جيدًا بالمعادلات التفاضلية العادية ، ومن ثم تصبح ديناميكيات النسبية العامة محلية وتذبذبية بشكل فعلي. ويُقارب التطور الزمني للحقول عند كل نقطة مكانية جيدًا بواسطة النماذج الكونية المتجانسة في تصنيف بيانكي.

بفصل مشتقات الزمان والمكان في معادلات أينشتاين، على سبيل المثال، بالطريقة المستخدمة في تصنيف الفضاءات المتجانسة ، ثم جعل الحدود التي تحتوي على مشتقات مكانية تساوي صفرًا، يمكن تعريف ما يُسمى بالنظرية المبتورة للنظام (المعادلات المبتورة). [ 20 ] عندئذٍ، يمكن جعل فرضية BKL أكثر تحديدًا:

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

تخمين قوي : مع اقتراب التفرد، تقترب معادلات أينشتاين من معادلات النظرية المقتطعة، بالإضافة إلى أن حلول المعادلات الكاملة يتم تقريبها بشكل جيد بواسطة حلول المعادلات المقتطعة.

في البداية، بدت فرضية BKL معتمدة على الإحداثيات وغير معقولة إلى حد ما. فعلى سبيل المثال، أشار بارو وتيبلر [ 21 ] [ 22 ] ، ضمن عشرة انتقادات لدراسات BKL، إلى عدم ملاءمة (برأيهما) اختيار الإطار المتزامن كوسيلة لفصل مشتقات الزمان والمكان. وقد أُعيدت صياغة فرضية BKL أحيانًا في الأدبيات العلمية على أنها تنص على أن مشتقات الزمان فقط هي المهمة بالقرب من التفرد. إلا أن هذا الادعاء، إذا أُخذ حرفيًا، خاطئ أو مضلل في أحسن الأحوال، إذ كما هو موضح في تحليل BKL نفسه، لا يمكن إهمال التدرجات المكانية لموتر القياس في الحلول العامة لجاذبية أينشتاين الخالصة في أربعة أبعاد مكانية-زمانية، بل إنها تلعب دورًا حاسمًا في ظهور النظام التذبذبي. مع ذلك، توجد صياغات جديدة لنظرية أينشتاين بدلالة متغيرات جديدة تتضمن التدرجات ذات الصلة، كما في متغيرات أشتيكار، والتي يكون فيها القول حول الدور المهيمن للمشتقات الزمنية صحيحًا. [ 20 ] صحيح أنه عند كل نقطة مكانية نحصل على وصف فعال للتفرد بدلالة نظام ديناميكي محدود الأبعاد موصوف بمعادلات تفاضلية عادية بالنسبة للزمن، لكن التدرجات المكانية تدخل هذه المعادلات بشكل غير بديهي.

أظهرت تحليلات لاحقة أجراها عدد كبير من الباحثين إمكانية تحديد فرضية BKL بدقة، ويوجد الآن كمٌّ هائل من الأدلة العددية والتحليلية التي تدعمها. [ 23 ] من الإنصاف القول إننا ما زلنا بعيدين عن إثبات هذه الفرضية القوية. ولكن تم إحراز تقدم ملحوظ في النماذج الأبسط. على وجه الخصوص، بيّن بيرغر، وغارفينكل، ومونكريف، وإيزنبرغ، وويفر، وآخرون أنه في فئة من النماذج، مع اقترابنا من التفرد، تقترب حلول معادلات أينشتاين الكاملة للمجال من الحلول "المُهيمن عليها حد السرعة" (المُقتطعة) التي تم الحصول عليها بإهمال المشتقات المكانية. [ 23 ] [ 24 ] [ 25 ] [ 26 ] [ 27 ] أظهر أندرسون وريندال [ 28 ] أنه بالنسبة للجاذبية المقترنة بحقل قياسي عديم الكتلة أو مائع صلب، فإنه لكل حل للمعادلات المقتطعة يوجد حل لمعادلات الحقل الكاملة يتقارب مع الحل المقتطع عند الاقتراب من التفرد، حتى في غياب التناظرات. وقد عُممت هذه النتائج لتشمل أيضًا حقول القياس من النوع p . [ 29 ] في هذه النماذج المقتطعة، تكون الديناميكيات أبسط، مما يسمح بصياغة دقيقة للفرضية التي يمكن إثباتها. في الحالة العامة، يأتي أقوى دليل حتى الآن من التطورات العددية. بدأ بيرغر ومونكريف [ 30 ] برنامجًا لتحليل التفردات الكونية العامة. بينما ركز العمل الأولي على حالات التناظر المختزل، [ 31 ] قام غارفينكل مؤخرًا [ 32 ] بإجراء تطور عددي لزمكانات بلا تناظرات، حيث يظهر سلوك الميكس ماستر بوضوح. أخيرًا، جاء دعم إضافي للفرضية من دراسة عددية لسلوك حقول الاختبار بالقرب من تفرد ثقب أسود شوارزشيلد. [ 33 ]

حل كاسنر

ديناميكيات مقاييس كاسنر ( المعادلة 2) في الإحداثيات الكروية نحو التفرد. معامل ليفشيتز-خالاتنيكوف هو u = 2 (1/ u = 0.5)، والإحداثي r هو 2πα ( 1 / u )τ، حيث τ هو الزمن اللوغاريتمي: τ = ln t . [ ملاحظة 3 ] الانكماش على طول المحاور خطي وغير متناحٍ (لا فوضوي).

يبدأ منهج BKL للفضاءات المتجانسة غير المتناحية (على عكس المتناحية) بتعميم حل خاص دقيق اشتقه كاسنر [ 34 ] لحقل في الفراغ، حيث يكون الفضاء متجانسًا وله مقياس إقليدي يعتمد على الزمن وفقًا لمقياس كاسنر.

( dl هو عنصر الخط ؛ dx و dy و dz هي إزاحات متناهية الصغر في الأبعاد المكانية الثلاثة ، و t هو الفترة الزمنية المنقضية منذ لحظة ابتدائية t₀ = 0 ). هنا، p₁ و p₂ و p₃ هي أي ثلاثة أعداد تحقق شروط كاسنر التالية

بسبب هذه العلاقات، يكون عدد واحد فقط من الأعداد الثلاثة مستقلاً ( معادلتان بثلاثة مجاهيل ). لا تتطابق الأعداد الثلاثة أبدًا؛ يتطابق عددان فقط في مجموعات القيم.(-13،23،23){\textstyle (-{\frac {1}{3}},{\frac {2}{3}},{\frac {2}{3}})}و(0، 0، 1). [ ملاحظة 4 ] في جميع الحالات الأخرى، تكون الأرقام مختلفة، أحدها سالب والآخران موجبان. وقد تم إثبات ذلك جزئيًا بتربيع طرفي المعادلة 3 ( الشرط الأول ) ثم فك التربيع.

(ص1+ص2+ص3)2=(ص12+ص22+ص32)+(2ص1ص2+2ص2ص3+2ص1ص3)=1{\displaystyle \left(p_{1}+p_{2}+p_{3}\right)^{2}=\left(p_{1}^{2}+p_{2}^{2}+p_{3}^{2}\right)+\left(2p_{1}p_{2}+2p_{2}p_{3}+2p_{1}p_{3}\right)=1}

على المدى(ص12+ص22+ص32){\displaystyle \left(p_{1}^{2}+p_{2}^{2}+p_{3}^{2}\right)}يساوي 1 وفقًا للشرط الثاني في المعادلة وبالتالي يجب أن يكون الحد الذي يحتوي على نواتج مختلطة مساويًا للصفر. وهذا ممكن إذا كان أحد الحدود p1 أو p2 أو p3 على الأقل سالبًا.

إذا تم ترتيب الأرقام بترتيب تصاعدي، p 1 < p 2 < p 3 ، فإنها تتغير في الفترات (الشكل 4).

رسم بياني لـ p1 و p2 و p3 بدلالة المتغير 1/ u . تتزايد قيم p1 ( u ) و p3 ( u ) بشكل رتيب، بينما تتناقص قيمة p2 ( u ) بشكل رتيب مع u .

تُقابل معادلة كاسنر ( 2) فضاءً مسطحًا متجانسًا ولكنه غير متناحٍ، حيث تزداد جميع الأحجام مع الزمن بحيث تزداد المسافات الخطية على طول المحورين y و z بينما تتناقص المسافة على طول المحور x . تُسبب اللحظة t = 0 حالةً شاذةً في الحل؛ ولا يمكن تجنب هذه الحالة الشاذة في المقياس عند t = 0 بأي تحويل إطار مرجعي. عند هذه الحالة الشاذة، تؤول ثوابت موتر الانحناء رباعي الأبعاد إلى اللانهاية. ويُستثنى من ذلك الحالة p₁ = p₂ = 0، p₃ = 1؛ إذ تُقابل هذه القيم فضاءً زمكانيًا مسطحًا: يُحوّل التحويل t = sh z = ζ، t = ch z = τ مقياس كاسنر ( المعادلة 2 ) إلى مقياس غاليلي .

تُحدد BKL الأعداد p1 و p2 و p3 بدلالة مُعامل واحد مستقل ( حقيقي ) u (معامل Lifshitz-Khalatnikov [ 35 ] ) على النحو التالي

تبدو معلمات مؤشر كاسنر غامضة حتى نتأمل في القيدين المفروضين على المؤشرات (المعادلة 3) . يحدد هذان القيدان المقياس الكلي للمؤشرات بحيث لا تتغير إلا نسبها . من الطبيعي اختيار إحدى هذه النسب كمعلمة جديدة، وهو ما يمكن فعله بست طرق مختلفة. على سبيل المثال، باختيار u = u₃² = p₃ / p₂ ، يصبح من السهل التعبير عن جميع النسب الست الممكنة بدلالة هذه النسبة. بحذف p₃ = u₂² أولًا ، ثم باستخدام القيد الخطي لحذف p₁ = 1p₂u₂² = 1 − (1 + u ) p₂ ، يختزل القيد التربيعي إلى معادلة تربيعية في p₂

بجذور p₂ = 0 (بديهي) و p₂ = ( 1 + u) / (1 + u + u₂)، ومنها يتم الحصول على p₁ و p₃ بالتعويض العكسي . يمكن تعريف ستة معاملات من هذا النوع uₑ = pₐ / pₑ ، حيث pₑpₑ pₐ عندما يكون ( c , b , a ) تبديلاً دورياً لـ ( 1 , 2 , 3 ) . [ 36 ]

يتم الحصول على جميع القيم المختلفة لـ p1 و p2 و p3 ، مرتبة كما سبق ، عندما تكون قيمة u ضمن النطاق u ≥ 1. أما القيم u < 1 فتُوضع ضمن هذا النطاق وفقًا لـ

في الحل المعمم، ينطبق الشكل المقابل للمعادلة 2 فقط على المقياس التقاربي (المقياس القريب من النقطة الشاذة t = 0)، وعلى الحدود الرئيسية لتوسيع متسلسلته بدلالة قوى t . في إطار المرجع المتزامن، يُكتب على شكل المعادلة 1 مع عنصر مسافة مكانية .

أين

تُحدد المتجهات ثلاثية الأبعاد l و m و n الاتجاهات التي تتغير بها المسافة المكانية مع الزمن وفقًا لقوانين القوى في المعادلة 8. هذه المتجهات ، بالإضافة إلى الأعداد p<sub> l</sub> و p<sub> m</sub> و p<sub> n </sub>، المرتبطة، كما في السابق، بالمعادلة 3 ، هي دوال لإحداثيات المكان. لا تُرتَّب القوى p<sub> l</sub> و p<sub> m</sub> و p<sub> n</sub> بترتيب تصاعدي، مع الاحتفاظ بالرموز p <sub>1</sub> و p<sub> 2</sub> و p<sub> 3 </sub> للأعداد في المعادلة 5 التي تبقى مرتبة بترتيب تصاعدي. مُحدِّد المقياس في المعادلة 7 هو

حيث v = l [ mn ]. من الملائم إدخال الكميات التالية [ ملاحظة 5 ]

إن المقياس المكاني في المعادلة 7 غير متناحٍ لأن قوى t في المعادلة 8 لا يمكن أن تتساوى قيمها. عند الاقتراب من نقطة التفرد عند t = 0، تتناقص المسافات الخطية في كل عنصر مكاني في اتجاهين وتزداد في الاتجاه الثالث. ويتناقص حجم العنصر بما يتناسب مع t .

يتم إدخال مقياس كاسنر في معادلات أينشتاين عن طريق استبدال موتر المقياس γ αβ من المعادلة 7 دون تحديد مسبق لاعتماد a و b و c على t : [ ملاحظة 2 ]

ϰαβ=2أ˙ألαلβ+2ب˙بمαمβ+2ج˙جنαنβ{\displaystyle \varkappa _{\alpha }^{\beta }={\frac {2{\dot {a}}}{a}}l_{\alpha }l^{\beta }+{\frac {2{\dot {b}}}{b}}m_{\alpha }m^{\beta }+{\frac {2{\dot {c}}}{c}}n_{\alpha }n^{\beta }}

حيث تشير النقطة فوق الرمز إلى التفاضل بالنسبة للزمن. تأخذ معادلة أينشتاين (المعادلة 11) الشكل التالي:

جميع حدودها من الرتبة الثانية بالنسبة للكمية الكبيرة (عند t → 0) 1/ t . في معادلات أينشتاين (المعادلة 12) ، لا تظهر الحدود من هذه الرتبة إلا من الحدود المشتقة زمنيًا. إذا لم تتضمن مركبات P αβ حدودًا من رتبة أعلى من الثانية، فإن

حيث تشير المؤشرات l و m و n إلى مكونات الموتر في الاتجاهات l و m و n . [ 12 ] هذه المعادلات بالإضافة إلى المعادلة 14 تعطي التعبيرات المعادلة 8 مع قوى تحقق المعادلة 3 .

مع ذلك، فإن وجود قوة سالبة واحدة بين القوى الثلاث p <sub>l</sub> و p<sub> m</sub> و p<sub> n</sub> يؤدي إلى ظهور حدود من P <sub>αβ </sub> برتبة أكبر من t <sup>-2 </sup> . إذا كانت القوة السالبة هي p <sub>l </sub> ( حيث p <sub> l</sub> = p <sub>1 </sub> < 0)، فإن P <sub>αβ</sub> تحتوي على دالة الإحداثيات λ، وتصبح المعادلة 12 كما يلي:

هنا، الحدود الثانية من الرتبة t −2( p m + p np l ) حيث p m + p np l = 1 + 2 | p l | > 1. [ ملاحظة 6 ] لإزالة هذه الحدود واستعادة المعادلة المترية 7 ، من الضروري فرض الشرط λ = 0 على دوال الإحداثيات.

تحتوي معادلات أينشتاين الثلاث المتبقية (المعادلة 13) على مشتقات زمنية من الرتبة الأولى فقط لموتر القياس. وهي تُعطي ثلاث علاقات مستقلة عن الزمن يجب فرضها كشروط ضرورية على دوال الإحداثيات في المعادلة 7. هذا، بالإضافة إلى الشرط λ = 0، يُشكّل أربعة شروط. تربط هذه الشروط عشر دوال إحداثيات مختلفة: ثلاثة مركبات لكل من المتجهات l و m و n ، ودالة واحدة في قوى t (أي من الدوال p <sub>l</sub> و p<sub> m</sub> و p<sub> n</sub> ، والتي تُقيّدها الشروط في المعادلة 3 ). عند حساب عدد الدوال الفيزيائية العشوائية، يجب مراعاة أن النظام المتزامن المستخدم هنا يسمح بتحويلات عشوائية مستقلة عن الزمن للإحداثيات المكانية الثلاث. لذلك، يحتوي الحل النهائي على 10 - 4 - 3 = 3 دوال فيزيائية عشوائية، وهو أقل بواحدة مما هو مطلوب للحل العام في الفراغ.

لا ينقص إدخال المادة من درجة العمومية التي تم التوصل إليها عند هذه النقطة؛ إذ تُدرج المادة في المعادلة المترية 7 ، وتُضيف أربع دوال إحداثية جديدة ضرورية لوصف التوزيع الأولي لكثافتها ومكونات سرعتها الثلاثة. وهذا يُتيح تحديد تطور المادة انطلاقًا من قوانين حركتها في مجال جاذبية مُعطى مسبقًا، وهي المعادلات الهيدروديناميكية .

حيث uᵢ هي السرعة رباعية الأبعاد، و ε وσ هما كثافتا الطاقة والإنتروبيا للمادة (انظر [ 37 ] و [ 38 ] ، وكذلك [ 39 ] ، وللمزيد من التفاصيل انظر [ 40 ] ). بالنسبة لمعادلة الحالة النسبية الفائقة p = ε/3، فإن الإنتروبيا σ ~ ε 1/4 . الحدود الرئيسية في المعادلتين 17 و 18 هي تلك التي تحتوي على مشتقات زمنية . من المعادلة 17 والمكونات المكانية للمعادلة 18، ​​نحصل على

ت(-زu0ε34)=0، 4εuαت+uαεت=0،{\displaystyle {\frac {\partial }{\partial t}}\left({\sqrt {-g}}u_{0}\varepsilon ^{\frac {3}{4}}\right)=0,\ 4\varepsilon \cdot {\frac {\partial u_{\alpha }}{\partial t}}+u_{\alpha }\cdot {\frac {\partial \varepsilon }{\partial t}}=0,}

مما أدى إلى

حيث تمثل 'const' كميات مستقلة عن الزمن. بالإضافة إلى ذلك، من المتطابقة u i u i = 1، نحصل على (لأن جميع المكونات المتغيرة لـ u α من نفس الرتبة)

u02uنuن=uن2ج2،{\displaystyle u_{0}^{2}\approx u_{n}u^{n}={\frac {u_{n}^{2}}{c^{2}}},}

حيث يمثل u <sub>n</sub> مركبة السرعة على طول اتجاه n المرتبطة بأعلى قوة (موجبة) لـ t (بافتراض أن p <sub>n</sub> = p<sub> 3 </sub>). من العلاقات المذكورة أعلاه، يتبين أن

أو

يمكن استخدام المعادلات المذكورة أعلاه للتأكد من أن مكونات موتر الإجهاد والطاقة والزخم للمادة الموجودة في الجانب الأيمن من المعادلات

R00=تي00-12تي، Rαβ=تيαβ-12دلتاαβتي،{\displaystyle R_{0}^{0}=T_{0}^{0}-{\frac {1}{2}}T,\ R_{\alpha }^{\beta }=T_{\alpha }^{\beta }-{\frac {1}{2}}\delta _{\alpha }^{\beta }T,}

في الواقع، هي أقل رتبة بمقدار 1/ t من الحدود الرئيسية في طرفيها الأيسر. في المعادلاتRα0=تيα0{\displaystyle R_{\alpha }^{0}=T_{\alpha }^{0}}إن وجود المادة لا يؤدي إلا إلى تغيير العلاقات المفروضة على وظائف الإحداثيات المكونة لها. [ 12 ]

إن حقيقة أن قيمة ε تصبح لانهائية وفقًا للمعادلة 21 تؤكد أنه في حل المعادلة 7، نتعامل مع حالة شاذة فيزيائية عند جميع قيم القوى p1 و p2 و p3 باستثناء ( 0 ، 0 ، 1 ) . بالنسبة لهذه القيم الأخيرة، تكون الحالة الشاذة غير فيزيائية ويمكن إزالتها بتغيير الإطار المرجعي.

تنشأ النقطة الشاذة الوهمية الموافقة للقوى (0، 0، 1) نتيجةً لعبور إحداثيات خط الزمن فوق سطح بؤري ثنائي الأبعاد . وكما أشير في المرجع [ 12 ] ، يمكن دائمًا اختيار إطار مرجعي متزامن بحيث يحدث هذا العبور الحتمي لخط الزمن تحديدًا على هذا السطح (بدلًا من سطح كاوستيك ثلاثي الأبعاد). لذلك، يجب أن يوجد حلٌّ يتضمن هذه النقطة الشاذة الوهمية المتزامنة للفضاء بأكمله، مع مجموعة كاملة من الدوال الاختيارية اللازمة للحل العام. بالقرب من النقطة t = 0، يسمح ذلك بتوسع منتظم بقوى صحيحة لـ t . لتحليل هذه الحالة، انظر المرجع [ 41 ] .

نمط التذبذب نحو التفرد

الحل العام، بحكم تعريفه، مستقر تمامًا؛ وإلا لما وُجد الكون. أي اضطراب يُعادل تغييرًا في الشروط الابتدائية في لحظة زمنية معينة؛ ولأن الحل العام يسمح بشروط ابتدائية اختيارية، فإن الاضطراب لا يستطيع تغيير طبيعته. من هذا المنظور، فإن الشروط الأربعة المفروضة على دوال الإحداثيات في معادلة الحل 7 هي من أنواع مختلفة: ثلاثة شروط تنشأ من المعادلات.Rα0{\displaystyle R_{\alpha }^{0}}تُعتبر القيم λ = 0 "طبيعية"؛ فهي نتيجة لبنية معادلات أينشتاين. مع ذلك، فإن الشرط الإضافي λ = 0، الذي يُسبب فقدان دالة مشتقة واحدة، هو من نوع مختلف تمامًا: إذ يُمكن أن يُخلّ عدم الاستقرار الناتج عن الاضطرابات بهذا الشرط. ويجب أن يُؤدي تأثير هذا الاضطراب إلى نقل النموذج إلى نمط آخر أكثر عمومية. ولا يُمكن اعتبار هذا الاضطراب صغيرًا: فالانتقال إلى نمط جديد يتجاوز نطاق الاضطرابات الصغيرة جدًا.

يُظهر تحليل سلوك النموذج تحت تأثير الاضطراب، الذي أجراه BKL، نمطًا تذبذبيًا معقدًا عند الاقتراب من نقطة التفرد. [ 2 ] [ 42 ] [ 43 ] [ 44 ] لم يتمكنوا من تقديم جميع تفاصيل هذا النمط في الإطار العام. ومع ذلك، فقد شرح BKL أهم خصائص وطبيعة الحل في نماذج محددة، مما يسمح بدراسة تحليلية معمقة.

تعتمد هذه النماذج على مقياس فضاء متجانس من نوع معين. إن افتراض تجانس الفضاء دون أي تناظر إضافي يتيح حرية كبيرة في اختيار المقياس. تُصنَّف جميع الفضاءات المتجانسة الممكنة (وإن كانت غير متناحية) وفقًا لبيانكي إلى عدة أنواع من بيانكي (من النوع الأول إلى التاسع) . [ 45 ] (انظر أيضًا الحل المتجانس المعمم ). يقتصر بحث BKL على فضاءات النوعين الثامن والتاسع من بيانكي.

إذا كان المقياس على شكل المعادلة 7 ، فإنه لكل نوع من أنواع الفضاءات المتجانسة توجد علاقة دالية بين متجهات الإسناد l و m و n وإحداثيات الفضاء. لا يهم الشكل المحدد لهذه العلاقة. المهم هو أنه بالنسبة للفضاءات من النوعين الثامن والتاسع، فإن الكميات λ و μ و ν (المعادلة 10) ثوابت، بينما جميع نواتج الضرب "المختلطة" l rot m و l rot n و m rot l ، إلخ ، تساوي أصفارًا. بالنسبة للفضاءات من النوع التاسع، فإن الكميات λ و μ و ν لها نفس الإشارة، ويمكن كتابة λ = μ = ν = 1 (لا يؤثر تغيير إشارة الثوابت الثلاثة في آن واحد). أما بالنسبة للفضاءات من النوع الثامن، فإن لثابتين إشارة معاكسة لإشارة الثابت الثالث؛ ويمكن كتابة، على سبيل المثال، λ = − 1 و μ = ν = 1. [ ملاحظة 7 ]

وبالتالي، تقتصر دراسة تأثير الاضطراب على "نمط كاسنر" على دراسة تأثير الحدود التي تحتوي على λ في معادلات أينشتاين. وتُعدّ فضاءات النوعين الثامن والتاسع أنسب النماذج لهذه الدراسة. ولأنّ الكميات الثلاث λ و μ و ν في هذين النوعين من فضاءات بيانكي لا تساوي الصفر، فإنّ الشرط λ = 0 لا يتحقق بغض النظر عن اتجاه اعتماد l و m و n على الزمن وفقًا لقانون القوة السالب .

معادلات أينشتاين لنماذج الفضاء من النوع الثامن والنوع التاسع هي [ 46 ] [ ملاحظة 2 ]

(المكونات المتبقية)Rل0{\displaystyle R_{l}^{0}}،Rم0{\displaystyle R_{m}^{0}}،Rن0{\displaystyle R_{n}^{0}}،Rلم{\displaystyle R_{l}^{m}}،Rلن{\displaystyle R_{l}^{n}}،Rمن{\displaystyle R_{m}^{n}}(أصفار متطابقة). تحتوي هذه المعادلات على دوال للزمن فقط؛ وهذا شرط يجب تحقيقه في جميع الفضاءات المتجانسة. هنا، المعادلتان 22 و 23 دقيقتان، ولا تعتمد صحتهما على مدى قربنا من نقطة التفرد عند t = 0. [ ملاحظة 8 ]

تأخذ المشتقات الزمنية في المعادلتين 22 و 23 شكلاً أبسط إذا تم استبدال а و b و c بلوغاريتماتها α و β و γ:

باستبدال المتغير t بالمتغير τ وفقًا لما يلي:

ثم (تشير الرموز السفلية إلى التفاضل بواسطة τ):

بجمع المعادلات 26 واستبدال المجموع (α + β + γ) τ τ في الطرف الأيسر وفقًا للمعادلة 27 ، نحصل على معادلة تحتوي فقط على المشتقات الأولى وهي التكامل الأول للنظام 26 :

تؤدي هذه المعادلة دور شرط مُلزم مفروض على الحالة الابتدائية للمعادلة 26. يُعد نمط كاسنر (المعادلة 8) حلاً للمعادلة 26 عند إهمال جميع الحدود في الطرف الأيمن. لكن لا يمكن أن يستمر هذا الوضع (عند t → 0) إلى ما لا نهاية، لأنه من بين هذه الحدود، توجد دائمًا حدود تتزايد. بالتالي، إذا كانت القوة السالبة في الدالة a ( t ) (حيث p<sub> l </sub> = p<sub> 1 </sub>)، فإن اضطراب نمط كاسنر سينشأ من الحدود λ<sub> 2</sub> a<sub> 4</sub> ؛ أما بقية الحدود فستتناقص مع تناقص t . إذا اقتصرت الحدود المتزايدة فقط على الطرف الأيمن من المعادلة 26 ، نحصل على النظام التالي:

(قارن المعادلة 16 ؛ حيث تم استبدال λ 2 = 1 أدناه). يجب أن يصف حل هذه المعادلات تطور المقياس من الحالة الابتدائية، والتي يتم وصفها بالمعادلة 8 مع مجموعة معينة من القوى (مع p l < 0)؛ لنفترض أن p l = p 1 ، و p m = p 2 ، و p n = p 3 بحيث

ثم

حيث Λ ثابت. تُعاد تعريف الشروط الابتدائية للمعادلة 29 على النحو التالي:

يمكن تكامل المعادلات 29 بسهولة؛ والحل الذي يحقق الشرط في المعادلة 32 هو

حيث أن b 0 و c 0 ثابتان آخران.

من السهل ملاحظة أن السلوك التقاربي للدوال في المعادلة 33 عند t → 0 هو المعادلة 30. والتعبيرات التقاربية لهذه الدوال والدالة t (τ) عند τ → −∞ هي [ ملاحظة 9 ]

أهـ-Λص1τ، بهـΛ(ص2+2ص1)τ، جهـΛ(ص3+2ص1)τ، تهـΛ(1+2ص1)τ.{\displaystyle a\sim e^{-\Lambda p_{1}\tau },\ b\sim e^{\Lambda (p_{2}+2p_{1})\tau },\ c\sim e^{\Lambda (p_{3}+2p_{1})\tau },\ t\sim e^{\Lambda (1+2p_{1})\tau }.}

عند التعبير عن a و b و c كدوال لـ t ، نحصل على

أين

ثم

يوضح ما سبق أن الاضطراب يعمل بطريقة تُغير نمط كاسنر بنمط آخر، وفي هذه العملية ينقلب الأس السالب لـ t من الاتجاه l إلى الاتجاه m : فإذا كان p'l < 0 سابقًا، فإنه يصبح الآن p'm < 0. خلال هذا التغيير ، تمر الدالة a ( t ) بقيمة عظمى، وتمر الدالة b ( t ) بقيمة صغرى؛ b ، التي كانت متناقصة سابقًا، تزداد الآن؛ و a، التي كانت متزايدة ، تصبح متناقصة؛ و c ( t ) المتناقصة تستمر في التناقص. الاضطراب نفسه (λ²a⁴α في المعادلة 29 )، الذي كان متزايدًا سابقًا، يبدأ الآن في التناقص والتلاشي. وبالمثل ، يؤدي التطور اللاحق إلى زيادة في الاضطراب من الحدود التي تحتوي على μ² ( بدلاً من λ² ) في المعادلة 26 ، ثم تغيير نمط كاسنر، وهكذا.

من الملائم كتابة قاعدة استبدال القوة (المعادلة 35) بمساعدة المعادلة 5 الخاصة بالمعاملات :

تبقى القوة الأكبر من بين القوتين الموجبتين موجبة.

يُطلق BKL على هذا الانقلاب في القوة السالبة بين الاتجاهين اسم حقبة كاسنر . ويكمن مفتاح فهم طبيعة تطور المقياس عند الاقتراب من التفرد في عملية تناوب حقبة كاسنر هذه تحديدًا، مع انقلاب القوى p <sub>l</sub> و p <sub> m</sub> و p<sub> n</sub> وفقًا للقاعدة المذكورة في المعادلة 37 .

تستمر التناوبات المتتالية الموضحة في المعادلة 37، مع انعكاس القوة السالبة p1 بين الاتجاهين l و m (حقب كاسنر)، باستنزاف الجزء الكامل من القيمة الابتدائية u حتى اللحظة التي تصبح فيها u < 1. تتحول القيمة u < 1 إلى u > 1 وفقًا للمعادلة 6 ؛ في هذه اللحظة ، تكون القوة السالبة إما pl أو pm ، بينما تصبح pn أصغر عددين موجبين ( pn = p2 ). ثم تقوم سلسلة حقب كاسنر التالية بانعكاس القوة السالبة بين الاتجاهين n و l أو بين n و m . عند قيمة ابتدائية عشوائية ( غير نسبية ) لـ تستمر عملية التناوب هذه بلا حدود. [ ملاحظة 10 ]

في الحل الدقيق لمعادلات أينشتاين، تفقد القوى p <sub>l</sub> و p <sub> m</sub> و p<sub> n</sub> معناها الدقيق الأصلي. يُدخل هذا الأمر نوعًا من "الغموض" في تحديد هذه الأعداد (وبالتالي، في تحديد المعامل u )، وهو ما يجعل تحليل أي قيم محددة (مثل القيم النسبية ) لـ u عديم الجدوى، على الرغم من صغره. لذلك، فإن القوانين التي تتعلق بالقيم غير النسبية العشوائية لـ u هي فقط التي تحمل معنىً خاصًا.

تُسمى الفترات الأطول التي تتذبذب فيها مقاييس المسافات المكانية على طول محورين بينما تتناقص المسافات على طول المحور الثالث بشكل رتيب، بالعصور ؛ وتتناقص الأحجام وفقًا لقانون قريب من t . عند الانتقال من عصر إلى آخر، ينعكس اتجاه تناقص المسافات بشكل رتيب من محور إلى آخر. يكتسب ترتيب هذه الانتقالات طابعًا تقاربيًا لعملية عشوائية . ويُعدّ الترتيب العشوائي نفسه سمة مميزة لتناوب أطوال العصور المتتالية (يقصد BKL بطول العصر عدد حقب كاسنر التي يحتويها العصر، وليس فترة زمنية).

لكل حقبة ( الحقبة s ) سلسلة من قيم المعامل u بدءًا من القيمة الأكبر،uالأعلى(s){\displaystyle u_{\max }^{(s)}}ومن خلال القيمuالأعلى(s){\displaystyle u_{\max }^{(s)}}-1،uالأعلى(s){\displaystyle u_{\max }^{(s)}}-2، ...، وصولاً إلى أصغرها،uمين(s){\displaystyle u_{\min }^{(s)}}< 1. ثم

أي أن k ( s ) = [uالأعلى(s){\displaystyle u_{\max }^{(s)}}حيث تشير الأقواس إلى الجزء الكامل من القيمة. يمثل الرقم k ( s ) طول الحقبة، ويُقاس بعدد حقب كاسنر التي تحتويها الحقبة. بالنسبة للحقبة التالية

في سلسلة الأرقام اللانهائية u ، المؤلفة بهذه القواعد ، توجد قيم صغيرة للغاية (ولكنها لا تساوي الصفر أبدًا) x ( s ) وأطوال كبيرة لا نهائية مقابلة k ( s ) .

تصبح سلسلة الحقبة أكثر كثافة عند الاقتراب من t = 0. ومع ذلك، فإن المتغير الطبيعي لوصف المسار الزمني لهذا التطور ليس الوقت العالمي t ، بل لوغاريتمه، ln t ، والذي من خلاله يتم تمديد العملية الكاملة للوصول إلى التفرد إلى −∞.

وفقًا للمعادلة 33 ، فإن إحدى الدوال a و b و c التي تمر بقيمة عظمى خلال الانتقال بين حقب كاسنر، عند ذروة قيمتها العظمى هي

حيث يُفترض أن قيمة a max كبيرة مقارنةً بـ b 0 و c 0 ؛ في المعادلة 38، u هي قيمة المعامل في حقبة كاسنر قبل الانتقال. يتضح من هنا أن قمم القيم العظمى المتتالية خلال كل حقبة تنخفض تدريجيًا. في الواقع، في حقبة كاسنر التالية، تكون قيمة هذا المعامل u' = u − 1، ويتم استبدال Λ وفقًا للمعادلة 36 بالمعادلة Λ' = Λ(1 − 2| p 1 ( u )|). لذلك، فإن نسبة قيمتين عظميين متتاليتين هي

أالأعلىأالأعلى=[ص1(u-1)ص1(u)(1-2|ص1(u)|)]12؛{\displaystyle {\frac {a'_{\max }}{a_{\max }}}=\left[{\frac {p_{1}(u-1)}{p_{1}(u)}}\left(1-2|p_{1}(u)|\right)\right]^{\frac {1}{2}};}

وأخيراً

ما سبق هو حلول لمعادلات أينشتاين في الفراغ. أما بالنسبة لنمط كاسنر النقي، فإن المادة لا تُغير الخصائص النوعية لهذا الحل، ويمكن كتابته ضمنه دون الأخذ في الاعتبار تأثيرها على المجال. مع ذلك، إذا طُبِّق هذا على النموذج قيد البحث، باعتباره حلاً دقيقاً لمعادلات أينشتاين، فإن الصورة الناتجة لتطور المادة لن تكون عامة، بل ستكون خاصة بالتناظر العالي الذي يميز هذا النموذج. رياضياً، ترتبط هذه الخصوصية بحقيقة أنه بالنسبة لهندسة الفضاء المتجانسة التي نناقشها هنا، فإن مركبات موتر ريتشيRα0{\displaystyle R_{\alpha }^{0}}تكون هذه القيم أصفارًا متطابقة، وبالتالي فإن معادلات أينشتاين لا تسمح بحركة المادة (التي تُنتج مركبات غير صفرية لموتر طاقة-زخم الإجهاد).تيα0{\displaystyle T_{\alpha }^{0}}بمعنى آخر ، يجب أن يكون الإطار المتزامن متحركًا أيضًا مع المادة. إذا استبدلنا في المعادلة 19 = 0 و u0 = 1، فإنها تصبح ε ~ ( abc ) −4/3 ~ t −4/3 .

يمكن تجنب هذه الصعوبة إذا اقتصر النموذج على تضمين الحدود الرئيسية للمقياس الحدي (عند t → 0) وكتابة مادة ذات توزيع أولي عشوائي للكثافات والسرعات. عندئذٍ، يتحدد مسار تطور المادة بقوانين حركتها العامة ، المعادلتين 17 و 18، والتي تُنتج المعادلة 21. خلال كل حقبة من حقب كاسنر، تزداد الكثافة وفقًا للقانون.

حيث أن p 3 هو، كما سبق، أكبر الأعداد p 1 و p 2 و p 3. تزداد كثافة المادة بشكل رتيب خلال جميع مراحل التطور نحو التفرد.

تطور المقاييس

تتوافق قيم u الكبيرة جدًا مع قوى كاسنر

وهي قريبة من القيم (0، 0، 1). القيمتان القريبتان من الصفر تكونان متقاربتين أيضًا، وبالتالي فإن التغيرات في نوعين من أنواع "الاضطرابات" الثلاثة (الحدود التي تحتوي على λ وμ وν في الطرف الأيمن من المعادلة 26 ) تكون متشابهة جدًا. إذا كانت هذه الحدود متقاربة جدًا في قيمها المطلقة في بداية حقبة طويلة كهذه، لحظة الانتقال بين حقبتين من حقب كاسنر (أو تم جعلها متقاربة بشكل مصطنع عن طريق تحديد الشروط الابتدائية)، فإنها ستظل متقاربة خلال معظم مدة الحقبة بأكملها. في هذه الحالة (يُطلق عليها BKL حالة التذبذبات الصغيرة )، يصبح التحليل القائم على تأثير نوع واحد من الاضطرابات غير صحيح؛ إذ يجب مراعاة التأثير المتزامن لنوعين من الاضطرابات.

اضطرابان

لنفترض حقبة زمنية طويلة، تشهد خلالها دالتان a و b و c (لنفترض أنهما a و b ) تذبذبات طفيفة، بينما تتناقص الدالة الثالثة ( c ) بشكل رتيب. سرعان ما تصبح الدالة الأخيرة صغيرة؛ لننظر في الحل ضمن المنطقة التي يمكن فيها إهمال c مقارنةً بـ a و b . تُجرى الحسابات أولًا لنموذج الفضاء من النوع التاسع بالتعويض وفقًا لذلك λ = μ = ν = 1. [ 43 ]

بعد إهمال الدالة c ، تعطي المعادلتان الأوليان (المعادلة 26)

ويمكن استخدام المعادلة 28 كمعادلة ثالثة، تأخذ الشكل التالي:

يُكتب حل المعادلة 44 على الصورة التالية:

α+β=2أ02ξ0(τ-τ0)+2lnأ0،{\displaystyle \alpha +\beta ={\frac {2a_{0}^{2}}{\xi _{0}}}\left(\tau -\tau _{0}\right)+2\ln a_{0},}

حيث α₀ و ξ₀ ثابتان موجبان، و τ₀ هي الحد الأعلى للعصر للمتغير τ. ومن الملائم إدخال متغير جديد (بدلاً من τ).

ثم

يتم تحويل المعادلتين 45 و 46 عن طريق إدخال المتغير χ = α − β:

انخفاض قيمة τ من τ₀ إلى −∞ يقابله انخفاض قيمة ξ من ξ₀ إلى 0. وتُحقق الفترة الطويلة ذات القيمتين a و b المتقاربتين (أي ذات قيمة χ صغيرة)، التي نتناولها هنا، إذا كانت ξ₀ قيمة كبيرة جدًا. في الواقع، عند قيم ξ الكبيرة، يكون حل المعادلة 49 في التقريب الأول باستخدام 1/ξ هو

حيث A ثابت؛ المضاعف1ξ{\displaystyle {\tfrac {1}{\sqrt {\xi }}}}يجعل χ كمية صغيرة بحيث يمكن استبدالها في المعادلة 49 بـ sh 2χ ≈ 2χ. [ ملاحظة 11 ]

من المعادلة 50 نحصل على

γξ=14ξ(2χξ2+χ2)=أ2، γ=أ2(ξ-ξ0)+جoنsت.{\displaystyle \gamma _{\xi }={\frac {1}{4}}\xi \left(2\chi _{\xi }^{2}+\chi ^{2}\right)=A^{2},\ \gamma =A^{2}\left(\xi -\xi _{0}\right)+\mathrm {const} .}

بعد تحديد α و β من المعادلة 48 والمعادلة 51 وتوسيع e α و e β في سلسلة وفقًا للتقريب أعلاه، نحصل في النهاية على: [ ملاحظة 12 ]

تُستنتج العلاقة بين المتغير ξ والزمن t من خلال تكامل التعريف dt = abc d τ والذي يعطي

يجب أن تكون قيمة الثابت c₀ ( قيمة c عند ξ = ξ₀ ) الآن c₀{\displaystyle \ll }α 0 ·

فضاء بيانكي من النوع الثامن (المفتوح) يخضع لديناميكيات فوضوية من نوع BKL (Mixmaster) قريبة من التفرد وفقًا للقواعد الواردة في المعادلة 35 مع بدايةu=13{\displaystyle u={\sqrt {13}}}. تقع النقطة الشاذة عند نقطة الانضغاط المركزية لسطح القطع الزائد.

لننظر الآن في المجال ξ{\displaystyle \ll }1. فيما يلي أهم الحدود في حل المعادلة 49 :

χ=α-β=كlnξ+جoنsت،{\displaystyle \chi =\alpha -\beta =k\ln \xi +\mathrm {const} ,\,}

حيث k ثابت في النطاق −1 < k < 1؛ يضمن هذا الشرط أن يكون الحد الأخير في المعادلة 49 صغيرًا (sh 2χ يحتوي على ξ 2 k و ξ −2 k ). بعد ذلك، وبعد تحديد α و β و t ، نحصل على

هذا نمط كاسنر آخر مع وجود الأس السالب t في الدالة c ( t ). [ ملاحظة 13 ]

تُصوّر هذه النتائج تطورًا مشابهًا نوعيًا لما وُصف أعلاه. خلال فترة زمنية طويلة تتوافق مع انخفاض كبير في قيمة ξ، تتذبذب الدالتان a و b ، وتبقى قيمتاهما متقاربتين في المقدار.أ-بأ1ξ{\displaystyle {\tfrac {a-b}{a}}\sim {\tfrac {1}{\sqrt {\xi }}}}في الوقت نفسه، تتحرك كلتا الدالتين أ و ب ببطء (ξ{\displaystyle \sim {\sqrt {\xi }}}تتناقص الدالة. دورة التذبذبات ثابتة بالنسبة للمتغير ξ  : Δξ = 2π (أو، وهو نفسه، بدورة ثابتة بالنسبة للزمن اللوغاريتمي: Δ ln t = 2π Α 2 ). الدالة الثالثة، c ، تتناقص بشكل رتيب وفقًا لقانون قريب من c = c 0 t / t 0 .

يستمر هذا التطور حتى تصبح قيمة ξ ≈ 1، وعندها لا تعود المعادلتان 52 و 53 قابلتين للتطبيق. وتتوافق مدته الزمنية مع تغير قيمة t من t₀ إلى القيمة t₁ ، المرتبطة بـ ξ₀ وفقًا لـ

يمكن تمثيل العلاقة بين ξ و t خلال هذه الفترة بالشكل التالي:

بعد ذلك، وكما هو موضح في المعادلة 55 ، تبدأ الدالة المتناقصة c بالزيادة بينما تبدأ الدالتان a و b بالتناقص. تستمر هذه المرحلة من كاسنر حتى تصبح الحدود / a²b² في المعادلة 22 قريبة من ، وتبدأ سلسلة جديدة من التذبذبات.

يتم الحصول على قانون تغير الكثافة خلال الحقبة الطويلة قيد المناقشة عن طريق استبدال المعادلة 52 في المعادلة 20 :

عندما تتغير قيمة ξ من ξ 0 إلى ξ ≈1، تزداد الكثافةξ02{\displaystyle \xi _{0}^{2}}مرات.

يجب التأكيد على أنه على الرغم من أن الدالة c ( t ) تتغير وفقًا لقانون قريب من c ~ t ، فإن المقياس في المعادلة 52 لا يتوافق مع مقياس كاسنر ذي القوى (0، 0، 1). يتوافق هذا الأخير مع حل دقيق وجده تاوب [ 47 ] وهو مسموح به في المعادلتين 26 و 27 ، وفيه

حيث p و δ1 و δ2 ثوابت . في المنطقة التقاربية τ → −∞، يمكن الحصول من هنا على a = b = ثابت، c = ثابت. t بعد التعويض е рτ = t . في هذا المقياس، تكون النقطة الشاذة عند t = 0 غير فيزيائية.

لنصف الآن الدراسة المماثلة لنموذج النوع الثامن، مع استبدال المعادلات 2628 λ = −1، μ = ν = 1. [ 44 ]

إذا كانت الدالة المتناقصة بشكل رتيب هي a خلال الفترة الطويلة، فلن يتغير شيء في التحليل السابق: بتجاهل في الطرف الأيمن من المعادلتين 26 و 28 ، نعود إلى المعادلتين 49 و 50 نفسيهما (مع تغيير في الترميز). مع ذلك، تحدث بعض التغييرات إذا كانت الدالة المتناقصة بشكل رتيب هي b أو c ؛ لنفترض أنها c .

كما في السابق، لدينا المعادلة 49 بنفس الرموز، وبالتالي، التعبيرات السابقة في المعادلة 52 للدالتين a (ξ) و b (ξ)، ولكن يتم استبدال المعادلة 50 بـ

يصبح الحد الرئيسي عند ξ الكبير الآن

γξ18ξ2،γ18(ξ2-ξ02)،{\displaystyle \gamma _{\xi }\approx {\frac {1}{8}}\xi \cdot 2,\quad \gamma \approx {\frac {1}{8}}\left(\xi ^{2}-\xi _{0}^{2}\right),}

لهذا السبب.

قيمة c كدالة للزمن t هي c = c₀t / t₀ ، لكن اعتماد ξ على الزمن يتغير . يعتمد طول الحقبة الطويلة على ξ₀ وفقًا لـ

من جهة أخرى، تحدد القيمة ξ₀ عدد تذبذبات الدالتين a و b خلال حقبة زمنية (تساوي ξ₀ / 2π). وبالنظر إلى طول الحقبة الزمنية باللوغاريتم الطبيعي (أي بنسبة t₀ / t₁ المعطاة )، فإن عدد التذبذبات للنوع الثامن سيكون، عمومًا، أقل من النوع التاسع. أما بالنسبة لدورة التذبذبات، فنحصل الآن على Δln t = πξ/2؛ وعلى عكس النوع التاسع، فإن الدورة ليست ثابتة طوال الحقبة الزمنية الطويلة، بل تتناقص تدريجيًا مع تناقص قيمة ξ .

المجال الزمني الصغير

تُخالف الحقب الطويلة المسار "المنتظم" للتطور، مما يُصعّب دراسة تطور الفترات الزمنية الممتدة عبر عدة حقب. ومع ذلك، يمكن إثبات أن هذه الحالات "الشاذة" تظهر في التطور التلقائي للنموذج إلى نقطة فريدة في الأزمنة الصغيرة t عند مسافات كبيرة كافية من نقطة بداية ذات شروط ابتدائية عشوائية. حتى في الحقب الطويلة، تظل الدالتان التذبذبيتان خلال الانتقالات بين حقب كاسنر مختلفتين لدرجة أن الانتقال يحدث تحت تأثير اضطراب واحد فقط. جميع النتائج في هذا القسم تنطبق بالتساوي على نماذج النوعين الثامن والتاسع. [ 48 ]

خلال كل حقبة من حقب كاسنر، يكون abc = Λ t ، أي α + β + γ = ln Λ + ln t . عند الانتقال من حقبة (بقيمة معينة للمعامل u ) إلى الحقبة التالية، يُضرب الثابت Λ في 1 + 2 p 1 = (1 – u + u 2 )/(1 + u + u 2 ) < 1. وبالتالي، يحدث انخفاض منتظم في Λ. ولكن من الضروري أن تكون القيمة المتوسطة (بالنسبة لأطوال الحقب k ) للتغير الكلي لـ ln Λ خلال الحقبة محدودة. في الواقع، قد يكون تباعد القيمة المتوسطة ناتجًا فقط عن زيادة سريعة جدًا لهذا التغير مع زيادة k . بالنسبة لقيمة كبيرة للمعامل u ، يكون ln(1 + 2 p 1 ) ≈ −2/ u . بالنسبة لقيمة كبيرة لـ k، تكون القيمة القصوى u (max) = k + x ≈ k. وبالتالي، فإن التغير الكلي لـ ln Λ خلال حقبة زمنية معينة يُعطى بمجموع على الصورة التالية

ln(1+2ص1)=+1ك-2+1ك-1+1ك{\displaystyle \sum \ln \left(1+2p_{1}\right)=\dots +{\frac {1}{k-2}}+{\frac {1}{k-1}}+{\frac {1}{k}}}

مع كتابة الحدود التي تُقابل قيمًا كبيرة لـ u فقط. عندما تزداد قيمة k، يزداد هذا المجموع بمقدار ln k . لكن احتمال ظهور حقبة ذات طول كبير k يتناقص بمقدار 1/ وفقًا للمعادلة 76 ؛ وبالتالي فإن متوسط ​​قيمة المجموع أعلاه محدود. ونتيجة لذلك، فإن التغير المنتظم للكمية ln Λ على مدى عدد كبير من الحقب سيكون متناسبًا مع هذا العدد. ولكن يُلاحظ في المعادلة 85 أنه مع t → 0، يزداد العدد s بمقدار ln |ln t | فقط. وهكذا، في الحد التقاربي لـ t الصغيرة جدًا ، يمكن بالفعل إهمال الحد ln Λ مقارنةً بـ ln t . في هذا التقريب [ ملاحظة 14 ]

حيث تشير Ω إلى "الزمن اللوغاريتمي"

ويمكن اعتبار عملية انتقال الحقب بمثابة سلسلة من ومضات زمنية قصيرة. كما تخضع مقادير القيم القصوى لدوال المقياس المتذبذبة لتغير منتظم. من المعادلة 39، بالنسبة لـ u ≫ 1، يتبين أنأالأعلى-أالأعلى-1/2u{\displaystyle a_{\max }^{\prime }-a_{\max }\approx -1/2u}بنفس الطريقة التي تم بها سابقًا بالنسبة للكمية ln Λ، يمكن استنتاج أن متوسط ​​انخفاض ارتفاع القيم العظمى خلال حقبة زمنية محددة محدود، وأن الانخفاض الكلي على مدى عدد كبير من الحقب الزمنية يزداد مع اقتراب t من الصفر تمامًا مثل ln Ω. في الوقت نفسه، فإن انخفاض القيم الصغرى، وبالتالي زيادة سعة التذبذبات ، يتناسب ( المعادلة 77 ) مع Ω. وفقًا للتقريب المعتمد، يُهمل انخفاض القيم العظمى مقارنةً بزيادة السعات، لذا فإن α max = 0، β max = 0، γ max = 0 للقيم العظمى لجميع الدوال المتذبذبة، وتمر الكميات α، β، γ بقيم سالبة فقط، وترتبط هذه القيم ببعضها البعض في كل لحظة زمنية بالعلاقة الموضحة في المعادلة 63 .

تغيرات α و β و γ كدوال للزمن اللوغاريتمي Ω خلال حقبة واحدة. تشير الخطوط المتقطعة الرأسية إلى تغيرات حقب كاسنر، المقابلة للأجزاء الخطية من المنحنيات. في الأعلى، تظهر قيم المعامل u التي تحدد أسس كاسنر. تكون مدة الحقبة الأخيرة أطول إذا كانت x صغيرة. في الحقبة الأولى من الحقبة التالية، تبدأ γ بالزيادة وتصبح α دالة متناقصة بشكل رتيب.

بالنظر إلى هذا التغير الفوري في الحقب، تُهمل فترات الانتقال لصغرها مقارنةً بطول الحقبة؛ وهذا الشرط مُتحقق بالفعل. [ ملاحظة 15 ] يتطلب استبدال القيم القصوى لـ α و β و γ بأصفار أن تكون الكميات ln (| p 1 |Λ) صغيرة مقارنةً بسعات تذبذبات الدوال المعنية. كما ذُكر سابقًا ، خلال الانتقالات بين الحقب، قد تصبح قيم | p 1 | صغيرة جدًا، بينما لا يرتبط مقدارها واحتمالية حدوثها بسعات التذبذب في اللحظة المعنية. لذلك، من حيث المبدأ، من الممكن الوصول إلى قيم صغيرة جدًا لـ | p 1 | بحيث يُنتهك الشرط المذكور أعلاه (القيم القصوى الصفرية). يمكن أن يؤدي هذا الانخفاض الحاد في α max إلى حالات خاصة مختلفة يصبح فيها الانتقال بين حقب كاسنر وفقًا للقاعدة المذكورة في المعادلة 37 غير صحيح (بما في ذلك الحالات الموصوفة أعلاه ). قد تُخالف هذه الحالات "الخطيرة" القوانين المستخدمة في التحليل الإحصائي أدناه. كما ذكرنا، فإن احتمال حدوث مثل هذه الانحرافات يتقارب بشكل مقارب إلى الصفر؛ وسيتم مناقشة هذه المسألة أدناه.

لنفترض حقبة زمنية تحتوي على k من حقب كاسنر مع وجود مُعامل u يمر عبر القيم

ولنفترض أن α و β هما الدالتان المتذبذبتان خلال هذه الحقبة (الشكل 4). [ ملاحظة 16 ]

اللحظات الابتدائية لحقب كاسنر ذات المعاملات u n هي Ω n . في كل لحظة ابتدائية، تكون إحدى القيمتين α أو β تساوي صفرًا، بينما تكون الأخرى عند قيمة دنيا. القيمتان α أو β في القيم الدنيا المتتالية، أي في اللحظات Ω هما

(دون تمييز بين القيمتين الصغريتين α و β). تتراوح قيم δn التي تقيس هاتين القيمتين الصغريتين بوحدات Ωn المقابلة بين 0 و1. تتناقص الدالة γ بشكل رتيب خلال هذه الفترة؛ ووفقًا للمعادلة 63، فإن قيمتها في اللحظة Ωn هي

خلال الحقبة التي تبدأ عند اللحظة Ω n وتنتهي عند اللحظة Ω n +1 ، تزداد إحدى الدالتين α أو β من −δ n Ω n إلى الصفر بينما تنخفض الأخرى من 0 إلى −δ n +1 Ω n +1 وفقًا لقوانين خطية، على التوالي:

جoنsت+|ص1(uن)|Ω{\displaystyle \mathrm {const} +|p_{1}(u_{n})|\Omega \,}وجoنsت-ص2(uن)Ω{\displaystyle \mathrm {const} -p_{2}(u_{n})\Omega \,}

مما يؤدي إلى علاقة تكرارية

وبالنسبة لطول الحقبة اللوغاريتمية

حيث، باختصار، f ( u ) = 1 + u + . ويُحسب مجموع أطوال n حقبة زمنية باستخدام الصيغة التالية :

يتضح من المعادلة 68 أن |α n+1 | > |α n |، أي أن سعة تذبذب الدالتين α و β تزداد خلال الحقبة بأكملها، حتى وإن كانت العوامل δ n صغيرة. إذا كان الحد الأدنى في بداية الحقبة عميقًا، فلن تصبح الحدود الدنيا التالية أقل عمقًا؛ بعبارة أخرى، يبقى الباقي |α  — β| عند لحظة الانتقال بين حقب كاسنر كبيرًا. لا يعتمد هذا على طول الحقبة لأن الانتقالات بين الحقب تُحدد بالقاعدة العامة في المعادلة 37 حتى بالنسبة للحقب الطويلة.

ترتبط سعة التذبذب الأخير للدالتين α أو β في حقبة زمنية معينة بسعة التذبذب الأول بالعلاقة |α k −1 | = |α 0 | ( k + x ) / (1 + x ). حتى عند قيم k الصغيرة جدًا، يمكن إهمال x مقارنةً بـ بحيث تصبح الزيادة في سعة تذبذب α و β متناسبة مع طول الحقبة الزمنية. بالنسبة للدالتين a = e α و b = e β، هذا يعني أنه إذا كانت سعة تذبذبهما في بداية حقبة زمنية ما هي A 0 ، فستصبح في نهاية هذه الحقبة A 0.أ0ك/(1+x){\displaystyle A_{0}^{k/(1+x)}}.

يزداد طول حقب كاسنر (بالزمن اللوغاريتمي) أيضًا داخل الحقبة الواحدة؛ ومن السهل حساب ذلك من المعادلة 69، حيث أن Δn + 1 > Δn . [ ملاحظة 17 ] إجمالي طول الحقبة هو

(ينشأ الحد الذي يحتوي على 1/ x من الحقبة الأخيرة، الحقبة k ، التي يكون طولها كبيرًا عند قيم x الصغيرة ؛ انظر الشكل 2). اللحظة Ωn التي تنتهي عندها الحقبة k من حقبة معينة هي في نفس اللحظة Ω'0 لبداية الحقبة التالية.

في حقبة كاسنر الأولى من العصر الجديد ، تكون الدالة γ هي أول ما يرتفع من القيمة الدنيا γk = −Ωk (1 − δk ) التي بلغتها في الحقبة السابقة؛ وتلعب هذه القيمة دور سعة البداية δ'0 Ω'0 لسلسلة التذبذبات الجديدة. ومن السهل استنتاج ما يلي:

من الواضح أن δ' 0 Ω' 0 > δ 0 Ω 0. حتى عند قيم k غير الكبيرة جدًا، يكون ازدياد السعة كبيرًا جدًا: تبدأ الدالة c = e γ بالتذبذب من السعةأ0أ0ك2{\displaystyle A_{0}'\sim A_{0}^{k^{2}}}. يتم تأجيل مسألة الحالات "الخطيرة" المذكورة أعلاه لخفض حد التذبذب العلوي بشكل كبير في الوقت الحالي.

وفقًا للمعادلة 40، فإن الزيادة في كثافة المادة خلال الحقب الأولى ( k - 1) تُعطى بالصيغة التالية:

ln(εن+1εن)=2[1-ص3(uن)]Δن+1.{\displaystyle \ln \left({\frac {\varepsilon _{n+1}}{\varepsilon _{n}}}\right)=2\left[1-p_{3}(u_{n})\right]\Delta _{n+1}.}

بالنسبة لآخر k حقبة من حقبة معينة، عند u = x < 1 ، تكون أعلى قيمة للقوة هي p2 ( x ) (وليس p3 ( x )). لذلك، بالنسبة لزيادة الكثافة على مدار الحقبة بأكملها، نحصل على

لذلك، حتى عند قيم k غير الكبيرة جدًا ،ε0/ε0أ02ك{\displaystyle \varepsilon _{0}'/\varepsilon _{0}\sim A_{0}^{2k}}خلال الحقبة التالية (بطول k ' ) ستزداد الكثافة بشكل أسرع بسبب زيادة السعة الابتدائية A 0 ':ε0"/ε0أ02ك"أ02ك2ك{\displaystyle \varepsilon _{0}''/\varepsilon _{0}'\sim A_{0}'^{2k''}\sim A_{0}^{2k^{2}k'}}إلخ. توضح هذه الصيغ الزيادة الحادة في كثافة المادة.

التحليل الإحصائي بالقرب من نقطة التفرد

تكتسب متتالية أطوال الحقب k ( s ) ، المقاسة بعدد حقب كاسنر التي تحتويها، صفة العملية العشوائية بشكل تقاربي. وينطبق الأمر نفسه على متتالية تبادل أزواج الدوال المتذبذبة عند الانتقال من حقبة إلى أخرى (يعتمد ذلك على ما إذا كانت الأعداد k ( s ) زوجية أم فردية). أحد مصادر هذه العشوائية هو القاعدة الموضحة في المعادلتين 41 و 42 ، والتي تنص على أن الانتقال من حقبة إلى أخرى يتحدد في متتالية عددية لانهائية من قيم u . بعبارة أخرى، تنص هذه القاعدة على أنه إذا بدأت المتتالية اللانهائية بأكملها بقيمة ابتدائية معينة، فإن u = 0.uالأعلى(0)=ك(0)+x(0){\displaystyle u_{\max }^{(0)}=k^{(0)}+x^{(0)}}إذاً، فإن أطوال الفترات الزمنية k (0) ، k (1) ، ...، هي الأعداد الموجودة في مفكوك الكسر المستمر البسيط.

يتوافق هذا التوسع مع تحويل الإسقاط للفترة [0, 1] على نفسها بالصيغة Tx = {1/ x }، أي x s +1 = {1/ x s }. ينتمي هذا التحويل إلى ما يُسمى بالتحويلات التوسعية للفترة [0, 1]، أي التحويلات xf ( x ) حيث | f′ ( x )| > 1. تتميز هذه التحويلات بخاصية عدم الاستقرار الأسي: إذا أخذنا نقطتين متقاربتين في البداية، فإن المسافة بينهما تزداد أُسّيًا مع تكرار التحويلات. من المعروف أن عدم الاستقرار الأسي يؤدي إلى ظهور خصائص عشوائية قوية.

من الممكن الانتقال إلى وصف احتمالي لمثل هذا التسلسل من خلال النظر ليس إلى قيمة ابتدائية محددة x (0)، بل إلى القيم x (0) = x الموزعة في الفترة من 0 إلى 1 وفقًا لقانون توزيع احتمالي معين w₀ ( x ). عندئذٍ ، ستكون لقيم x (s) التي تُنهي كل حقبة توزيعات تتبع قوانين معينة wₛ ( x ) . لنفترض أن wₛ (x)dx هو احتمال أن تنتهي الحقبة s بالقيمة xₛ (x) .uالأعلى(s)=x{\displaystyle u_{\max }^{(s)}=x}تقع ضمن نطاق محدد dx .

يمكن أن تنتج القيمة x (s) = x ، التي تُنهي الحقبة s ، من القيم الأولية (لهذه الحقبة).uالأعلى(s)=x+ك{\displaystyle u_{\max }^{(s)}=x+k}، حيث k = 1، 2، ...؛ هذه القيم منuالأعلى(s){\displaystyle u_{\max }^{(s)}}تُقابل هذه القيم x ( s -1) = 1/( k + x ) للفترة السابقة. وبناءً على ذلك، يمكن كتابة العلاقة التكرارية التالية، التي تُعبّر عن توزيع الاحتمالات w ( s ) بدلالة التوزيع w (s -1) ( x ):

ws(x)دx=ك=1ws-1(1ك+x)|د1ك+x|{\displaystyle w_{s}(x)dx=\sum _{k=1}^{\infty }w_{s-1}\left({\frac {1}{k+x}}\right)\left\vert d{\frac {1}{k+x}}\right\vert }

أو

إذا كان التوزيع w s ( x ) يميل مع ازدياد s إلى توزيع حدّي ثابت (مستقل عن s ) w ( x )، فإن هذا التوزيع الأخير يجب أن يحقق معادلة مستنتجة من المعادلة 73c بحذف مؤشرات الدالتين w s −1 ( x ) و w s ( x ). لهذه المعادلة حل

(مُعَيَّر إلى الوحدة ومُقَسَّم إلى الرتبة الأولى من x ). [ ملاحظة 18 ]

لكي يكون طول الحقبة s هو k ، يجب أن تنتهي الحقبة السابقة برقم x يقع بين 1/( k + 1) و 1/ k . وبالتالي، فإن احتمال أن يكون طول الحقبة k يساوي (في حالة الاستقرار)

عند قيم كبيرة لـ k

عند ربط الخصائص الإحصائية للنموذج الكوني بالخصائص الإرجودية للتحويل x s +1 = {1/ x s }، تجدر الإشارة إلى نقطة مهمة. في متتالية لانهائية من الأعداد x المُنشأة وفقًا لهذه القاعدة، ستُلاحظ قيم صغيرة جدًا (لكنها لا تتلاشى أبدًا) لـ x تُقابل أطوالًا كبيرة جدًا k. قد تُؤدي هذه الحالات (ليس بالضرورة!) إلى ظهور مواقف مُحددة يفقد فيها مفهوم العصور، كتسلسلات حقب كاسنر المُتبادلة وفقًا للقاعدة في المعادلة 37 ، معناه (على الرغم من استمرار النمط التذبذبي لتطور النموذج). يُمكن أن يتجلى هذا الوضع "الشاذ"، على سبيل المثال، في ضرورة الاحتفاظ في الجانب الأيمن من المعادلة. 26 مصطلحًا ليس فقط مع إحدى الدوال a و b و c (على سبيل المثال، a 4 )، كما هو الحال في التبادل "العادي" لحقب كاسنر، ولكن في نفس الوقت مع اثنين منها (على سبيل المثال، a 4 و b 4 و a 2 b 2 ).

عند الخروج من سلسلة تذبذبات "شاذة"، تعود دورة من الفترات المنتظمة. يؤكد التحليل الإحصائي لسلوك النموذج، القائم كليًا على التكرارات المنتظمة للتحويلات (المعادلة 42)، نظريةً مهمة: أن احتمال ظهور الحالات الشاذة يؤول تقاربًا إلى الصفر عندما يؤول عدد التكرارات s إلى اللانهاية (أي عندما يؤول الزمن t إلى الصفر)، وهو ما تم إثباته في نهاية هذا القسم. وتعود صحة هذه النظرية إلى حد كبير إلى معدل الزيادة السريع جدًا في سعة التذبذبات خلال كل فترة، وخاصةً عند الانتقال من فترة إلى أخرى.

إن عملية استرخاء النموذج الكوني إلى النظام الإحصائي "الثابت" (مع t → 0 بدءًا من "لحظة أولية" معينة) أقل إثارة للاهتمام، مع ذلك، من خصائص هذا النظام نفسه مع مراعاة القوانين الملموسة لتغير الخصائص الفيزيائية للنموذج خلال العصور المتعاقبة.

يمكن استخلاص فكرة عن معدل استقرار التوزيع من المثال التالي. لنفترض أن القيم الأولية x (0) موزعة في فترة ضيقة عرضها δx ( 0) حول عدد محدد. من العلاقة التكرارية في المعادلة 73c (أو مباشرة من مفكوك المعادلة 73a ) ، من السهل استنتاج أن عرض التوزيعات ws ( x ) (حول أعداد محددة أخرى) سيكون مساويًا لـ

(هذا التعبير صالح فقط طالما أنه يحدد الكميات δ x (s) ≪ 1).

القيمة المتوسطةك¯{\displaystyle {\bar {k}}}، المحسوبة من هذا التوزيع، تتباعد لوغاريتميًا. بالنسبة لمتتالية مقطوعة عند عدد كبير جدًا، ولكنه لا يزال محدودًا N ، يكون لديناك¯lnشمال{\displaystyle {\bar {k}}\sim \ln N}إن فائدة المتوسط ​​في هذه الحالة محدودة للغاية بسبب عدم استقراره: فبسبب التناقص البطيء لـ W ( k )، تتباعد تقلبات k أسرع من متوسطها. ومن الخصائص الأنسب لهذه المتتالية احتمال أن ينتمي عدد مختار عشوائيًا منها إلى حقبة طولها حيث K قيمة كبيرة. هذا الاحتمال هو ln K / ln N. ويكون صغيرًا إذا1كشمال{\displaystyle 1\ll K\ll N}. في هذا الصدد، يمكن القول أن الرقم المختار عشوائياً من التسلسل المعطى ينتمي إلى الحقبة الطويلة باحتمالية عالية.

من الملائم حساب متوسط ​​التعبيرات التي تعتمد في آنٍ واحد على k ( s ) و x ( s ) . بما أن هاتين الكميتين مشتقتان من الكمية نفسها x ( s -1) (التي تُنهي الفترة السابقة)، وفقًا للصيغة k ( s ) + x ( s ) = 1/ x ( s -1) ، فلا يمكن اعتبار توزيعاتهما الإحصائية مستقلة. يمكن الحصول على التوزيع المشترك W s ( k , x ) dx لكلا الكميتين من التوزيع w s -1 ( x ) dx عن طريق استبدال x بـ 1/( x + k ) في الأخير. بعبارة أخرى، تُعطى الدالة W s ( k , x ) بالتعبير الموجود أسفل علامة الجمع في الطرف الأيمن من المعادلة 73c . ​​في حالة الاستقرار، بأخذ w من المعادلة 74 ، نحصل على

إن جمع هذا التوزيع على k يعيدنا إلى المعادلة 74 ، والتكامل بالنسبة إلى dx إلى المعادلة 75 .

تُعاد كتابة الصيغ التكرارية التي تُحدد الانتقالات بين العصور، حيث يُشير الفهرس s إلى العصور المتتالية (وليس حقب كاسنر في عصر مُعين!)، بدءًا من عصر ما ( s = 0) يُعرَّف بأنه عصر ابتدائي. يُمثل Ω ( s ) و ε ( s ) على التوالي اللحظة الابتدائية وكثافة المادة الابتدائية في العصر s ؛ بينما يُمثل δ ( s ) Ω ( s ) سعة التذبذب الابتدائية لزوج الدوال α، β، γ، الذي يتذبذب في العصر المُعين؛ و k ( s ) هو طول العصر s ، و x ( s ) يُحدد طول (عدد حقب كاسنر) العصر التالي وفقًا للمعادلة k ( s + 1) = [1/ x ( s ) ] . وفقًا للمعادلات 71-73

( تم إدخال ξ ( s ) في المعادلة 77 لاستخدامها لاحقًا).

تتمتع الكميات δ ( s ) بتوزيع إحصائي ثابت P (δ) وقيمة متوسطة ثابتة (تقلبات نسبية صغيرة). ولتحديدها، استخدم KL [ 48 ] بالاشتراك مع إيليا ليفشيتز ، شقيق يفغيني ليفشيتز، (مع بعض التحفظات) طريقة تقريبية مبنية على افتراض الاستقلال الإحصائي للكمية العشوائية δ ( s ) والكميتين العشوائيتين k ( s ) و x ( s ) . بالنسبة للدالة P (δ)، تم وضع معادلة تكاملية تعبر عن حقيقة أن الكميتين δ ( s +1) و δ ( s ) المرتبطتين بالعلاقة المذكورة في المعادلة 78 لهما نفس التوزيع؛ وقد تم حل هذه المعادلة عدديًا. في دراسة لاحقة، أظهر خالاتنيكوف وآخرون [ 49 ] أنه يمكن إيجاد التوزيع P (δ) بدقة باستخدام طريقة تحليلية (انظر الشكل 5 ).

بالنسبة للخصائص الإحصائية في حالة الاستقرار، من المنطقي إدخال ما يُسمى بالامتداد الطبيعي للتحويل Tx = {1/ x } عن طريق تمديده بلا حدود إلى المؤشرات السالبة. بعبارة أخرى، هذا انتقال من متتالية لانهائية أحادية الجانب من الأعداد ( x₀ , x₁ , x₂ , ...)، المرتبطة بالمساواة Tx = {1/ x }، إلى متتالية "لانهائية مزدوجة" X = (..., x⁻¹, x₀, x₁, x₂ , ... ) من الأعداد المرتبطة بنفس المساواة لجميع قيم s < –∞ < . بالطبع ، هذا الامتداد ليس فريدًا بالمعنى الحرفي للكلمة (لأن xₛ⁻¹ لا يُحدد بشكل فريد بواسطة xₛ ) ، ولكن جميع الخصائص الإحصائية للمتتالية الممتدة تكون منتظمة على طولها بالكامل، أي أنها ثابتة بالنسبة لأي إزاحة (ويفقد x₀ معناه كشرط "ابتدائي") . المتتالية X مكافئة لمتتالية من الأعداد الصحيحة K = (..., k −1 , k 0 , k 1 , k 2 , ...)، والتي تُنشأ وفقًا للقاعدة k s = [1/ x s –1 ]. وبالعكس، يُحدد كل عدد من X بواسطة الأعداد الصحيحة في K ككسر مستمر لانهائي

(سهولة إدخال الترميز)xs+1+{\displaystyle x_{s+1}^{+}}(مع إزاحة الدليل بمقدار 1، سيتضح ذلك فيما يلي). وللاختصار، يُشار إلى الكسر المستمر ببساطة عن طريق تعداد مقاماته (بين قوسين مربعين)؛ ثم تعريفxs+{\displaystyle x_{s}^{+}}يمكن كتابتها على النحو التالي

تُعرَّف الكميات العكسية بكسر مستمر ذي مقامات رجعية (في اتجاه تناقص الأسس).

يتم تحويل علاقة التكرار في المعادلة 78 بإدخال الرمز η s = (1 − δ s )/δ s مؤقتًا . وبذلك ، يمكن إعادة كتابة المعادلة 78 على النحو التالي:

ηs+1=1ηsxs-1+كs{\displaystyle \eta _{s+1}={\frac {1}{\eta _{s}x_{s-1}+k_{s}}}}

يتم الحصول على كسر متصل لانهائي عن طريق التكرار

ηs+1xs=[كs،كs-1،...]=xs+1-{\displaystyle \eta _{s+1}x_{s}=\left[k_{s},k_{s-1},\dots \right]=x_{s+1}^{-}}

لذلكηs=xs-/xs+{\displaystyle \eta _{s}=x_{s}^{-}/x_{s}^{+}}وأخيراً

يحتوي هذا التعبير عن δ s على كميتين عشوائيتين فقط (بدلاً من الثلاث في [ 48 ] ).xs+{\displaystyle x_{s}^{+}}وxs-{\displaystyle x_{s}^{-}}، كل منها يأخذ قيمًا في الفترة [0، 1].

ويترتب على تعريف المعادلة 79c ما يلي:1/xs-=xs-+كs=xs-+[1/xs+]{\displaystyle 1/x_{s}^{-}=x_{s}^{-}+k_{s}=x_{s}^{-}+\left[1/x_{s}^{+}\right]}وبالتالي، فإن إزاحة المتتالية X بأكملها خطوة واحدة إلى اليمين تعني تحويلاً مشتركاً للكميات.xs+{\displaystyle x_{s}^{+}}وxs-{\displaystyle x_{s}^{-}}وفق

هذا تحويل أحادي في المربع الواحدي . وبالتالي، لدينا الآن تحويل أحادي لكميتين بدلاً من تحويل غير أحادي Tx = {1/ x } لكمية واحدة.

الكمياتxs+{\displaystyle x_{s}^{+}}وxs-{\displaystyle x_{s}^{-}}لدينا توزيع مشترك ثابت P ( x + , x- ). بما أن المعادلة 79e هي تحويل أحادي، فإن شرط ثبات التوزيع يُعبَّر عنه ببساطة بمعادلة دالة .

حيث J هي جاكوبيان التحويل.

يؤدي تحريك المتتالية X بمقدار خطوة واحدة إلى التحويل التالي T للمربع الواحد:

x=1x،y=11x+y{\displaystyle x^{\prime }={\frac {1}{x}},\quad y^{\prime }={\frac {1}{{\frac {1}{x}}+y}}}

(معxx0+{\displaystyle x\equiv x_{0}^{+}}،yx0-{\displaystyle y\equiv x_{0}^{-}}انظر المعادلة 79هـ ). تُعرّف دالة الكثافة P ( x , y ) المقياس الثابت لهذا التحويل. من الطبيعي افتراض أن P ( x , y ) دالة متناظرة لـ x و y . هذا يعني أن المقياس ثابت بالنسبة للتحويل S ( x , y ) = ( y , x )، وبالتالي بالنسبة لحاصل الضرب ST مع ST ( x , y ) = ( x , y ) و

x"=11x+y،y"=1x{\displaystyle x''={\frac {1}{{\frac {1}{x}}+y}},\quad y''={\frac {1}{x}}}

من الواضح أن ST لها تكامل أول H = 1/( x + y) . على الخط H = ثابت c ، يكون للتحويل الشكل التالي:

1x"=[1x]+y=[1x]+ج-1x=ج-{1x}{\displaystyle {\frac {1}{x''}}=\left[{\frac {1}{x}}\right]+y=\left[{\frac {1}{x}}\right]+c-{\frac {1}{x}}=c-\left\{{\frac {1}{x}}\right\}}

وبالتالي، يجب أن تكون كثافة القياس الثابت لـ ST على الشكل التالي

و(ج) دج د1x=و(1x+y)1x2دx دy{\displaystyle f(c)\ dc\ d{\frac {1}{x}}=f\left({\frac {1}{x}}+y\right){\frac {1}{x^{2}}}dx\ dy}

وبأخذ التناظر P ( x , y )= P ( y , x ) في الاعتبار، يصبح هذا f ( c )= c −2 وبالتالي (بعد التطبيع)

(يؤدي تكاملها على x + أو x إلى الدالة w ( x ) المعادلة 74 ). وقد استخدم تشيرنوف وبارو [ 50 ] اختزال التحويل إلى دالة أحادية، وحصلا على صيغة من شكل المعادلة 79g ولكن لمتغيرات أخرى؛ ولا تتضمن ورقتهما تطبيقات على المشكلات التي تناولها خالاتنيكوف وآخرون [ 49 ].

يمكن التحقق من صحة المعادلة 79g أيضًا عن طريق حساب مباشر؛ مصفوفة جاكوبي للتحويل في المعادلة 79e هي

ج=(xs+1+،xs+1-)(xs+،xs-)=xs+1+xs+xs+1-xs-=(xs+1+xs+)2{\displaystyle J={\frac {\partial \left(x_{s+1}^{+},x_{s+1}^{-}\right)}{\partial \left(x_{s}^{+},x_{s}^{-}\right)}}={\frac {\partial x_{s+1}^{+}}{\partial x_{s}^{+}}}{\frac {\partial x_{s+1}^{-}}{\partial x_{s}^{-}}}=\left({\frac {x_{s+1}^{+}}{x_{s}^{+}}}\right)^{2}}

(يجب ملاحظة ذلك في حساباته)[1/xs+]+{1/xs+}=1/xs+{\displaystyle \left[1/x_{s}^{+}\right]+\left\{1/x_{s}^{+}\right\}=1/x_{s}^{+}}).

دالة توزيع الاحتمال P (δ). الخط الأحمر: الدالة الدقيقة (المعادلة 79h ). الخط الأزرق: الحل التقريبي للمعادلة التكاملية [ 48 ] . يبدو المنحنيان متشابهين بشكل ملحوظ، ومتوسط ​​كلا التوزيعين يساوي 0.50. [ ملاحظة 19 ]

بما أن δ s مُعبر عنها في المعادلة 79d بدلالة الكميتين العشوائيتين x + و x ، فإن معرفة توزيعهما المشترك تُتيح حساب التوزيع الإحصائي P (δ) بتكامل P ( x + , x ) على أحد المتغيرين عند قيمة ثابتة لـ δ. ونظرًا لتناظر الدالة في المعادلة 79g بالنسبة للمتغيرين x + و x ، فإن P (δ) = P (1 − δ)، أي أن الدالة P (δ) متناظرة حول النقطة δ = 1/2. إذن

P(دلتا) ددلتا=ددلتا01P(x+،x+دلتا1-دلتا)(x-دلتا)x+دx+{\displaystyle P(\delta )\ d\delta =d\delta \int _{0}^{1}P\left(x^{+},{\frac {x^{+}\delta }{1-\delta }}\right)\left({\frac {\partial x^{-}}{\partial \delta }}\right)_{x^{+}}dx^{+}}

عند تقييم هذا التكامل (لـ 0 ≤ δ ≤ 1/2 ثم باستخدام التناظر المذكور أعلاه)، أخيرًا

القيمة المتوسطةدلتا¯{\displaystyle {\bar {\delta }}}= 1/2 بالفعل نتيجة لتناظر الدالة P (δ). وبالتالي، فإن متوسط ​​قيمة السعة الأولية (في كل حقبة) لتذبذبات الدوال α، β، γ يزداد بمقدار Ω/2.

يتم إيجاد العلاقة الإحصائية بين الفترات الزمنية الكبيرة Ω وعدد الحقب s التي تحتويها من خلال التطبيق المتكرر للمعادلة 77 :

مع ذلك، فإن حساب المتوسط ​​المباشر لهذه المعادلة غير منطقي: نظرًا للتناقص البطيء للدالة W ( k ) في المعادلة 76 ، فإن القيم المتوسطة للكمية expξ ( s ) غير مستقرة بالمعنى المذكور أعلاه - إذ تزداد التقلبات بسرعة أكبر من القيمة المتوسطة نفسها مع ازدياد نطاق حساب المتوسط. ويمكن التخلص من هذا التذبذب بأخذ اللوغاريتم: الفترة الزمنية "اللوغاريتمية المزدوجة".

يُعبَّر عنه بمجموع الكميات ξ ( p ) التي لها توزيع إحصائي مستقر. القيمة المتوسطة لـ τ هيτ¯=sξ¯{\displaystyle {\bar {\tau }}=s{\bar {\xi }}}لحسابξ¯{\displaystyle {\bar {\xi }}}لاحظ أنه يمكن إعادة كتابة المعادلة 77 على النحو التالي

بالنسبة للتوزيع الثابتlnxs¯=lnxs-1¯{\displaystyle {\overline {\ln x_{s}}}={\overline {\ln x_{s-1}}}}وبفضل تناظر الدالة P (δ) أيضًاlnدلتاs¯=ln(دلتاs+1)¯{\displaystyle {\overline {\ln \delta _{s}}}={\overline {\ln \left(\delta _{s+1}\right)}}}. لذلك

ξ¯=-2lnx¯=-201w(x)lnx دx=π26ln2=2.37{\displaystyle {\bar {\xi }}=-2{\overline {\ln x}}=-2\int _{0}^{1}w(x)\ln x\ dx={\frac {\pi ^{2}}{6\ln 2}}=2.37}

( w ( x ) من المعادلة 74 ). وبالتالي

والذي يحدد متوسط ​​الفترة الزمنية اللوغاريتمية المزدوجة التي تحتوي على s حقب متتالية.

بالنسبة لقيم s الكبيرة، يكون عدد الحدود في مجموع المعادلة 81 كبيرًا، ووفقًا للنظريات العامة لنظرية الإرجودية، تتوزع قيم τ حولτs¯{\displaystyle {\overline {\tau _{s}}}}وفقًا لقانون جاوس مع الكثافة

يُعد حساب التباين D τ أكثر تعقيدًا، لأنه لا يتطلب معرفة فقطξ¯{\displaystyle {\bar {\xi }}}وξ2¯{\displaystyle {\overline {\xi ^{2}}}}هناك حاجة إلى ذلك، ولكن أيضًا إلى الارتباطاتξصξص¯{\displaystyle {\overline {\xi _{p}\xi _{p\prime }}}}يمكن تبسيط الحساب بإعادة ترتيب الحدود في معادلة المجموع 81. باستخدام المعادلة 81أ، يمكن إعادة كتابة المجموع على النحو التالي:

ص=1sξص=lnص=1sدلتاص(1-دلتاص+1)xصxص-1=lnص=1sدلتاص(1-دلتاص)xص-12+lnx0xs+ln1-دلتا11-دلتاs+1{\displaystyle \sum _{p=1}^{s}\xi _{p}=\ln \prod _{p=1}^{s}{\frac {\delta _{p}}{\left(1-\delta _{p+1}\right)x_{p}x_{p-1}}}=\ln \prod _{p=1}^{s}{\frac {\delta _{p}}{\left(1-\delta _{p}\right)x_{p-1}^{2}}}+\ln {\frac {x_{0}}{x_{s}}}+\ln {\frac {1-\delta _{1}}{1-\delta _{s+1}}}}

لا يزداد الحدّان الأخيران مع ازدياد قيمة s ؛ ويمكن حذف هذين الحدّين لأن قوانين الحدّ للقيم الكبيرة s هي السائدة.

(يؤخذ التعبير 79d لـ δ p في الاعتبار). وبنفس الدقة (أي حتى الحدود التي لا تزداد مع s ) المساواة

صحيح. في الواقع، بموجب المعادلة 79هـ

xص+1++1xص+1-=1xص++xص-{\displaystyle x_{p+1}^{+}+{\frac {1}{x_{p+1}^{-}}}={\frac {1}{x_{p}^{+}}}+x_{p}^{-}}

وبالتالي

ln(1+xص+1+xص+1-)-lnxص+1-=ln(1+xص+xص-)-lnxص+{\displaystyle \ln \left(1+x_{p+1}^{+}x_{p+1}^{-}\right)-\ln x_{p+1}^{-}=\ln \left(1+x_{p}^{+}x_{p}^{-}\right)-\ln x_{p}^{+}}

بجمع هذه المتطابقة على p ، نحصل على المعادلة 82c . وأخيرًا، بنفس الدقةxص+{\displaystyle x_{p}^{+}}يتم تغيير x p تحت علامة الجمع، وبالتالي يتم تمثيل τ s على النحو التالي

تباين هذا المجموع في حالة القيم الكبيرة لـ s هو

يُؤخذ في الاعتبار أنه بفضل التجانس الإحصائي للتسلسل فإن الارتباطاتηصηص¯{\displaystyle {\overline {\eta _{p}\eta _{p\prime }}}}يعتمد فقط على الفروق | pp |. القيمة المتوسطةη¯=ξ¯{\displaystyle {\bar {\eta }}={\bar {\xi }}}متوسط ​​المربع

η2¯=401w(x)ln2x دx=6ξ(3)ln2=10.40{\displaystyle {\overline {\eta ^{2}}}=4\int _{0}^{1}w(x)\ln ^{2}x\ dx={\frac {6\xi (3)}{\ln 2}}=10.40}

مع الأخذ في الاعتبار أيضًا قيم الارتباطاتη0ηص¯{\displaystyle {\overline {\eta _{0}\eta _{p}}}}مع p = 1، 2، 3 (محسوبة عدديًا) يتم الحصول على النتيجة النهائية D τ s = (3.5 ± 0.1) s .

مع زيادة قيمة يزداد التذبذب النسبيدτs/τs¯{\displaystyle D_{{\tau }_{s}}/{\overline {\tau _{s}}}}يؤول إلى الصفر عندما s −1/2 . بعبارة أخرى، تصبح العلاقة الإحصائية في المعادلة 82 شبه مؤكدة عند قيم s الكبيرة . وهذا يُتيح عكس العلاقة، أي تمثيلها على أنها اعتماد متوسط ​​عدد الفترات s τ التي يتم تبادلها في فترة زمنية معينة τ من الزمن اللوغاريتمي المزدوج.

إن التوزيع الإحصائي للقيم الدقيقة لـ s τ حول متوسطها هو أيضًا توزيع غاوسي بتباين σ

دsτ=3.5sτ¯3τ2=0.26τ{\displaystyle D_{s_{\tau }}=3.5{\frac {{\overline {s_{\tau }}}^{3}}{\tau ^{2}}}=0.26\tau }

يتم إعطاء التوزيع الإحصائي المعني بواسطة نفس التوزيع الغاوسي الذي يكون فيه المتغير العشوائي الآن s τ عند قيمة τ معينة:

من وجهة النظر هذه، فإن مصدر السلوك الإحصائي هو التعسف في اختيار نقطة البداية للفترة τ المتراكبة على التسلسل اللانهائي للحقب المتبادلة.

بالنسبة لكثافة المادة، يمكن إعادة كتابة المعادلة 79 مع مراعاة المعادلة 80 بالشكل التالي:

lnlnε(s+1)ε(s)=ηs+ص=0s-1ξص،ηs=ln[2دلتا(s)(ك(s)+x(s)-1)Ω(0)]{\displaystyle \ln \ln {\frac {\varepsilon ^{(s+1)}}{\varepsilon ^{(s)}}}=\eta _{s}+\sum _{p=0}^{s-1}\xi _{p},\quad \eta _{s}=\ln \left[2\delta ^{(s)}\left(k^{(s)}+x^{(s)}-1\right)\Omega ^{(0)}\right]}

ثم، بالنسبة لتغير الطاقة الكلي خلال حقبة زمنية معينة ،

يُشكّل الحدّ الذي يحتوي على المجموع في p المساهمة الرئيسية في هذا التعبير لأنه يتضمن أسًا ذا قوة كبيرة. وبترك هذا الحدّ فقط وأخذ متوسط ​​المعادلة 87 ، نحصل في طرفها الأيمن على التعبير التاليsξ¯{\displaystyle s{\bar {\xi }}}وهذا يتطابق مع المعادلة 82 ؛ جميع الحدود الأخرى في المجموع (بما في ذلك الحدود التي تحتوي على η s في أسسها) تؤدي فقط إلى تصحيحات من رتبة نسبية 1/ s . لذلك،

بفضل الطبيعة شبه المؤكدة للعلاقة بين τ s و s ، يمكن كتابة المعادلة 88 على النحو التالي

lnln(ετ/ε(0))¯=τأوlnln(ε(s)/ε(0))¯=2.1s،{\displaystyle {\overline {\ln \ln \left(\varepsilon _{\tau }/\varepsilon ^{(0)}\right)}}=\tau \quad {\text{or}}\quad {\overline {\ln \ln \left(\varepsilon ^{(s)}/\varepsilon ^{(0)}\right)}}=2.1s,}

والذي يحدد قيمة اللوغاريتم المزدوج لزيادة الكثافة المتوسطة بواسطة فترات زمنية لوغاريتمية مزدوجة معينة τ أو بواسطة عدد معين من الحقب s .

توجد هذه العلاقات الإحصائية المستقرة تحديدًا لفترات زمنية لوغاريتمية مزدوجة ولزيادة الكثافة. أما بالنسبة لخصائص أخرى، مثل ln (ε ( s )(0) ) أو Ω (s) / Ω (0) = exp τ فإن التذبذب النسبي يزداد أُسّيًا مع زيادة نطاق المتوسط، مما يُفقد مصطلح القيمة المتوسطة معناه المستقر.

يمكن تتبع أصل العلاقة الإحصائية في المعادلة 88 من القانون الأولي الذي يحكم تغير الكثافة خلال حقب كاسنر الفردية. وفقًا للمعادلة 21 ، خلال التطور بأكمله لدينا

lnlnε(ت)=ثابت+lnΩ+ln2(1-ص3(ت))،{\displaystyle \ln \ln \varepsilon (t)={\text{const}}+\ln \Omega +\ln 2(1-p_{3}(t)),}

مع تغير قيمة 1 − p 3 ( t ) من حقبة إلى أخرى، مرورًا بقيم في الفترة من 0 إلى 1. يزداد الحد ln Ω = ln ln (1/ t ) بشكل مطرد؛ من ناحية أخرى، لا يمكن للحد ln2(1 − p 3 ) أن يأخذ قيمًا كبيرة (مقارنة بـ ln Ω) إلا عندما تظهر قيم p 3 قريبة جدًا من الواحد (أي، قيمة صغيرة جدًا لـ | p 1 |). هذه هي تحديدًا الحالات "الخطيرة" التي تعيق المسار المنتظم للتطور المعبر عنه بالعلاقات المتكررة في المعادلات من 77 إلى 79 .

يبقى أن نُثبت أن مثل هذه الحالات لا تظهر فعليًا في نظام الحد التقاربي. يبدأ التطور التلقائي للنموذج عند لحظة معينة تُحدد فيها شروط ابتدائية محددة بطريقة اعتباطية. وعليه، يُقصد بـ"التقاربي" نظامًا بعيدًا بما فيه الكفاية عن اللحظة الابتدائية المختارة.

تُعدّ الحالات الخطيرة تلك التي تظهر فيها قيم صغيرة للغاية للمعامل u = x (وبالتالي | p 1 | ≈ x ) في نهاية حقبة زمنية. ويُعدّ معيار اختيار هذه الحالات هو المتباينة

حيث | α ( s ) | هو عمق الحد الأدنى الأولي للوظائف التي تتذبذب في العصر s (سيكون من الأنسب اختيار السعة النهائية، لكن ذلك من شأنه أن يعزز معيار الاختيار فقط).

تُحدد الشروط الابتدائية قيمة x (0) في الحقبة الأولى. وتُعتبر القيم في الفترة δx (0) ~ exp(−|α (0) |) خطيرة، وكذلك في الفترات التي قد تُؤدي إلى حالات خطيرة في الحقب اللاحقة. ولكي تقع x ( s ) في الفترة الخطيرة δx ( s ) ~ exp(−|α ( s ) |)، يجب أن تقع القيمة الابتدائية x (0) في فترة عرضها δx ( 0) ~ δx ( s ) / k ( 1)^2 ... k ( s )^2 . [ 51 ] لذلك، من فترة الوحدة التي تضم جميع القيم الممكنة لـ x (0) ، ستظهر الحالات الخطيرة في الأجزاء λ من هذه الفترة.

(يُؤخذ المجموع الداخلي على جميع القيم k (1) ، k (2) ، ...، k ( s ) من 1 إلى ∞). من السهل إثبات أن هذه الحقبة تتقارب إلى القيمة λ{\displaystyle \ll }1، الذي يُحدد مقداره بالحد الأول في المعادلة 90. ويمكن إثبات ذلك من خلال هيمنة قوية للحقبة التي يُستبدل بها |α ( s ) | = (s + 1) |α (0) |، بغض النظر عن أطوال الحقب k (1) ، k (2) ، ... (في الواقع، تزداد قيمة |α ( s ) | بشكل أسرع بكثير؛ حتى في أسوأ الحالات k (1) = k (2) = ... = 1، تزداد قيم |α ( s ) | مع q s(0) | حيث q > 1). مع ملاحظة أن

ك1ك(1)2ك(2)2...ك(s)2=(π2/6)s{\displaystyle \sum _{k}{\frac {1}{k^{(1)^{2}}k^{(2)^{2}}\dots k^{(s)^{2}}}}=\left(\pi ^{2}/6\right)^{s}}

يحصل المرء

λ=خبرة(|-α(0)|)s=0[(π2/6)خبرة(|-α(0)|)]sخبرة(|-α(0)|).{\displaystyle \lambda =\exp \left(\left|-\alpha ^{(0)}\right|\right)\sum _{s=0}^{\infty }\left[\left(\pi ^{2}/6\right)\exp \left(\left|-\alpha ^{(0)}\right|\right)\right]^{s}\approx \exp \left(\left|-\alpha ^{(0)}\right|\right).}

إذا كانت القيمة الابتدائية لـ x (0) تقع خارج المنطقة الخطرة λ، فلن تحدث أي حالات خطرة. أما إذا كانت تقع داخل هذه المنطقة، فستحدث حالات خطرة، ولكن بعد انتهائها، يستأنف النموذج تطوره "المنتظم" بقيمة ابتدائية جديدة قد تدخل المنطقة الخطرة من حين لآخر (باحتمالية λ). وتتكرر الحالات الخطرة باحتمالات λ₂ ، λ₃ ، ...، تتقارب تقاربًا مقاربًا إلى الصفر.

حل عام مع تذبذبات صغيرة

في النماذج المذكورة أعلاه، دُرست تطورات المقياس بالقرب من نقطة التفرد باستخدام مثال مقاييس الفضاء المتجانسة. ويتضح من خصائص هذا التطور أن البناء التحليلي للحل العام لنقطة تفرد من هذا النوع يجب أن يُجرى بشكل منفصل لكل مكون من مكونات التطور الأساسية: بالنسبة لحقب كاسنر، ولعملية الانتقالات بين الحقب الناتجة عن "الاضطرابات"، وللفترات الطويلة التي يحدث فيها اضطرابان في آن واحد. خلال حقبة كاسنر (أي عند الاضطرابات الصغيرة)، يُعطى المقياس بالمعادلة 7 دون الشرط λ = 0.

قام BKL بتطوير نموذج مستقل عن توزيع المادة (متجانس أو غير متجانس) لفترات طويلة مع تذبذبات صغيرة. وقد تبين أن التغير الزمني لهذا الحل مشابه جدًا للتغير في حالة النماذج المتجانسة؛ ويمكن الحصول على هذا الأخير من النموذج المستقل عن التوزيع عن طريق اختيار خاص للدوال العشوائية التي يحتويها. [ 52 ]

مع ذلك، من الملائم بناء الحل العام في نظام إحداثيات يختلف نوعًا ما عن الإطار المرجعي المتزامن: g₀α = 0 كما في الإطار المتزامن، ولكن بدلًا من g₀₀ = 1 ، أصبح g₀₀ = −g₃₃ . بتعريف موتر قياس الفضاء γαβ = −gαβ مرة أخرى ، نحصل بالتالي على

يُكتب إحداثي الفضاء الخاص على الصورة x₃ = ويُكتب إحداثي الزمن على الصورة x₀ = ξ (باعتباره مختلفًا عن الزمن الخاص t ). سيُبين أن ξ يُقابل المتغير نفسه المُعرَّف في النماذج المتجانسة. يُشار إلى التفاضل بالنسبة إلى ξ و z ، على التوالي، بالنقطة والشرطة. تأخذ المؤشرات اللاتينية a و b و c القيم 1 و 2، والتي تُقابل إحداثيات الفضاء x₁ و x₂ ، والتي ستُكتب أيضًا على الصورة x و y . بالتالي، فإن المقياس هو

يجب أن يحقق الحل المطلوب المتباينات

(تحدد هذه الشروط أن إحدى الدوال a 2 ، b 2 ، c 2 صغيرة مقارنة بالدالتين الأخريين، وهو ما كان عليه الحال أيضًا مع النماذج المتجانسة).

تعني المتباينة 94 أن المكونات γ a 3 صغيرة بمعنى أنه عند أي نسبة بين الإزاحتين dx a و dz ، يمكن حذف الحدود التي تحتوي على حاصل ضرب dx a dz في مربع عنصر الطول المكاني dl 2. لذلك، فإن التقريب الأول للحل هو مقياس معادلة 92 حيث γ a 3 = 0: [ ملاحظة 20 ]

يمكن للمرء أن يقتنع بسهولة عن طريق حساب مكونات موتر ريتشيR00{\displaystyle R_{0}^{0}}،R30{\displaystyle R_{3}^{0}}،R33{\displaystyle R_{3}^{3}}،Rأب{\displaystyle R_{a}^{b}}باستخدام المعادلة المترية 95 والشرط الوارد في المعادلة 93، والذي ينص على أن جميع الحدود التي تحتوي على مشتقات بالنسبة للإحداثيات x a صغيرة مقارنةً بالحدود التي تحتوي على مشتقات بالنسبة للإحداثيات ξ و z (نسبتها ~ γ 33 / γ ab ). بعبارة أخرى، للحصول على معادلات التقريب الرئيسي، يجب اشتقاق γ 33 و γ ab في المعادلة 95 كما لو أنهما لا يعتمدان على x a .

نحصل على المعادلات التالية: [ ملاحظة 21 ]

يتم رفع وخفض المؤشر هنا بمساعدة γ ab . الكمياتϰ{\displaystyle \varkappa }وλ هما الانقباضاتϰأأ{\displaystyle \varkappa _{a}^{a}}وλأأ{\displaystyle \lambda _{a}^{a}}وبذلك

أما بالنسبة لمكونات موتر ريتشيRأ0{\displaystyle R_{a}^{0}}،Rأ3{\displaystyle R_{a}^{3}}وبناءً على هذه الحسابات، فإن قيمتيهما تساويان صفرًا. وفي التقريب التالي (أي مع مراعاة صغر قيمة γ a 3 والمشتقات بالنسبة إلى x و y )، يتم تحديد الكميتين γ a 3 من خلال القيمتين γ 33 و γ ab المعروفتين مسبقًا .

يؤدي اختصار المعادلة 97 إلىجي+جي¨=0{\displaystyle G^{\prime \prime }+{\ddot {G}}=0}وبالتالي،

تختلف الحالات الممكنة باختلاف المتغير G. في الحالة المذكورة أعلاه ، g 00 = γ 33{\displaystyle \gg }γ ab وشمالز٠٠(جي˙)2-γ33(جي)2=4γ33و˙1و˙2{\displaystyle N\approx g^{00}\left({\dot {G}}\right)^{2}-\gamma ^{33}\left(G^{\prime }\right)^{2}=4\gamma ^{33}{\dot {f}}_{1}{\dot {f}}_{2}}الحالة N > 0 (الكمية N متغيرة مع الزمن) تؤدي إلى حالات تفرد زمنية مهمة. بالتعويض في المعادلة 101، f₁ = ½ (ξ + z ) sin y و f₂ = ½ (ξ − z ) sin y ، نحصل على G من النوع

لا يُقلل هذا الاختيار من عمومية الاستنتاجات؛ إذ يُمكن إثبات إمكانية العمومية (في التقريب الأول) بمجرد النظر إلى التحويلات المسموح بها المتبقية للمتغيرات. عند N < 0 (حيث N كمية مكانية)، يُمكن استبدال G = z، مما يُعمم مقياس أينشتاين-روزن المعروف . [ 53 ] عند N = 0، نصل إلى مقياس موجة روبنسون-بوندي [ 54 ] الذي يعتمد فقط على ξ + z أو فقط على ξ − z (انظر [ 55 ] ). تم وضع العامل sin y في المعادلة 102 للمقارنة المريحة مع النماذج المتجانسة. مع الأخذ في الاعتبار المعادلة 102 ، تُصبح المعادلات من 97 إلى 99 كما يلي:

المعادلات الرئيسية هي المعادلة 103 التي تحدد مكونات γ ab ؛ ثم يتم إيجاد الدالة ψ من خلال تكامل بسيط للمعادلة 104 - المعادلة 105 .

تتراوح قيم المتغير ξ من 0 إلى ∞. يُنظر في حل المعادلة 103 عند حدين، ξ{\displaystyle \gg }1 و{\displaystyle \ll }1. عند قيم ξ الكبيرة، يمكن البحث عن حل يأخذ شكل تحليل 1 / ξ :

وبذلك

( تتطلب المعادلة 107 تحقق الشرط 102 ). بالتعويض بالمعادلة 103 في المعادلة 106 ، نحصل في الرتبة الأولى على

حيث تشكل الكميات a و ac مصفوفة معكوسة للمصفوفة a و ac . حل المعادلة 108 له الشكل التالي :

حيث أن l a و m a و ρ هي دوال اختيارية للإحداثيات x و y مقيدة بالشرط eq. 110 المشتق من eq. 107 .

لإيجاد الحدود العليا لهذا التحليل، من الملائم كتابة مصفوفة الكميات المطلوبة γ ab بالشكل التالي:

حيث يرمز الرمز ~ إلى منقولة المصفوفة. المصفوفة H متناظرة وأثرها يساوي صفرًا. تضمن المعادلة 111 تناظر γ ab وتحقيق الشرط الوارد في المعادلة 102. إذا استُبدل exp H بالعدد 1، نحصل من المعادلة 111 على γ ab = ξ a ab ، حيث a ab مأخوذة من المعادلة 109. بعبارة أخرى، يتوافق الحد الأول من تحليل γ ab مع H = 0؛ أما الحدود الأعلى فتُحسب بتحليل القوى للمصفوفة H التي تُعتبر مكوناتها صغيرة.

تُكتب المكونات المستقلة للمصفوفة H على النحو σ و φ بحيث

بإدخال المعادلة 111 في المعادلة 103 مع الإبقاء فقط على الحدود الخطية بواسطة H ، نحصل على σ و φ

σ¨+ξ-1σ˙-σ=0،{\displaystyle {\ddot {\sigma }}+\xi ^{-1}{\dot {\sigma }}-\sigma ^{\prime \prime }=0,}

إذا حاولنا إيجاد حل لهذه المعادلات على شكل متسلسلة فورييه بدلالة الإحداثي z ، فسنحصل على معادلات بيسل كدوال لمعاملات المتسلسلة في ξ. أما الحدود التقاربية الرئيسية للحل عند قيم ξ الكبيرة فهي [ ملاحظة 22 ].

σ=1ξن=-(أ1نهـأنانωξ+ب1نهـ-أنانωξ)هـأنانωz،{\displaystyle \sigma ={\frac {1}{\sqrt {\xi }}}\sum _{n=-\infty }^{\infty }\left(A_{1n}e^{in\omega \xi }+B_{1n}e^{-in\omega \xi }\right)e^{in\omega z},}
ωن2=ن2ω2+4ρ2.{\displaystyle \omega _{n}^{2}=n^{2}\omega ^{2}+4\rho ^{2}.}

المعاملان A و B دالتان عقديتان اختياريتان للإحداثيات x و y ، وتحققان الشروط اللازمة لقيمتي σ و φ الحقيقيتين؛ التردد الأساسي ω دالة حقيقية اختيارية للإحداثيات x و y . الآن، من المعادلتين 104 و 105، من السهل الحصول على الحد الأول من الدالة ψ:

(هذا الحد يختفي إذا كانت ρ = 0؛ في هذه الحالة يكون الحد الرئيسي هو الحد الخطي لـ ξ من التفكيك: ψ = ξ q ( x , y ) حيث q دالة موجبة [ 56 ] ).

لذلك، عند قيم ξ الكبيرة، تتذبذب مركبات موتر القياس γ ab مع انخفاض ξ على خلفية انخفاض بطيء ناتج عن انخفاض عامل ξ في المعادلة 111. ينخفض ​​المركب γ 33 = e ψ بسرعة وفقًا لقانون قريب من exp (ρ 2 ξ 2 )؛ وهذا ما يجعل الشرط في المعادلة 93 ممكنًا . [ ملاحظة 23 ]

بعد ذلك، ستنظر BKL في الحالة ξ{\displaystyle \ll }1. يتم إيجاد التقريب الأول لحل المعادلة 103 من خلال الافتراض (الذي أكدته النتيجة) أنه في هذه المعادلات يمكن حذف الحدود التي تحتوي على مشتقات بالنسبة للإحداثيات:

هذه المعادلة بالإضافة إلى الشرط الوارد في المعادلة 102 تعطي

حيث أن λ a و μ a و s 1 و s 2 هي دوال اختيارية لجميع الإحداثيات الثلاثة x و y و z ، والتي ترتبط بشروط أخرى

تعطينا المعادلتان 104 و 105 الآن

المشتقاتλأب{\displaystyle {\lambda _{a}^{b}}^{\prime }}تحتوي القيم المحسوبة باستخدام المعادلة 118 على الحدود ~ ξ 4 s 1 − 2 و ~ ξ 4 s 2 − 2 ، بينما الحدود المتبقية في المعادلة 117 هي ~ ξ −2 . لذلك، يُسمح بتطبيق المعادلة 103 بدلاً من المعادلة 117 في حالة الشرطين s 1 > 0 و s 2 > 0؛ ومن ثم 1 −s12-s22{\displaystyle s_{1}^{2}-s_{2}^{2}}> 0.

وبالتالي، عند قيم ξ الصغيرة، تتوقف تذبذبات الدالة γ ab بينما تبدأ الدالة γ 33 في الازدياد عند قيم ξ المتناقصة. هذا نمط كاسنر، وعند مقارنة γ 33 بـ γ ab ، لا ينطبق التقريب المذكور أعلاه.

وللتحقق من مدى توافق هذا التحليل، قامت شركة BKL بدراسة المعادلاتRα0{\displaystyle R_{\alpha }^{0}}= 0,Rα3{\displaystyle R_{\alpha }^{3}}= 0، وبحساب المكونات γ a 3 منها ، تأكدت صحة المتباينة 94. وقد أظهرت هذه الدراسة [ 52 ] أن المكونات γ a 3 في كلتا المنطقتين التقاربيتين كانت ≈ γ 33. لذا، فإن صحة المتباينة 93 تستلزم مباشرةً صحة المتباينة 94 .

يحتوي هذا الحل، كما ينبغي أن يكون عليه الحال في الحالة العامة لحقل في الفراغ، على أربع دوال اختيارية للإحداثيات المكانية الثلاثة x و y و z . في المنطقة ξ{\displaystyle \ll }1 هذه الدوال هي، على سبيل المثال، lect 1 , lect 2 , μ 1 , s 1 . في المنطقة ξ{\displaystyle \gg }1 يتم تعريف الدوال الأربع بواسطة متسلسلة فورييه بواسطة الإحداثي z من المعادلة 115 بمعاملات هي دوال لـ x و y ؛ على الرغم من أن تحليل متسلسلة فورييه (أو التكامل؟) يميز فئة خاصة من الدوال، إلا أن هذه الفئة كبيرة بما يكفي لتشمل أي مجموعة فرعية محدودة من مجموعة جميع الشروط الأولية الممكنة.

يحتوي الحل أيضًا على عدد من الدوال الاختيارية الأخرى للإحداثيات x و y . تظهر هذه الدوال الاختيارية ثنائية الأبعاد ، بشكل عام، لأن العلاقات بين الدوال ثلاثية الأبعاد في حلول معادلات أينشتاين تفاضلية (وليست جبرية)، بغض النظر عن المشكلة الأعمق المتعلقة بالمعنى الهندسي لهذه الدوال. لم يحسب BKL عدد الدوال ثنائية الأبعاد المستقلة لأنه في هذه الحالة يصعب استخلاص استنتاجات قاطعة نظرًا لأن الدوال ثلاثية الأبعاد تُعرَّف بواسطة مجموعة من الدوال ثنائية الأبعاد (انظر [ 52 ] لمزيد من التفاصيل). [ ملاحظة 24 ]

وأخيراً، يواصل BKL إظهار أن الحل العام يحتوي على الحل الخاص الذي تم الحصول عليه أعلاه للنماذج المتجانسة.

باستبدال متجهات الأساس للفضاء المتجانس من النوع التاسع لبيانكي في المعادلة 7، يأخذ مقياس الزمكان لهذا النموذج الشكل التالي

عندما c 2{\displaystyle \ll }في المعادلة a² و، يمكن تجاهل في كل مكان باستثناء الحد dz² . للانتقال من الإطار المتزامن المستخدم في المعادلة 121 إلى إطار بشروط المعادلة 91 ، يتم إجراء التحويل dt = cd ξ/2 والتعويض zz /2. بافتراض أن χ ≡ ln( a / b ){\displaystyle \ll }1، نحصل من المعادلة 121 في التقريب الأول:

وبالمثل، باستخدام متجهات الأساس للفضاء المتجانس من النوع الثامن لبيانكي، نحصل على

وفقًا لتحليل الفضاءات المتجانسة أعلاه، في كلتا الحالتين ab = ξ (تبسيطًا)أ02{\displaystyle a_{0}^{2}}= ξ 0 ) و χ من المعادلة 51 ؛ يتم إعطاء الدالة c (ξ) بواسطة الصيغتين 53 و 61 ، على التوالي ، لنماذج النوعين التاسع والثامن.

يتم الحصول على مقياس مماثل للنوع الثامن من المعادلات 112 و 115 و 116 باختيار متجهين ثنائيي الأبعاد l و a و ma على النحو التالي :

واستبدال

للحصول على مقياس النوع التاسع، ينبغي استبدال

(لحساب c (ξ) فإن التقريب في المعادلة 116 غير كافٍ ويتم حساب الحد في ψ الخطي بواسطة ξ [ 56 ] )

أُجري هذا التحليل للفراغ. إن تضمين المادة لا يُقلل من عمومية الحل ولا يُغير من خصائصه النوعية. [ 56 ] [ 52 ]

من أهم القيود التي تعيق الحل العام أن جميع الدوال ثلاثية الأبعاد الواردة في المعادلتين 122 و 123 يجب أن يكون لها فاصل تغير مميز واحد ومشترك. هذا الشرط وحده يسمح بتقريب جميع مشتقات المكونات المكانية للمقياس في معادلات أينشتاين بضرب هذه المكونات ببساطة في أعداد موجية مميزة، مما ينتج عنه معادلات تفاضلية عادية من النوع المُستنتج للنموذج المتجانس من النوع التاسع. وهذا هو سبب التطابق بين الحلول المتجانسة والعامة.

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

الاستنتاجات

تصف نظرية BKL حالات التفرد في الحل الكوني لمعادلات أينشتاين، والتي تتسم بطابع تذبذبي معقد. ورغم أن هذه الحالات قد دُرست بشكل أساسي على نماذج متجانسة مكانيًا، إلا أن هناك أسبابًا وجيهة للافتراض بأن حالات التفرد في الحل العام لمعادلات أينشتاين لها نفس الخصائص؛ وهذا ما يجعل نموذج BKL ذا أهمية بالغة في علم الكونيات.

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

تم تأجيل دراسة المشكلة المتعلقة بالقيود المحتملة على هندسة الفضاء الناتجة عن التفرد. ومع ذلك، من الواضح منذ البداية أن نموذج BKL الأصلي قابل للتطبيق على كل من الفضاء المحدود وغير المحدود؛ ويتجلى ذلك في وجود نماذج تفرد تذبذبية لكل من الزمكان المغلق والمفتوح.

يُضفي النمط التذبذبي للاقتراب من التفرد بُعدًا جديدًا على مصطلح "محدودية الزمن". فبين أي لحظة زمنية محدودة في العالم t واللحظة t = 0، يوجد عدد لا نهائي من التذبذبات. وبهذا المعنى، تكتسب العملية طابعًا لانهائيًا. وبدلًا من الزمن t ، يُعدّ ln t متغيرًا أكثر ملاءمة لوصفها، حيث تُوسّع العملية لتشمل-{\displaystyle -\infty }.

تُراعي نظرية BKL تطور المقياس في اتجاه تناقص الزمن. معادلات أينشتاين متناظرة بالنسبة لإشارة الزمن، لذا فإن تطور المقياس في اتجاه تزايد الزمن ممكنٌ بنفس القدر. مع ذلك، يختلف هذان الحالتان اختلافًا جوهريًا، لأن الماضي والمستقبل ليسا متكافئين بالمعنى الفيزيائي. لا يمكن أن يكون للتفرد المستقبلي معنى فيزيائي إلا إذا كان ممكنًا في ظل ظروف ابتدائية عشوائية موجودة في لحظة سابقة. لا يتطابق توزيع المادة والحقول في لحظة ما من تطور الكون بالضرورة مع الشروط المحددة اللازمة لوجود حل خاص معين لمعادلات أينشتاين.

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

في هذا السياق، تبرز أهمية خاصية نموذج "ميكس ماستر" الذي عرضه ميسنر [ 58 ] والمتعلقة بانتشار الإشارات الضوئية. في النموذج المتناحي، يوجد "أفق ضوئي"، أي أنه في كل لحظة زمنية، توجد مسافة قصوى يستحيل عندها تبادل الإشارات الضوئية، وبالتالي، وجود علاقة سببية: لا يمكن للإشارة أن تصل إلى هذه المسافات منذ اللحظة t = 0.

يتم تحديد انتشار الإشارة بواسطة المعادلة ds = 0. في النموذج المتناحي بالقرب من نقطة التفرد t = 0، يكون عنصر الفاصل الزمني هودs2=دت2-2تدل¯2{\displaystyle ds^{2}=dt^{2}-2td{\bar {l}}^{2}}، أيندل¯2{\displaystyle d{\bar {l}}^{2}}هي صيغة تفاضلية مكانية مستقلة عن الزمن. [ 59 ] بالاستبدالت=η2/2{\displaystyle t=\eta ^{2}/2}العائد

"المسافة"Δل¯{\displaystyle \Delta {\bar {l}}}تصل الإشارة إلى

بما أن η، مثل t ، تتغير قيمها بدءًا من 0، فإن الإشارات حتى "اللحظة" η لا يمكنها الانتشار إلا على مسافةΔل¯η{\displaystyle \Delta {\bar {l}}\leq \eta }والذي يحدد أبعد مسافة إلى الأفق.

يُشكّل وجود أفق ضوئي في النموذج المتناحي إشكاليةً في فهم أصل التناحي الملحوظ حاليًا في الإشعاع المتبقي. فبحسب هذا النموذج، يعني التناحي الملحوظ خصائص متناحية للإشعاع القادم إلى الراصد من مناطق في الفضاء لا يمكن ربطها ببعضها البعض بشكل سببي. أما في نموذج التطور التذبذبي بالقرب من التفرد، فقد يختلف الوضع.

على سبيل المثال، في النموذج المتجانس للفضاء من النوع التاسع، تنتشر الإشارة في اتجاه تتغير فيه المقاييس لفترة طويلة وفقًا لقانون قريب من t . مربع عنصر المسافة في هذا الاتجاه هو dl² = ​​t²ل¯2{\displaystyle {\bar {l}}^{2}}، والعنصر المقابل في الفترة رباعية الأبعاد هودs2=دت2-ت2ل¯2{\displaystyle ds^{2}=dt^{2}-t^{2}{\bar {l}}^{2}}الاستبدالت=هـη{\displaystyle t=e^{\eta }}يضع هذا في النموذج

وبالنسبة لانتشار الإشارة، لدينا معادلة من النوع 128 مرة أخرى. والفرق المهم هو أن المتغير η يمر الآن بقيم تبدأ من-{\displaystyle -\infty }(إذا كانت المعادلة المترية 129 صالحة لجميع قيم t بدءًا من t = 0).

لذلك، بالنسبة لكل "لحظة" معينة η، يتم إيجاد فترات وسيطة Δη كافية لتغطية الإشارة لكل مسافة محدودة.

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

مع مرور الوقت والابتعاد عن نقطة التفرد، يتزايد تأثير المادة على تطور الأبعاد، والذي كان ضئيلاً في المراحل الأولى من التطور، تدريجياً حتى يصبح هو المهيمن. ومن المتوقع أن يؤدي هذا التأثير إلى "تجانس" تدريجي للفضاء، ما يجعل خصائصه أقرب إلى نموذج فريدمان الذي يصف بدقة الحالة الراهنة للكون.

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

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

ملحوظات

  1. يمكن العثور على محاكاة متحركة مماثلة من تصميم ديفيد غارفينكل في [ 1 ] .
  2. ١ ٢ ٣ يستخدم نظام BKL نفس الاصطلاح الوارد في كتاب لاندو وليفشيتز (١٩٨٨) . تُشير المؤشرات اللاتينية إلى القيم ٠، ١، ٢، ٣؛ بينما تُشير المؤشرات اليونانية إلى القيم المكانية ١، ٢، ٣. يأخذ المقياس g ik الإشارة (+ − − −)؛ γ αβ = − g αβ هو موتر المقياس المكاني ثلاثي الأبعاد. يستخدم نظام BKL نظام وحدات تكون فيه سرعة الضوء وثابت أينشتاين للجاذبية مساويين لـ ١.
  3. يتم اشتقاق التعبير عن r عن طريق لوغاريتم معاملات القوة في المقياس: ln [ t 2 p α (1/ u ) ] = 2 p α (1/ u ) ln t .
  4. عندما ( p 1 , p 2 , p 3 ) = (0, 0, 1) فإن مقياس الزمكان المعادلة 1 مع dl 2 من المعادلة 2 يتحول إلى مقياس غاليلي مع الاستبدال t sh z = ζ, t ch z = τ، أي أن التفرد وهمي والزمكان مسطح.
  5. هنا وفيما يلي، يجب فهم جميع رموز عمليات المتجهات (الضرب المتجهي، وعمليات الدوران، والتدرج، وما إلى ذلك) بطريقة رسمية للغاية على أنها عمليات على المكونات المتغيرة للمتجهات l و m و n بحيث يتم تنفيذها في الإحداثيات الديكارتية x 1 و x 2 و x 3 .
  6. باستثناء الحالة ( p 1 , p 2 , p 3 ) = (0, 0, 1)، حيث تكون التفردية المترية وهمية.
  7. الثوابت λ و μ و ν هي ما يسمى بالثوابت الهيكلية لمجموعة الحركة الفضائية.
  8. تحتوي معادلات أينشتاين للفضاء المتجانس، في صيغتها الدقيقة، عمومًا، على ست دوال زمنية مختلفة γ ab ( t ) في المقياس. إن الحصول على نظام متسق من المعادلات الدقيقة في هذه الحالة للمقياس الذي يحتوي على ثلاث دوال زمنية فقط (γ 11 = а 2 ، γ 22 = b 2 ، γ 33 = c 2 ) يرتبط بتناظر يؤدي إلى اختفاء ستة مكونات لموتر ريتشي.
  9. يمكن إيجادالقيم التقاربية لـ ατ و βτ و γτ عند τ → −∞ دون حل المعادلة 29 بالكامل . يكفي ملاحظة أن المعادلة الأولى تأخذ شكل "جسيم" يتحرك في بُعد واحد في مجال جدار جهد أُسّي، حيث يلعب α دور الثابت. في هذا التشبيه، يشير نمط كاسنر إلى حركة حرة بسرعة ثابتة ατ = Λp / 1 . بعد الانعكاس عن الجدار، يتحرك الجسيم بحرية بسرعة ατ = −Λp / 1 . مع ملاحظة أنه من المعادلة... 29 α τ + β τ = const، و α τ + γ τ = const، يمكن للمرء أن يرى أن β τ و γ τ يأخذان القيم β τ = Λ( p 2 − 2 p 1 )، γ τ = Λ( p 3 − 2 p 1 ).
  10. إن إدخال المكونات غير القطرية لـ γ ab ( t ) يضفي بعض الميزات الجديدة على نموذج BKL: دوران المحاور المقابلة لقوى حقبة كاسنر؛ تمت دراسة هذه المشكلة في Belinsky و Khalatnikov و Lifshitz (1971).
  11. إن الثابت في وسيط الجيب، بالطبع، ليس بالضرورة هو نفسه ξ 0 في المعادلة 47 والمعادلة 48 ؛ ومع ذلك، فإن جعلهما متطابقين لا يغير بأي حال من الأحوال طبيعة الحل.
  12. في حساب أكثر دقة، يظهر مصطلح لوغاريتمي يتغير ببطء في وسيط الجيب، ويظهر مضاعف أمام الأس في التعبير عن с (ξ)، انظر Belinsky، Khalatnikov & Lifshitz 1970 ، الملحق ب .
  13. إذا استبدلنا sh 2χ بـ 2χ في المعادلة 49 وحلناها لجميع قيم ξ، نحصل على χ = c 1 J 0 (ξ) + c 2 N 0 (ξ) حيث J 0 و N 0 دالتان من دوال بيسل من النوع الأول والثاني على التوالي. هذا الحل يربط بين الحالتين الحدّيتين، ويسمح بربط المعاملات الثابتة في المعادلتين 52 و 55 بترتيب تقريبي.
  14. بما أن a و b و c لها بُعد الطول، فإن لوغاريتماتها تُعرَّف فقط حتى ثابت إضافي يعتمد على اختيار وحدات الطول؛ وبهذا المعنى ، فإن المعادلة 63 لها معنى مشروط يتوافق مع اختيار معين للقيمة الصفرية لـ α و β و γ.
  15. وفقًا للمعادلة 32 ، تكون الانتقالات كبيرة مع قيمة صغيرة لـ | p1 | ( أي قيمة كبيرة لـ u ) وتكون ≈1/| p1 | ~ u . ولكن حتى في هذه الحالة Δn ~ u n | αn |{\displaystyle \gg }غير
  16. يُعدّ تحديد حدود الحقبة وفقًا للمعادلة 64 أمرًا ذا مغزى، لأنه في هذه الحالة، تشمل الحقبة جميع الفترات التي تتناقص فيها الدالة الثالثة، γ( t )، بشكل رتيب. إذا عُرّفت الحقبة بتسلسل قيم u من k + x إلى 1 + x ، فإن التناقص الرتيب لـ γ( t ) سيستمر خلال الفترة الأولى من الحقبة التالية.
  17. أطوال الحقب الزمنية كبيرة مقارنةً بالانتقالات بينها. وفقًا للمعادلة 33، تكون أطوال الانتقالات كبيرة عند قيم صغيرة لـ | p1 | (أي عند قيم كبيرة لـ u ) وتتناسب طرديًا مع 1/| p1 | ∝ u . ولكن حتى في هذه الحالة، فإن Δn u n | αn |{\displaystyle \gg }يو إن .
  18. ↑ كانت المعادلة 74 معروفةً بالفعل لغوس ، وقد تناول روديون كوزمين معادلةً من النوع 73c في هذا السياق (انظر توزيع غاوس-كوزمين ). للمزيد من المعلومات حول السلوك الفوضوي والإنتروبيا للكسور المستمرة، انظر: ليناس فيبستاس، 2008، إنتروبيا الكسور المستمرة (إنتروبيا غاوس-كوزمين).
  19. إن رسم الدالة P (δ) في الشكل 2 من كتاب ليفشيتز، ليفشيتز، وخالاتنيكوف 1970 غير صحيح لعدة أسباب. يبدو أنه تم إدخال بعض الأخطاء في إعداد برنامج الحل العددي للمعادلة التكاملية. كما تم إجراء اختزال "قسري" للقيمتين P (0) و P (1) نظرًا للحاشية غير الصحيحة في كتاب ليفشيتز، ليفشيتز ، وخالاتنيكوف 1970 ، القسم 4. إن الاحتمالية المحدودة للقيمة δ = 0 لا تعني إمكانية أن تصبح السعة الأولية للتذبذب صفرًا (وهو ما يتعارض مع المسار المنتظم للتطور الموضح في الشكل 4). من المعادلة 78، يؤولδ( s+1) إلى الصفر عندما x( s) → 0، وهو ما يتناسب مع x( s) . لكن السعة تعطى بواسطة حاصل ضرب δ s+1 Ω s+1 ، والذي يؤول إلى حد نهائي لأن التعبير eq 77 يحتوي على حد مع 1/ x s .
  20. لاحظ أن هذا المقياس يسمح بتحويلات عشوائية من النوع ξ′ + z ″ = f 1 (ξ + z ), ξ′ − z ′ = f 2 (ξ − z ), x a = f a ( x 1 , x 2 ).
  21. المعادلةR00+R33=0{\displaystyle R_{0}^{0}+R_{3}^{3}=0}هي نتيجة مباشرة للمعادلة 97 - المعادلة 99 إذاجي˙0{\displaystyle {\dot {G}}\neq 0}أوجي0{\displaystyle G^{\prime }\neq 0}القضيةجي˙=جي=0{\displaystyle {\dot {G}}=G^{\prime }=0}لا يتطلب الأمر معالجة خاصة: يمكن إثبات أن مقياس الزمكان في هذه الحالة يتقارب (في التقريب الأول) إلى مقياس غاليليو.
  22. من الممكن البحث عن حل في صورة تكاملات فورييه؛ لم تُدرس هذه المسألة بالتفصيل. لذلك، لا تشترط نظرية BKL تحليل متسلسلة فورييه كشرط إلزامي لاعتماد الدالتين σ و φ على الإحداثيات.
  23. تؤدي الحدود المربعة H في المعادلة 103 إلى تصحيحات صغيرة فقط (≈1/ξ) في σ و φ. أما الحساب باستخدام الحدود التكعيبية فيؤدي إلى ظهور اعتماد ضعيف لـ A و B على ξ، والذي يمكن تمثيله بظهور أطوار لوغاريتمية في العوامل المتذبذبة في المعادلة 115. هذه الحسابات لحالة ρ = 0 واردة في Belinsky & Khalatnikov (1970 ، الملحق ب) (انظر الحالة المماثلة للنماذج المتجانسة، Belinsky, Khalatnikov & Lifshitz (1970 ، الملحق ب) ).
  24. يحتوي التفكيك المنتظم للحل العام لمعادلات أينشتاين (بالإضافة إلى الدوال الأربع ثلاثية الأبعاد) على ثلاث دوال مستقلة لإحداثيين (انظر بتروف 1969 ، الفصل 40 ؛ ليفشيتز وخالاتنيكوف (1963 ، الملحق أ) )

مراجع

  1. غارفينكل، ديفيد (2007). "عن التفردات وصناعة الخبز" . أينشتاين أونلاين . المجلد 3. معهد ماكس بلانك لفيزياء الجاذبية . 03-1014 . تاريخ الاسترجاع: 15 أكتوبر 2020 .
  2. 1 2 3 بيلينسكي، خالاتنيكوف وليفشيتز 1970
  3. ديماريه، هينو وسبينديل 1985 .
  4. ديماريه وآخرون 1986 .
  5. ^ ديماريت، دي روب وهينو 1989 .
  6. دامور وهينو 2000 .
  7. دامور وآخرون 2001 .
  8. ^ دامور وهينو ونيكولاي 2003 .
  9. Kac 1983 .
  10. دامور 2015 .
  11. ^ هينو وبيرسون وسبيندل 2008 .
  12. 1 2 3 4 5 6 ليفشيتز وخالاتنيكوف 1963
  13. 1 2 لاندو وليفشيتز 1988 ، الفصل 97
  14. ^ ليفشيتز وخالاتنيكوف 1961 أ .
  15. ^ ليفشيتز و خالاتنيكوف 1961 ب .
  16. 1 2 ليفشيتز، سوداكوف وخالاتنيكوف 1961
  17. هوكينج 1965 .
  18. هوكينج وإيليس 1968 .
  19. جيروش 1966 .
  20. 1 2 أشتيكار، هندرسون وسلون 2011
  21. بارو وتيبلر 1979 .
  22. بارو وتيبلر 1981 .
  23. 1 2 بيرغر 2002
  24. غارفينكل 2004 .
  25. بيرغر ومونكريف 1993 .
  26. بيرغر وآخرون 1998 .
  27. ويفر، إيزنبرغ وبيرغر 1998 .
  28. ^ أندرسون وريندال 2001 .
  29. دامور وآخرون 2002 .
  30. بيرغر ومونكريف 1998 .
  31. بيرغر ومونكريف 2000 .
  32. غارفينكل 2007 .
  33. ساوتومي، أخوري وغارفينكل 2010 .
  34. كاسنر 1921 .
  35. روغ 1994 .
  36. ^ بيني وشيروبيني وجانتزين 2007 .
  37. Landau & Lifshitz 1988 ، الفصل 117 ، المسألة 3.
  38. Landau & Lifshitz 1987 ، الفصل 134 ، المعادلة 134.15.
  39. ميسنر، ثورن وويلر 1973 ، ص 564 . 
  40. نيلسون 1981 .
  41. بيلينسكي وخالاتنيكوف 1966 .
  42. خالاتنيكوف وليفشيتز 1970 .
  43. 1 2 بيلينسكي وخالاتنيكوف 1969أ
  44. 1 2 ليفشيتز وخالاتنيكوف 1970
  45. بيلينسكي، خالاتنيكوف وليفشيتز 1970 ، الملحق ج.
  46. ليفشيتز وخالاتنيكوف 1963 ، الملحق ج.
  47. تاوب 1951 .
  48. 1 2 3 4 ليفشيتز، ليفشيتز وخالاتنيكوف 1970
  49. 1 2 خالاتنيكوف وآخرون 1985
  50. تشيرنوف وبارو 1983 .
  51. بيلينسكي، خالاتنيكوف وليفشيتز 1970 ، الملحق أ.
  52. 1 2 3 4 بيلينسكي وخالاتنيكوف 1970
  53. أينشتاين وروزن 1937 .
  54. ^ بوندي وبيراني وروبنسون 1959 .
  55. Landau & Lifshitz 1988 ، الفصل 109.
  56. 1 2 3 بيلينسكي وخالاتنيكوف 1969ب
  57. بيلينسكي 1992 .
  58. ميسنر 1969 .
  59. Landau & Lifshitz 1988 ، الفصل 103–105.

فهرس