تفرد BKL

تُعدّ نقطة بيلينسكي -خالاتنيكوف-ليفشيتز ( 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 ] معادلة أينشتاين
| المعادلة 1 |
إن كتابة المعادلات في إطار متزامن تعطي نتيجة يصبح فيها المحدد المتري g صفراً حتماً في زمن محدود بغض النظر عن أي افتراضات حول توزيع المادة. [ 12 ] [ 13 ]
مع ذلك، تم التخلي عن الجهود المبذولة لإيجاد تفرد فيزيائي عام بعد أن اتضح أن التفرد المذكور أعلاه مرتبط بخاصية هندسية محددة للإطار المتزامن: تقاطع إحداثيات خط الزمن . يحدث هذا التقاطع على بعض الأسطح الفائقة المحيطة ، وهي نظائر رباعية الأبعاد للأسطح الكاوية في البصريات الهندسية ؛ وتصبح قيمة g صفرًا عند هذا التقاطع تحديدًا. [ 16 ] لذلك، على الرغم من أن هذا التفرد عام، إلا أنه وهمي وليس فيزيائيًا؛ إذ يختفي عند تغيير الإطار المرجعي. ويبدو أن هذا الأمر قد ثبط عزيمة الباحثين عن مواصلة البحث في هذا الاتجاه.
مرت عدة سنوات قبل أن يعود الاهتمام بهذه المشكلة للظهور مجدداً عندما نشر بنروز ( 1965 ) نظرياته التي ربطت وجود حالة شاذة ذات طبيعة غير معروفة ببعض الافتراضات العامة جداً التي لا علاقة لها باختيار إطار مرجعي. وقد توصل هوكينغ [ 17 ] [ 18 ] وجيروش [ 19 ] لاحقاً إلى نظريات مماثلة (انظر نظريات بنروز-هوكينغ حول الحالات الشاذة ). وقد أدى ذلك إلى إحياء الاهتمام بالبحث عن حلول شاذة.
حل متجانس معمّم
في فضاء متجانس ومتناحٍ، تُحدد المترية بشكل كامل، ولا يبقى سوى إشارة الانحناء. بافتراض تجانس الفضاء فقط دون أي تناظر إضافي كالتناحٍ، تتوفر حرية أكبر بكثير في اختيار المترية. يتعلق ما يلي بالجزء المكاني من المترية عند لحظة زمنية معينة t، بافتراض إطار مرجعي متزامن بحيث يكون 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 ]
حل كاسنر

يبدأ منهج BKL للفضاءات المتجانسة غير المتناحية (على عكس المتناحية) بتعميم حل خاص دقيق اشتقه كاسنر [ 34 ] لحقل في الفراغ، حيث يكون الفضاء متجانسًا وله مقياس إقليدي يعتمد على الزمن وفقًا لمقياس كاسنر.
| المعادلة 2 |
( dl هو عنصر الخط ؛ dx و dy و dz هي إزاحات متناهية الصغر في الأبعاد المكانية الثلاثة ، و t هو الفترة الزمنية المنقضية منذ لحظة ابتدائية t₀ = 0 ). هنا، p₁ و p₂ و p₃ هي أي ثلاثة أعداد تحقق شروط كاسنر التالية
| المعادلة 3 |
بسبب هذه العلاقات، يكون عدد واحد فقط من الأعداد الثلاثة مستقلاً ( معادلتان بثلاثة مجاهيل ). لا تتطابق الأعداد الثلاثة أبدًا؛ يتطابق عددان فقط في مجموعات القيم.و(0، 0، 1). [ ملاحظة 4 ] في جميع الحالات الأخرى، تكون الأرقام مختلفة، أحدها سالب والآخران موجبان. وقد تم إثبات ذلك جزئيًا بتربيع طرفي المعادلة 3 ( الشرط الأول ) ثم فك التربيع.
على المدىيساوي 1 وفقًا للشرط الثاني في المعادلة 3، وبالتالي يجب أن يكون الحد الذي يحتوي على نواتج مختلطة مساويًا للصفر. وهذا ممكن إذا كان أحد الحدود p1 أو p2 أو p3 على الأقل سالبًا.
إذا تم ترتيب الأرقام بترتيب تصاعدي، p 1 < p 2 < p 3 ، فإنها تتغير في الفترات (الشكل 4).
| المعادلة 4 |

تُقابل معادلة كاسنر ( 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 ] ) على النحو التالي
| المعادلة 5 |
تبدو معلمات مؤشر كاسنر غامضة حتى نتأمل في القيدين المفروضين على المؤشرات (المعادلة 3) . يحدد هذان القيدان المقياس الكلي للمؤشرات بحيث لا تتغير إلا نسبها . من الطبيعي اختيار إحدى هذه النسب كمعلمة جديدة، وهو ما يمكن فعله بست طرق مختلفة. على سبيل المثال، باختيار u = u₃² = p₃ / p₂ ، يصبح من السهل التعبير عن جميع النسب الست الممكنة بدلالة هذه النسبة. بحذف p₃ = u₂² أولًا ، ثم باستخدام القيد الخطي لحذف p₁ = 1 − p₂ − u₂² = 1 − (1 + u ) p₂ ، يختزل القيد التربيعي إلى معادلة تربيعية في p₂
| المعادلة 5أ |
بجذور 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 فتُوضع ضمن هذا النطاق وفقًا لـ
| المعادلة 6 |
في الحل المعمم، ينطبق الشكل المقابل للمعادلة 2 فقط على المقياس التقاربي (المقياس القريب من النقطة الشاذة t = 0)، وعلى الحدود الرئيسية لتوسيع متسلسلته بدلالة قوى t . في إطار المرجع المتزامن، يُكتب على شكل المعادلة 1 مع عنصر مسافة مكانية .
| المعادلة 7 |
أين
| المعادلة 8 |
تُحدد المتجهات ثلاثية الأبعاد 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 هو
| المعادلة 9 |
حيث v = l [ mn ]. من الملائم إدخال الكميات التالية [ ملاحظة 5 ]
| المعادلة 10 |
إن المقياس المكاني في المعادلة 7 غير متناحٍ لأن قوى t في المعادلة 8 لا يمكن أن تتساوى قيمها. عند الاقتراب من نقطة التفرد عند t = 0، تتناقص المسافات الخطية في كل عنصر مكاني في اتجاهين وتزداد في الاتجاه الثالث. ويتناقص حجم العنصر بما يتناسب مع t .
يتم إدخال مقياس كاسنر في معادلات أينشتاين عن طريق استبدال موتر المقياس γ αβ من المعادلة 7 دون تحديد مسبق لاعتماد a و b و c على t : [ ملاحظة 2 ]
حيث تشير النقطة فوق الرمز إلى التفاضل بالنسبة للزمن. تأخذ معادلة أينشتاين (المعادلة 11) الشكل التالي:
| المعادلة 14 |
جميع حدودها من الرتبة الثانية بالنسبة للكمية الكبيرة (عند t → 0) 1/ t . في معادلات أينشتاين (المعادلة 12) ، لا تظهر الحدود من هذه الرتبة إلا من الحدود المشتقة زمنيًا. إذا لم تتضمن مركبات P αβ حدودًا من رتبة أعلى من الثانية، فإن
| المعادلة 15 |
حيث تشير المؤشرات 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 كما يلي:
| المعادلة 16 |
هنا، الحدود الثانية من الرتبة t −2( p m + p n − p l ) حيث p m + p n − p 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 ، وتُضيف أربع دوال إحداثية جديدة ضرورية لوصف التوزيع الأولي لكثافتها ومكونات سرعتها الثلاثة. وهذا يُتيح تحديد تطور المادة انطلاقًا من قوانين حركتها في مجال جاذبية مُعطى مسبقًا، وهي المعادلات الهيدروديناميكية .
| المعادلة 17 |
| المعادلة 18 |
حيث uᵢ هي السرعة رباعية الأبعاد، و ε وσ هما كثافتا الطاقة والإنتروبيا للمادة (انظر [ 37 ] و [ 38 ] ، وكذلك [ 39 ] ، وللمزيد من التفاصيل انظر [ 40 ] ). بالنسبة لمعادلة الحالة النسبية الفائقة p = ε/3، فإن الإنتروبيا σ ~ ε 1/4 . الحدود الرئيسية في المعادلتين 17 و 18 هي تلك التي تحتوي على مشتقات زمنية . من المعادلة 17 والمكونات المكانية للمعادلة 18، نحصل على
مما أدى إلى
| المعادلة 19 |
حيث تمثل 'const' كميات مستقلة عن الزمن. بالإضافة إلى ذلك، من المتطابقة u i u i = 1، نحصل على (لأن جميع المكونات المتغيرة لـ u α من نفس الرتبة)
حيث يمثل u <sub>n</sub> مركبة السرعة على طول اتجاه n المرتبطة بأعلى قوة (موجبة) لـ t (بافتراض أن p <sub>n</sub> = p<sub> 3 </sub>). من العلاقات المذكورة أعلاه، يتبين أن
| المعادلة 20 |
أو
| المعادلة 21 |
يمكن استخدام المعادلات المذكورة أعلاه للتأكد من أن مكونات موتر الإجهاد والطاقة والزخم للمادة الموجودة في الجانب الأيمن من المعادلات
في الواقع، هي أقل رتبة بمقدار 1/ t من الحدود الرئيسية في طرفيها الأيسر. في المعادلاتإن وجود المادة لا يؤدي إلا إلى تغيير العلاقات المفروضة على وظائف الإحداثيات المكونة لها. [ 12 ]
إن حقيقة أن قيمة ε تصبح لانهائية وفقًا للمعادلة 21 تؤكد أنه في حل المعادلة 7، نتعامل مع حالة شاذة فيزيائية عند جميع قيم القوى p1 و p2 و p3 باستثناء ( 0 ، 0 ، 1 ) . بالنسبة لهذه القيم الأخيرة، تكون الحالة الشاذة غير فيزيائية ويمكن إزالتها بتغيير الإطار المرجعي.
تنشأ النقطة الشاذة الوهمية الموافقة للقوى (0، 0، 1) نتيجةً لعبور إحداثيات خط الزمن فوق سطح بؤري ثنائي الأبعاد . وكما أشير في المرجع [ 12 ] ، يمكن دائمًا اختيار إطار مرجعي متزامن بحيث يحدث هذا العبور الحتمي لخط الزمن تحديدًا على هذا السطح (بدلًا من سطح كاوستيك ثلاثي الأبعاد). لذلك، يجب أن يوجد حلٌّ يتضمن هذه النقطة الشاذة الوهمية المتزامنة للفضاء بأكمله، مع مجموعة كاملة من الدوال الاختيارية اللازمة للحل العام. بالقرب من النقطة t = 0، يسمح ذلك بتوسع منتظم بقوى صحيحة لـ t . لتحليل هذه الحالة، انظر المرجع [ 41 ] .
نمط التذبذب نحو التفرد
الحل العام، بحكم تعريفه، مستقر تمامًا؛ وإلا لما وُجد الكون. أي اضطراب يُعادل تغييرًا في الشروط الابتدائية في لحظة زمنية معينة؛ ولأن الحل العام يسمح بشروط ابتدائية اختيارية، فإن الاضطراب لا يستطيع تغيير طبيعته. من هذا المنظور، فإن الشروط الأربعة المفروضة على دوال الإحداثيات في معادلة الحل 7 هي من أنواع مختلفة: ثلاثة شروط تنشأ من المعادلات.تُعتبر القيم λ = 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 ]
| المعادلة 22 |
| المعادلة 23 |
(المكونات المتبقية)،،،،،(أصفار متطابقة). تحتوي هذه المعادلات على دوال للزمن فقط؛ وهذا شرط يجب تحقيقه في جميع الفضاءات المتجانسة. هنا، المعادلتان 22 و 23 دقيقتان، ولا تعتمد صحتهما على مدى قربنا من نقطة التفرد عند t = 0. [ ملاحظة 8 ]
تأخذ المشتقات الزمنية في المعادلتين 22 و 23 شكلاً أبسط إذا تم استبدال а و b و c بلوغاريتماتها α و β و γ:
| المعادلة 24 |
باستبدال المتغير t بالمتغير τ وفقًا لما يلي:
| المعادلة 25 |
ثم (تشير الرموز السفلية إلى التفاضل بواسطة τ):
| المعادلة 26 |
| المعادلة 27 |
بجمع المعادلات 26 واستبدال المجموع (α + β + γ) τ τ في الطرف الأيسر وفقًا للمعادلة 27 ، نحصل على معادلة تحتوي فقط على المشتقات الأولى وهي التكامل الأول للنظام 26 :
| المعادلة 28 |
تؤدي هذه المعادلة دور شرط مُلزم مفروض على الحالة الابتدائية للمعادلة 26. يُعد نمط كاسنر (المعادلة 8) حلاً للمعادلة 26 عند إهمال جميع الحدود في الطرف الأيمن. لكن لا يمكن أن يستمر هذا الوضع (عند t → 0) إلى ما لا نهاية، لأنه من بين هذه الحدود، توجد دائمًا حدود تتزايد. بالتالي، إذا كانت القوة السالبة في الدالة a ( t ) (حيث p<sub> l </sub> = p<sub> 1 </sub>)، فإن اضطراب نمط كاسنر سينشأ من الحدود λ<sub> 2</sub> a<sub> 4</sub> ؛ أما بقية الحدود فستتناقص مع تناقص t . إذا اقتصرت الحدود المتزايدة فقط على الطرف الأيمن من المعادلة 26 ، نحصل على النظام التالي:
| المعادلة 29 |
(قارن المعادلة 16 ؛ حيث تم استبدال λ 2 = 1 أدناه). يجب أن يصف حل هذه المعادلات تطور المقياس من الحالة الابتدائية، والتي يتم وصفها بالمعادلة 8 مع مجموعة معينة من القوى (مع p l < 0)؛ لنفترض أن p l = p 1 ، و p m = p 2 ، و p n = p 3 بحيث
| المعادلة 30 |
ثم
| المعادلة 31 |
حيث Λ ثابت. تُعاد تعريف الشروط الابتدائية للمعادلة 29 على النحو التالي:
| المعادلة 32 |
يمكن تكامل المعادلات 29 بسهولة؛ والحل الذي يحقق الشرط في المعادلة 32 هو
| المعادلة 33 |
حيث أن b 0 و c 0 ثابتان آخران.
من السهل ملاحظة أن السلوك التقاربي للدوال في المعادلة 33 عند t → 0 هو المعادلة 30. والتعبيرات التقاربية لهذه الدوال والدالة t (τ) عند τ → −∞ هي [ ملاحظة 9 ]
عند التعبير عن a و b و c كدوال لـ t ، نحصل على
| المعادلة 34 |
أين
| المعادلة 35 |
ثم
| المعادلة 36 |
يوضح ما سبق أن الاضطراب يعمل بطريقة تُغير نمط كاسنر بنمط آخر، وفي هذه العملية ينقلب الأس السالب لـ t من الاتجاه l إلى الاتجاه m : فإذا كان p'l < 0 سابقًا، فإنه يصبح الآن p'm < 0. خلال هذا التغيير ، تمر الدالة a ( t ) بقيمة عظمى، وتمر الدالة b ( t ) بقيمة صغرى؛ b ، التي كانت متناقصة سابقًا، تزداد الآن؛ و a، التي كانت متزايدة ، تصبح متناقصة؛ و c ( t ) المتناقصة تستمر في التناقص. الاضطراب نفسه (λ²a⁴α في المعادلة 29 )، الذي كان متزايدًا سابقًا، يبدأ الآن في التناقص والتلاشي. وبالمثل ، يؤدي التطور اللاحق إلى زيادة في الاضطراب من الحدود التي تحتوي على μ² ( بدلاً من λ² ) في المعادلة 26 ، ثم تغيير نمط كاسنر، وهكذا.
من الملائم كتابة قاعدة استبدال القوة (المعادلة 35) بمساعدة المعادلة 5 الخاصة بالمعاملات :
| المعادلة 37 |
تبقى القوة الأكبر من بين القوتين الموجبتين موجبة.
يُطلق 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 . عند قيمة ابتدائية عشوائية ( غير نسبية ) لـ u، تستمر عملية التناوب هذه بلا حدود. [ ملاحظة 10 ]
في الحل الدقيق لمعادلات أينشتاين، تفقد القوى p <sub>l</sub> و p <sub> m</sub> و p<sub> n</sub> معناها الدقيق الأصلي. يُدخل هذا الأمر نوعًا من "الغموض" في تحديد هذه الأعداد (وبالتالي، في تحديد المعامل u )، وهو ما يجعل تحليل أي قيم محددة (مثل القيم النسبية ) لـ u عديم الجدوى، على الرغم من صغره. لذلك، فإن القوانين التي تتعلق بالقيم غير النسبية العشوائية لـ u هي فقط التي تحمل معنىً خاصًا.
تُسمى الفترات الأطول التي تتذبذب فيها مقاييس المسافات المكانية على طول محورين بينما تتناقص المسافات على طول المحور الثالث بشكل رتيب، بالعصور ؛ وتتناقص الأحجام وفقًا لقانون قريب من t . عند الانتقال من عصر إلى آخر، ينعكس اتجاه تناقص المسافات بشكل رتيب من محور إلى آخر. يكتسب ترتيب هذه الانتقالات طابعًا تقاربيًا لعملية عشوائية . ويُعدّ الترتيب العشوائي نفسه سمة مميزة لتناوب أطوال العصور المتتالية (يقصد BKL بطول العصر عدد حقب كاسنر التي يحتويها العصر، وليس فترة زمنية).
لكل حقبة ( الحقبة s ) سلسلة من قيم المعامل u بدءًا من القيمة الأكبر،ومن خلال القيم-1،-2، ...، وصولاً إلى أصغرها،< 1. ثم
| المعادلة 41 |
أي أن k ( s ) = [حيث تشير الأقواس إلى الجزء الكامل من القيمة. يمثل الرقم k ( s ) طول الحقبة، ويُقاس بعدد حقب كاسنر التي تحتويها الحقبة. بالنسبة للحقبة التالية
| المعادلة 42 |
في سلسلة الأرقام اللانهائية u ، المؤلفة بهذه القواعد ، توجد قيم صغيرة للغاية (ولكنها لا تساوي الصفر أبدًا) x ( s ) وأطوال كبيرة لا نهائية مقابلة k ( s ) .
تصبح سلسلة الحقبة أكثر كثافة عند الاقتراب من t = 0. ومع ذلك، فإن المتغير الطبيعي لوصف المسار الزمني لهذا التطور ليس الوقت العالمي t ، بل لوغاريتمه، ln t ، والذي من خلاله يتم تمديد العملية الكاملة للوصول إلى التفرد إلى −∞.
وفقًا للمعادلة 33 ، فإن إحدى الدوال a و b و c التي تمر بقيمة عظمى خلال الانتقال بين حقب كاسنر، عند ذروة قيمتها العظمى هي
| المعادلة 38 |
حيث يُفترض أن قيمة a max كبيرة مقارنةً بـ b 0 و c 0 ؛ في المعادلة 38، u هي قيمة المعامل في حقبة كاسنر قبل الانتقال. يتضح من هنا أن قمم القيم العظمى المتتالية خلال كل حقبة تنخفض تدريجيًا. في الواقع، في حقبة كاسنر التالية، تكون قيمة هذا المعامل u' = u − 1، ويتم استبدال Λ وفقًا للمعادلة 36 بالمعادلة Λ' = Λ(1 − 2| p 1 ( u )|). لذلك، فإن نسبة قيمتين عظميين متتاليتين هي
وأخيراً
| المعادلة 39 |
ما سبق هو حلول لمعادلات أينشتاين في الفراغ. أما بالنسبة لنمط كاسنر النقي، فإن المادة لا تُغير الخصائص النوعية لهذا الحل، ويمكن كتابته ضمنه دون الأخذ في الاعتبار تأثيرها على المجال. مع ذلك، إذا طُبِّق هذا على النموذج قيد البحث، باعتباره حلاً دقيقاً لمعادلات أينشتاين، فإن الصورة الناتجة لتطور المادة لن تكون عامة، بل ستكون خاصة بالتناظر العالي الذي يميز هذا النموذج. رياضياً، ترتبط هذه الخصوصية بحقيقة أنه بالنسبة لهندسة الفضاء المتجانسة التي نناقشها هنا، فإن مركبات موتر ريتشيتكون هذه القيم أصفارًا متطابقة، وبالتالي فإن معادلات أينشتاين لا تسمح بحركة المادة (التي تُنتج مركبات غير صفرية لموتر طاقة-زخم الإجهاد).بمعنى آخر ، يجب أن يكون الإطار المتزامن متحركًا أيضًا مع المادة. إذا استبدلنا في المعادلة 19 uα = 0 و u0 = 1، فإنها تصبح ε ~ ( abc ) −4/3 ~ t −4/3 .
يمكن تجنب هذه الصعوبة إذا اقتصر النموذج على تضمين الحدود الرئيسية للمقياس الحدي (عند t → 0) وكتابة مادة ذات توزيع أولي عشوائي للكثافات والسرعات. عندئذٍ، يتحدد مسار تطور المادة بقوانين حركتها العامة ، المعادلتين 17 و 18، والتي تُنتج المعادلة 21. خلال كل حقبة من حقب كاسنر، تزداد الكثافة وفقًا للقانون.
| المعادلة 40 |
حيث أن p 3 هو، كما سبق، أكبر الأعداد p 1 و p 2 و p 3. تزداد كثافة المادة بشكل رتيب خلال جميع مراحل التطور نحو التفرد.
تطور المقاييس
تتوافق قيم u الكبيرة جدًا مع قوى كاسنر
| المعادلة 43 |
وهي قريبة من القيم (0، 0، 1). القيمتان القريبتان من الصفر تكونان متقاربتين أيضًا، وبالتالي فإن التغيرات في نوعين من أنواع "الاضطرابات" الثلاثة (الحدود التي تحتوي على λ وμ وν في الطرف الأيمن من المعادلة 26 ) تكون متشابهة جدًا. إذا كانت هذه الحدود متقاربة جدًا في قيمها المطلقة في بداية حقبة طويلة كهذه، لحظة الانتقال بين حقبتين من حقب كاسنر (أو تم جعلها متقاربة بشكل مصطنع عن طريق تحديد الشروط الابتدائية)، فإنها ستظل متقاربة خلال معظم مدة الحقبة بأكملها. في هذه الحالة (يُطلق عليها BKL حالة التذبذبات الصغيرة )، يصبح التحليل القائم على تأثير نوع واحد من الاضطرابات غير صحيح؛ إذ يجب مراعاة التأثير المتزامن لنوعين من الاضطرابات.
اضطرابان
لنفترض حقبة زمنية طويلة، تشهد خلالها دالتان a و b و c (لنفترض أنهما a و b ) تذبذبات طفيفة، بينما تتناقص الدالة الثالثة ( c ) بشكل رتيب. سرعان ما تصبح الدالة الأخيرة صغيرة؛ لننظر في الحل ضمن المنطقة التي يمكن فيها إهمال c مقارنةً بـ a و b . تُجرى الحسابات أولًا لنموذج الفضاء من النوع التاسع بالتعويض وفقًا لذلك λ = μ = ν = 1. [ 43 ]
بعد إهمال الدالة c ، تعطي المعادلتان الأوليان (المعادلة 26)
| المعادلة 44 |
| المعادلة 45 |
ويمكن استخدام المعادلة 28 كمعادلة ثالثة، تأخذ الشكل التالي:
| المعادلة 46 |
يُكتب حل المعادلة 44 على الصورة التالية:
حيث α₀ و ξ₀ ثابتان موجبان، و τ₀ هي الحد الأعلى للعصر للمتغير τ. ومن الملائم إدخال متغير جديد (بدلاً من τ).
| المعادلة 47 |
ثم
| المعادلة 48 |
يتم تحويل المعادلتين 45 و 46 عن طريق إدخال المتغير χ = α − β:
| المعادلة 49 |
| المعادلة 50 |
انخفاض قيمة τ من τ₀ إلى −∞ يقابله انخفاض قيمة ξ من ξ₀ إلى 0. وتُحقق الفترة الطويلة ذات القيمتين a و b المتقاربتين (أي ذات قيمة χ صغيرة)، التي نتناولها هنا، إذا كانت ξ₀ قيمة كبيرة جدًا. في الواقع، عند قيم ξ الكبيرة، يكون حل المعادلة 49 في التقريب الأول باستخدام 1/ξ هو
| المعادلة 51 |
حيث A ثابت؛ المضاعفيجعل χ كمية صغيرة بحيث يمكن استبدالها في المعادلة 49 بـ sh 2χ ≈ 2χ. [ ملاحظة 11 ]
من المعادلة 50 نحصل على
بعد تحديد α و β من المعادلة 48 والمعادلة 51 وتوسيع e α و e β في سلسلة وفقًا للتقريب أعلاه، نحصل في النهاية على: [ ملاحظة 12 ]
| المعادلة 52 |
| المعادلة 53 |
تُستنتج العلاقة بين المتغير ξ والزمن t من خلال تكامل التعريف dt = abc d τ والذي يعطي
| المعادلة 54 |
يجب أن تكون قيمة الثابت c₀ ( قيمة c عند ξ = ξ₀ ) الآن c₀α 0 ·

لننظر الآن في المجال ξ1. فيما يلي أهم الحدود في حل المعادلة 49 :
حيث k ثابت في النطاق −1 < k < 1؛ يضمن هذا الشرط أن يكون الحد الأخير في المعادلة 49 صغيرًا (sh 2χ يحتوي على ξ 2 k و ξ −2 k ). بعد ذلك، وبعد تحديد α و β و t ، نحصل على
| المعادلة 55 |
هذا نمط كاسنر آخر مع وجود الأس السالب t في الدالة c ( t ). [ ملاحظة 13 ]
تُصوّر هذه النتائج تطورًا مشابهًا نوعيًا لما وُصف أعلاه. خلال فترة زمنية طويلة تتوافق مع انخفاض كبير في قيمة ξ، تتذبذب الدالتان a و b ، وتبقى قيمتاهما متقاربتين في المقدار.في الوقت نفسه، تتحرك كلتا الدالتين أ و ب ببطء (تتناقص الدالة. دورة التذبذبات ثابتة بالنسبة للمتغير ξ : Δξ = 2π (أو، وهو نفسه، بدورة ثابتة بالنسبة للزمن اللوغاريتمي: Δ ln t = 2π Α 2 ). الدالة الثالثة، c ، تتناقص بشكل رتيب وفقًا لقانون قريب من c = c 0 t / t 0 .
يستمر هذا التطور حتى تصبح قيمة ξ ≈ 1، وعندها لا تعود المعادلتان 52 و 53 قابلتين للتطبيق. وتتوافق مدته الزمنية مع تغير قيمة t من t₀ إلى القيمة t₁ ، المرتبطة بـ ξ₀ وفقًا لـ
| المعادلة 56 |
يمكن تمثيل العلاقة بين ξ و t خلال هذه الفترة بالشكل التالي:
| المعادلة 57 |
بعد ذلك، وكما هو موضح في المعادلة 55 ، تبدأ الدالة المتناقصة c بالزيادة بينما تبدأ الدالتان a و b بالتناقص. تستمر هذه المرحلة من كاسنر حتى تصبح الحدود c² / a²b² في المعادلة 22 قريبة من t² ، وتبدأ سلسلة جديدة من التذبذبات.
يتم الحصول على قانون تغير الكثافة خلال الحقبة الطويلة قيد المناقشة عن طريق استبدال المعادلة 52 في المعادلة 20 :
| المعادلة 58 |
عندما تتغير قيمة ξ من ξ 0 إلى ξ ≈1، تزداد الكثافةمرات.
يجب التأكيد على أنه على الرغم من أن الدالة c ( t ) تتغير وفقًا لقانون قريب من c ~ t ، فإن المقياس في المعادلة 52 لا يتوافق مع مقياس كاسنر ذي القوى (0، 0، 1). يتوافق هذا الأخير مع حل دقيق وجده تاوب [ 47 ] وهو مسموح به في المعادلتين 26 و 27 ، وفيه
| المعادلة 59 |
حيث p و δ1 و δ2 ثوابت . في المنطقة التقاربية τ → −∞، يمكن الحصول من هنا على a = b = ثابت، c = ثابت. t بعد التعويض е рτ = t . في هذا المقياس، تكون النقطة الشاذة عند t = 0 غير فيزيائية.
لنصف الآن الدراسة المماثلة لنموذج النوع الثامن، مع استبدال المعادلات 26 – 28 λ = −1، μ = ν = 1. [ 44 ]
إذا كانت الدالة المتناقصة بشكل رتيب هي a خلال الفترة الطويلة، فلن يتغير شيء في التحليل السابق: بتجاهل a² في الطرف الأيمن من المعادلتين 26 و 28 ، نعود إلى المعادلتين 49 و 50 نفسيهما (مع تغيير في الترميز). مع ذلك، تحدث بعض التغييرات إذا كانت الدالة المتناقصة بشكل رتيب هي b أو c ؛ لنفترض أنها c .
كما في السابق، لدينا المعادلة 49 بنفس الرموز، وبالتالي، التعبيرات السابقة في المعادلة 52 للدالتين a (ξ) و b (ξ)، ولكن يتم استبدال المعادلة 50 بـ
| المعادلة 60 |
يصبح الحد الرئيسي عند ξ الكبير الآن
لهذا السبب.
| المعادلة 61 |
قيمة c كدالة للزمن t هي c = c₀t / t₀ ، لكن اعتماد ξ على الزمن يتغير . يعتمد طول الحقبة الطويلة على ξ₀ وفقًا لـ
| المعادلة 62 |
من جهة أخرى، تحدد القيمة ξ₀ عدد تذبذبات الدالتين 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 Λ خلال حقبة زمنية معينة يُعطى بمجموع على الصورة التالية
مع كتابة الحدود التي تُقابل قيمًا كبيرة لـ u فقط. عندما تزداد قيمة k، يزداد هذا المجموع بمقدار ln k . لكن احتمال ظهور حقبة ذات طول كبير k يتناقص بمقدار 1/ k² وفقًا للمعادلة 76 ؛ وبالتالي فإن متوسط قيمة المجموع أعلاه محدود. ونتيجة لذلك، فإن التغير المنتظم للكمية ln Λ على مدى عدد كبير من الحقب سيكون متناسبًا مع هذا العدد. ولكن يُلاحظ في المعادلة 85 أنه مع t → 0، يزداد العدد s بمقدار ln |ln t | فقط. وهكذا، في الحد التقاربي لـ t الصغيرة جدًا ، يمكن بالفعل إهمال الحد ln Λ مقارنةً بـ ln t . في هذا التقريب [ ملاحظة 14 ]
| المعادلة 63 |
حيث تشير Ω إلى "الزمن اللوغاريتمي"
| المعادلة 64 |
ويمكن اعتبار عملية انتقال الحقب بمثابة سلسلة من ومضات زمنية قصيرة. كما تخضع مقادير القيم القصوى لدوال المقياس المتذبذبة لتغير منتظم. من المعادلة 39، بالنسبة لـ u ≫ 1، يتبين أنبنفس الطريقة التي تم بها سابقًا بالنسبة للكمية ln Λ، يمكن استنتاج أن متوسط انخفاض ارتفاع القيم العظمى خلال حقبة زمنية محددة محدود، وأن الانخفاض الكلي على مدى عدد كبير من الحقب الزمنية يزداد مع اقتراب t من الصفر تمامًا مثل ln Ω. في الوقت نفسه، فإن انخفاض القيم الصغرى، وبالتالي زيادة سعة التذبذبات ، يتناسب ( المعادلة 77 ) مع Ω. وفقًا للتقريب المعتمد، يُهمل انخفاض القيم العظمى مقارنةً بزيادة السعات، لذا فإن α max = 0، β max = 0، γ max = 0 للقيم العظمى لجميع الدوال المتذبذبة، وتمر الكميات α، β، γ بقيم سالبة فقط، وترتبط هذه القيم ببعضها البعض في كل لحظة زمنية بالعلاقة الموضحة في المعادلة 63 .

بالنظر إلى هذا التغير الفوري في الحقب، تُهمل فترات الانتقال لصغرها مقارنةً بطول الحقبة؛ وهذا الشرط مُتحقق بالفعل. [ ملاحظة 15 ] يتطلب استبدال القيم القصوى لـ α و β و γ بأصفار أن تكون الكميات ln (| p 1 |Λ) صغيرة مقارنةً بسعات تذبذبات الدوال المعنية. كما ذُكر سابقًا ، خلال الانتقالات بين الحقب، قد تصبح قيم | p 1 | صغيرة جدًا، بينما لا يرتبط مقدارها واحتمالية حدوثها بسعات التذبذب في اللحظة المعنية. لذلك، من حيث المبدأ، من الممكن الوصول إلى قيم صغيرة جدًا لـ | p 1 | بحيث يُنتهك الشرط المذكور أعلاه (القيم القصوى الصفرية). يمكن أن يؤدي هذا الانخفاض الحاد في α max إلى حالات خاصة مختلفة يصبح فيها الانتقال بين حقب كاسنر وفقًا للقاعدة المذكورة في المعادلة 37 غير صحيح (بما في ذلك الحالات الموصوفة أعلاه ). قد تُخالف هذه الحالات "الخطيرة" القوانين المستخدمة في التحليل الإحصائي أدناه. كما ذكرنا، فإن احتمال حدوث مثل هذه الانحرافات يتقارب بشكل مقارب إلى الصفر؛ وسيتم مناقشة هذه المسألة أدناه.
لنفترض حقبة زمنية تحتوي على k من حقب كاسنر مع وجود مُعامل u يمر عبر القيم
| المعادلة 65 |
ولنفترض أن α و β هما الدالتان المتذبذبتان خلال هذه الحقبة (الشكل 4). [ ملاحظة 16 ]
اللحظات الابتدائية لحقب كاسنر ذات المعاملات u n هي Ω n . في كل لحظة ابتدائية، تكون إحدى القيمتين α أو β تساوي صفرًا، بينما تكون الأخرى عند قيمة دنيا. القيمتان α أو β في القيم الدنيا المتتالية، أي في اللحظات Ω n، هما
| المعادلة 66 |
(دون تمييز بين القيمتين الصغريتين α و β). تتراوح قيم δn التي تقيس هاتين القيمتين الصغريتين بوحدات Ωn المقابلة بين 0 و1. تتناقص الدالة γ بشكل رتيب خلال هذه الفترة؛ ووفقًا للمعادلة 63، فإن قيمتها في اللحظة Ωn هي
| المعادلة 67 |
خلال الحقبة التي تبدأ عند اللحظة Ω n وتنتهي عند اللحظة Ω n +1 ، تزداد إحدى الدالتين α أو β من −δ n Ω n إلى الصفر بينما تنخفض الأخرى من 0 إلى −δ n +1 Ω n +1 وفقًا لقوانين خطية، على التوالي:
- و
مما يؤدي إلى علاقة تكرارية
| المعادلة 68 |
وبالنسبة لطول الحقبة اللوغاريتمية
| المعادلة 69 |
حيث، باختصار، f ( u ) = 1 + u + u² . ويُحسب مجموع أطوال n حقبة زمنية باستخدام الصيغة التالية :
| المعادلة 70 |
يتضح من المعادلة 68 أن |α n+1 | > |α n |، أي أن سعة تذبذب الدالتين α و β تزداد خلال الحقبة بأكملها، حتى وإن كانت العوامل δ n صغيرة. إذا كان الحد الأدنى في بداية الحقبة عميقًا، فلن تصبح الحدود الدنيا التالية أقل عمقًا؛ بعبارة أخرى، يبقى الباقي |α — β| عند لحظة الانتقال بين حقب كاسنر كبيرًا. لا يعتمد هذا على طول الحقبة k، لأن الانتقالات بين الحقب تُحدد بالقاعدة العامة في المعادلة 37 حتى بالنسبة للحقب الطويلة.
ترتبط سعة التذبذب الأخير للدالتين α أو β في حقبة زمنية معينة بسعة التذبذب الأول بالعلاقة |α k −1 | = |α 0 | ( k + x ) / (1 + x ). حتى عند قيم k الصغيرة جدًا، يمكن إهمال x مقارنةً بـ k، بحيث تصبح الزيادة في سعة تذبذب α و β متناسبة مع طول الحقبة الزمنية. بالنسبة للدالتين a = e α و b = e β، هذا يعني أنه إذا كانت سعة تذبذبهما في بداية حقبة زمنية ما هي A 0 ، فستصبح في نهاية هذه الحقبة A 0..
يزداد طول حقب كاسنر (بالزمن اللوغاريتمي) أيضًا داخل الحقبة الواحدة؛ ومن السهل حساب ذلك من المعادلة 69، حيث أن Δn + 1 > Δn . [ ملاحظة 17 ] إجمالي طول الحقبة هو
| المعادلة 71 |
(ينشأ الحد الذي يحتوي على 1/ x من الحقبة الأخيرة، الحقبة k ، التي يكون طولها كبيرًا عند قيم x الصغيرة ؛ انظر الشكل 2). اللحظة Ωn التي تنتهي عندها الحقبة k من حقبة معينة هي في نفس اللحظة Ω'0 لبداية الحقبة التالية.
في حقبة كاسنر الأولى من العصر الجديد ، تكون الدالة γ هي أول ما يرتفع من القيمة الدنيا γk = −Ωk (1 − δk ) التي بلغتها في الحقبة السابقة؛ وتلعب هذه القيمة دور سعة البداية δ'0 Ω'0 لسلسلة التذبذبات الجديدة. ومن السهل استنتاج ما يلي:
| المعادلة 72 |
من الواضح أن δ' 0 Ω' 0 > δ 0 Ω 0. حتى عند قيم k غير الكبيرة جدًا، يكون ازدياد السعة كبيرًا جدًا: تبدأ الدالة c = e γ بالتذبذب من السعة. يتم تأجيل مسألة الحالات "الخطيرة" المذكورة أعلاه لخفض حد التذبذب العلوي بشكل كبير في الوقت الحالي.
وفقًا للمعادلة 40، فإن الزيادة في كثافة المادة خلال الحقب الأولى ( k - 1) تُعطى بالصيغة التالية:
بالنسبة لآخر k حقبة من حقبة معينة، عند u = x < 1 ، تكون أعلى قيمة للقوة هي p2 ( x ) (وليس p3 ( x )). لذلك، بالنسبة لزيادة الكثافة على مدار الحقبة بأكملها، نحصل على
| المعادلة 73 |
لذلك، حتى عند قيم k غير الكبيرة جدًا ،خلال الحقبة التالية (بطول k ' ) ستزداد الكثافة بشكل أسرع بسبب زيادة السعة الابتدائية A 0 ':إلخ. توضح هذه الصيغ الزيادة الحادة في كثافة المادة.
التحليل الإحصائي بالقرب من نقطة التفرد
تكتسب متتالية أطوال الحقب k ( s ) ، المقاسة بعدد حقب كاسنر التي تحتويها، صفة العملية العشوائية بشكل تقاربي. وينطبق الأمر نفسه على متتالية تبادل أزواج الدوال المتذبذبة عند الانتقال من حقبة إلى أخرى (يعتمد ذلك على ما إذا كانت الأعداد k ( s ) زوجية أم فردية). أحد مصادر هذه العشوائية هو القاعدة الموضحة في المعادلتين 41 و 42 ، والتي تنص على أن الانتقال من حقبة إلى أخرى يتحدد في متتالية عددية لانهائية من قيم u . بعبارة أخرى، تنص هذه القاعدة على أنه إذا بدأت المتتالية اللانهائية بأكملها بقيمة ابتدائية معينة، فإن u = 0.إذاً، فإن أطوال الفترات الزمنية k (0) ، k (1) ، ...، هي الأعداد الموجودة في مفكوك الكسر المستمر البسيط.
| المعادلة 73أ |
يتوافق هذا التوسع مع تحويل الإسقاط للفترة [0, 1] على نفسها بالصيغة Tx = {1/ x }، أي x s +1 = {1/ x s }. ينتمي هذا التحويل إلى ما يُسمى بالتحويلات التوسعية للفترة [0, 1]، أي التحويلات x → f ( x ) حيث | f′ ( x )| > 1. تتميز هذه التحويلات بخاصية عدم الاستقرار الأسي: إذا أخذنا نقطتين متقاربتين في البداية، فإن المسافة بينهما تزداد أُسّيًا مع تكرار التحويلات. من المعروف أن عدم الاستقرار الأسي يؤدي إلى ظهور خصائص عشوائية قوية.
من الممكن الانتقال إلى وصف احتمالي لمثل هذا التسلسل من خلال النظر ليس إلى قيمة ابتدائية محددة x (0)، بل إلى القيم x (0) = x الموزعة في الفترة من 0 إلى 1 وفقًا لقانون توزيع احتمالي معين w₀ ( x ). عندئذٍ ، ستكون لقيم x (s) التي تُنهي كل حقبة توزيعات تتبع قوانين معينة wₛ ( x ) . لنفترض أن wₛ (x)dx هو احتمال أن تنتهي الحقبة s بالقيمة xₛ (x) .تقع ضمن نطاق محدد dx .
يمكن أن تنتج القيمة x (s) = x ، التي تُنهي الحقبة s ، من القيم الأولية (لهذه الحقبة).، حيث k = 1، 2، ...؛ هذه القيم منتُقابل هذه القيم x ( s -1) = 1/( k + x ) للفترة السابقة. وبناءً على ذلك، يمكن كتابة العلاقة التكرارية التالية، التي تُعبّر عن توزيع الاحتمالات w ( s ) بدلالة التوزيع w (s -1) ( x ):
أو
| المعادلة 73 ج |
إذا كان التوزيع w s ( x ) يميل مع ازدياد s إلى توزيع حدّي ثابت (مستقل عن s ) w ( x )، فإن هذا التوزيع الأخير يجب أن يحقق معادلة مستنتجة من المعادلة 73c بحذف مؤشرات الدالتين w s −1 ( x ) و w s ( x ). لهذه المعادلة حل
| المعادلة 74 |
(مُعَيَّر إلى الوحدة ومُقَسَّم إلى الرتبة الأولى من x ). [ ملاحظة 18 ]
لكي يكون طول الحقبة s هو k ، يجب أن تنتهي الحقبة السابقة برقم x يقع بين 1/( k + 1) و 1/ k . وبالتالي، فإن احتمال أن يكون طول الحقبة k يساوي (في حالة الاستقرار)
| المعادلة 75 |
عند قيم كبيرة لـ k
| المعادلة 76 |
عند ربط الخصائص الإحصائية للنموذج الكوني بالخصائص الإرجودية للتحويل 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 ) (حول أعداد محددة أخرى) سيكون مساويًا لـ
| المعادلة 76أ |
(هذا التعبير صالح فقط طالما أنه يحدد الكميات δ x (s) ≪ 1).
القيمة المتوسطة، المحسوبة من هذا التوزيع، تتباعد لوغاريتميًا. بالنسبة لمتتالية مقطوعة عند عدد كبير جدًا، ولكنه لا يزال محدودًا N ، يكون لديناإن فائدة المتوسط في هذه الحالة محدودة للغاية بسبب عدم استقراره: فبسبب التناقص البطيء لـ W ( k )، تتباعد تقلبات k أسرع من متوسطها. ومن الخصائص الأنسب لهذه المتتالية احتمال أن ينتمي عدد مختار عشوائيًا منها إلى حقبة طولها K، حيث K قيمة كبيرة. هذا الاحتمال هو ln K / ln 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 ، نحصل على
| المعادلة 76ب |
إن جمع هذا التوزيع على 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
| المعادلة 77 |
| المعادلة 78 |
| المعادلة 79 |
( تم إدخال ξ ( 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 ككسر مستمر لانهائي
| المعادلة 79أ |
(سهولة إدخال الترميز)(مع إزاحة الدليل بمقدار 1، سيتضح ذلك فيما يلي). وللاختصار، يُشار إلى الكسر المستمر ببساطة عن طريق تعداد مقاماته (بين قوسين مربعين)؛ ثم تعريفيمكن كتابتها على النحو التالي
| المعادلة 79ب |
تُعرَّف الكميات العكسية بكسر مستمر ذي مقامات رجعية (في اتجاه تناقص الأسس).
| المعادلة 79 ج |
يتم تحويل علاقة التكرار في المعادلة 78 بإدخال الرمز η s = (1 − δ s )/δ s مؤقتًا . وبذلك ، يمكن إعادة كتابة المعادلة 78 على النحو التالي:
يتم الحصول على كسر متصل لانهائي عن طريق التكرار
لذلكوأخيراً
| المعادلة 79د |
يحتوي هذا التعبير عن δ s على كميتين عشوائيتين فقط (بدلاً من الثلاث في [ 48 ] ).و، كل منها يأخذ قيمًا في الفترة [0، 1].
ويترتب على تعريف المعادلة 79c ما يلي:وبالتالي، فإن إزاحة المتتالية X بأكملها خطوة واحدة إلى اليمين تعني تحويلاً مشتركاً للكميات.ووفق
| المعادلة 79هـ |
هذا تحويل أحادي في المربع الواحدي . وبالتالي، لدينا الآن تحويل أحادي لكميتين بدلاً من تحويل غير أحادي Tx = {1/ x } لكمية واحدة.
الكمياتولدينا توزيع مشترك ثابت P ( x + , x- ). بما أن المعادلة 79e هي تحويل أحادي، فإن شرط ثبات التوزيع يُعبَّر عنه ببساطة بمعادلة دالة .
| المعادلة 79f |
حيث J هي جاكوبيان التحويل.
يؤدي تحريك المتتالية X بمقدار خطوة واحدة إلى التحويل التالي T للمربع الواحد:
(مع،انظر المعادلة 79هـ ). تُعرّف دالة الكثافة P ( x , y ) المقياس الثابت لهذا التحويل. من الطبيعي افتراض أن P ( x , y ) دالة متناظرة لـ x و y . هذا يعني أن المقياس ثابت بالنسبة للتحويل S ( x , y ) = ( y , x )، وبالتالي بالنسبة لحاصل الضرب ST مع ST ( x , y ) = ( x ″ , y ″ ) و
من الواضح أن ST لها تكامل أول H = 1/( x + y) . على الخط H = ثابت ≡ c ، يكون للتحويل الشكل التالي:
وبالتالي، يجب أن تكون كثافة القياس الثابت لـ ST على الشكل التالي
وبأخذ التناظر P ( x , y )= P ( y , x ) في الاعتبار، يصبح هذا f ( c )= c −2 وبالتالي (بعد التطبيع)
| ما يعادل 79 غرامًا |
(يؤدي تكاملها على x + أو x – إلى الدالة w ( x ) المعادلة 74 ). وقد استخدم تشيرنوف وبارو [ 50 ] اختزال التحويل إلى دالة أحادية، وحصلا على صيغة من شكل المعادلة 79g ولكن لمتغيرات أخرى؛ ولا تتضمن ورقتهما تطبيقات على المشكلات التي تناولها خالاتنيكوف وآخرون [ 49 ].
يمكن التحقق من صحة المعادلة 79g أيضًا عن طريق حساب مباشر؛ مصفوفة جاكوبي للتحويل في المعادلة 79e هي
(يجب ملاحظة ذلك في حساباته)).

بما أن δ s مُعبر عنها في المعادلة 79d بدلالة الكميتين العشوائيتين x + و x − ، فإن معرفة توزيعهما المشترك تُتيح حساب التوزيع الإحصائي P (δ) بتكامل P ( x + , x − ) على أحد المتغيرين عند قيمة ثابتة لـ δ. ونظرًا لتناظر الدالة في المعادلة 79g بالنسبة للمتغيرين x + و x − ، فإن P (δ) = P (1 − δ)، أي أن الدالة P (δ) متناظرة حول النقطة δ = 1/2. إذن
عند تقييم هذا التكامل (لـ 0 ≤ δ ≤ 1/2 ثم باستخدام التناظر المذكور أعلاه)، أخيرًا
| المعادلة 79ح |
القيمة المتوسطة= 1/2 بالفعل نتيجة لتناظر الدالة P (δ). وبالتالي، فإن متوسط قيمة السعة الأولية (في كل حقبة) لتذبذبات الدوال α، β، γ يزداد بمقدار Ω/2.
يتم إيجاد العلاقة الإحصائية بين الفترات الزمنية الكبيرة Ω وعدد الحقب s التي تحتويها من خلال التطبيق المتكرر للمعادلة 77 :
| المعادلة 80 |
مع ذلك، فإن حساب المتوسط المباشر لهذه المعادلة غير منطقي: نظرًا للتناقص البطيء للدالة W ( k ) في المعادلة 76 ، فإن القيم المتوسطة للكمية expξ ( s ) غير مستقرة بالمعنى المذكور أعلاه - إذ تزداد التقلبات بسرعة أكبر من القيمة المتوسطة نفسها مع ازدياد نطاق حساب المتوسط. ويمكن التخلص من هذا التذبذب بأخذ اللوغاريتم: الفترة الزمنية "اللوغاريتمية المزدوجة".
| المعادلة 81 |
يُعبَّر عنه بمجموع الكميات ξ ( p ) التي لها توزيع إحصائي مستقر. القيمة المتوسطة لـ τ هيلحسابلاحظ أنه يمكن إعادة كتابة المعادلة 77 على النحو التالي
| المعادلة 81أ |
بالنسبة للتوزيع الثابتوبفضل تناظر الدالة P (δ) أيضًا. لذلك
( w ( x ) من المعادلة 74 ). وبالتالي
| المعادلة 82 |
والذي يحدد متوسط الفترة الزمنية اللوغاريتمية المزدوجة التي تحتوي على s حقب متتالية.
بالنسبة لقيم s الكبيرة، يكون عدد الحدود في مجموع المعادلة 81 كبيرًا، ووفقًا للنظريات العامة لنظرية الإرجودية، تتوزع قيم τ حولوفقًا لقانون جاوس مع الكثافة
| المعادلة 82أ |
يُعد حساب التباين D τ أكثر تعقيدًا، لأنه لا يتطلب معرفة فقطوهناك حاجة إلى ذلك، ولكن أيضًا إلى الارتباطاتيمكن تبسيط الحساب بإعادة ترتيب الحدود في معادلة المجموع 81. باستخدام المعادلة 81أ، يمكن إعادة كتابة المجموع على النحو التالي:
لا يزداد الحدّان الأخيران مع ازدياد قيمة s ؛ ويمكن حذف هذين الحدّين لأن قوانين الحدّ للقيم الكبيرة s هي السائدة.
| المعادلة 82ب |
(يؤخذ التعبير 79d لـ δ p في الاعتبار). وبنفس الدقة (أي حتى الحدود التي لا تزداد مع s ) المساواة
| المعادلة 82 ج |
صحيح. في الواقع، بموجب المعادلة 79هـ
وبالتالي
بجمع هذه المتطابقة على p ، نحصل على المعادلة 82c . وأخيرًا، بنفس الدقةيتم تغيير x p تحت علامة الجمع، وبالتالي يتم تمثيل τ s على النحو التالي
| المعادلة 83 |
تباين هذا المجموع في حالة القيم الكبيرة لـ s هو
| المعادلة 84 |
يُؤخذ في الاعتبار أنه بفضل التجانس الإحصائي للتسلسل X، فإن الارتباطاتيعتمد فقط على الفروق | p − p ′ |. القيمة المتوسطةمتوسط المربع
مع الأخذ في الاعتبار أيضًا قيم الارتباطاتمع p = 1، 2، 3 (محسوبة عدديًا) يتم الحصول على النتيجة النهائية D τ s = (3.5 ± 0.1) s .
مع زيادة قيمة s، يزداد التذبذب النسبييؤول إلى الصفر عندما s −1/2 . بعبارة أخرى، تصبح العلاقة الإحصائية في المعادلة 82 شبه مؤكدة عند قيم s الكبيرة . وهذا يُتيح عكس العلاقة، أي تمثيلها على أنها اعتماد متوسط عدد الفترات s τ التي يتم تبادلها في فترة زمنية معينة τ من الزمن اللوغاريتمي المزدوج.
| المعادلة 85 |
إن التوزيع الإحصائي للقيم الدقيقة لـ s τ حول متوسطها هو أيضًا توزيع غاوسي بتباين σ
يتم إعطاء التوزيع الإحصائي المعني بواسطة نفس التوزيع الغاوسي الذي يكون فيه المتغير العشوائي الآن s τ عند قيمة τ معينة:
| المعادلة 86 |
من وجهة النظر هذه، فإن مصدر السلوك الإحصائي هو التعسف في اختيار نقطة البداية للفترة τ المتراكبة على التسلسل اللانهائي للحقب المتبادلة.
بالنسبة لكثافة المادة، يمكن إعادة كتابة المعادلة 79 مع مراعاة المعادلة 80 بالشكل التالي:
ثم، بالنسبة لتغير الطاقة الكلي خلال حقبة زمنية معينة ،
| المعادلة 87 |
يُشكّل الحدّ الذي يحتوي على المجموع في p المساهمة الرئيسية في هذا التعبير لأنه يتضمن أسًا ذا قوة كبيرة. وبترك هذا الحدّ فقط وأخذ متوسط المعادلة 87 ، نحصل في طرفها الأيمن على التعبير التاليوهذا يتطابق مع المعادلة 82 ؛ جميع الحدود الأخرى في المجموع (بما في ذلك الحدود التي تحتوي على η s في أسسها) تؤدي فقط إلى تصحيحات من رتبة نسبية 1/ s . لذلك،
| المعادلة 88 |
بفضل الطبيعة شبه المؤكدة للعلاقة بين τ s و s ، يمكن كتابة المعادلة 88 على النحو التالي
والذي يحدد قيمة اللوغاريتم المزدوج لزيادة الكثافة المتوسطة بواسطة فترات زمنية لوغاريتمية مزدوجة معينة τ أو بواسطة عدد معين من الحقب s .
توجد هذه العلاقات الإحصائية المستقرة تحديدًا لفترات زمنية لوغاريتمية مزدوجة ولزيادة الكثافة. أما بالنسبة لخصائص أخرى، مثل ln (ε ( s ) /ε (0) ) أو Ω (s) / Ω (0) = exp τ s، فإن التذبذب النسبي يزداد أُسّيًا مع زيادة نطاق المتوسط، مما يُفقد مصطلح القيمة المتوسطة معناه المستقر.
يمكن تتبع أصل العلاقة الإحصائية في المعادلة 88 من القانون الأولي الذي يحكم تغير الكثافة خلال حقب كاسنر الفردية. وفقًا للمعادلة 21 ، خلال التطور بأكمله لدينا
مع تغير قيمة 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 ) في نهاية حقبة زمنية. ويُعدّ معيار اختيار هذه الحالات هو المتباينة
| المعادلة 89 |
حيث | α ( 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) ، ستظهر الحالات الخطيرة في الأجزاء λ من هذه الفترة.
| المعادلة 90 |
(يُؤخذ المجموع الداخلي على جميع القيم k (1) ، k (2) ، ...، k ( s ) من 1 إلى ∞). من السهل إثبات أن هذه الحقبة تتقارب إلى القيمة λ1، الذي يُحدد مقداره بالحد الأول في المعادلة 90. ويمكن إثبات ذلك من خلال هيمنة قوية للحقبة التي يُستبدل بها |α ( s ) | = (s + 1) |α (0) |، بغض النظر عن أطوال الحقب k (1) ، k (2) ، ... (في الواقع، تزداد قيمة |α ( s ) | بشكل أسرع بكثير؛ حتى في أسوأ الحالات k (1) = k (2) = ... = 1، تزداد قيم |α ( s ) | مع q s |α (0) | حيث q > 1). مع ملاحظة أن
يحصل المرء
إذا كانت القيمة الابتدائية لـ x (0) تقع خارج المنطقة الخطرة λ، فلن تحدث أي حالات خطرة. أما إذا كانت تقع داخل هذه المنطقة، فستحدث حالات خطرة، ولكن بعد انتهائها، يستأنف النموذج تطوره "المنتظم" بقيمة ابتدائية جديدة قد تدخل المنطقة الخطرة من حين لآخر (باحتمالية λ). وتتكرر الحالات الخطرة باحتمالات λ₂ ، λ₃ ، ...، تتقارب تقاربًا مقاربًا إلى الصفر.
حل عام مع تذبذبات صغيرة
في النماذج المذكورة أعلاه، دُرست تطورات المقياس بالقرب من نقطة التفرد باستخدام مثال مقاييس الفضاء المتجانسة. ويتضح من خصائص هذا التطور أن البناء التحليلي للحل العام لنقطة تفرد من هذا النوع يجب أن يُجرى بشكل منفصل لكل مكون من مكونات التطور الأساسية: بالنسبة لحقب كاسنر، ولعملية الانتقالات بين الحقب الناتجة عن "الاضطرابات"، وللفترات الطويلة التي يحدث فيها اضطرابان في آن واحد. خلال حقبة كاسنر (أي عند الاضطرابات الصغيرة)، يُعطى المقياس بالمعادلة 7 دون الشرط λ = 0.
قام BKL بتطوير نموذج مستقل عن توزيع المادة (متجانس أو غير متجانس) لفترات طويلة مع تذبذبات صغيرة. وقد تبين أن التغير الزمني لهذا الحل مشابه جدًا للتغير في حالة النماذج المتجانسة؛ ويمكن الحصول على هذا الأخير من النموذج المستقل عن التوزيع عن طريق اختيار خاص للدوال العشوائية التي يحتويها. [ 52 ]
مع ذلك، من الملائم بناء الحل العام في نظام إحداثيات يختلف نوعًا ما عن الإطار المرجعي المتزامن: g₀α = 0 كما في الإطار المتزامن، ولكن بدلًا من g₀₀ = 1 ، أصبح g₀₀ = −g₃₃ . بتعريف موتر قياس الفضاء γαβ = −gαβ مرة أخرى ، نحصل بالتالي على
| المعادلة 91 |
يُكتب إحداثي الفضاء الخاص على الصورة x₃ = z، ويُكتب إحداثي الزمن على الصورة x₀ = ξ (باعتباره مختلفًا عن الزمن الخاص t ). سيُبين أن ξ يُقابل المتغير نفسه المُعرَّف في النماذج المتجانسة. يُشار إلى التفاضل بالنسبة إلى ξ و z ، على التوالي، بالنقطة والشرطة. تأخذ المؤشرات اللاتينية a و b و c القيم 1 و 2، والتي تُقابل إحداثيات الفضاء x₁ و x₂ ، والتي ستُكتب أيضًا على الصورة x و y . بالتالي، فإن المقياس هو
| المعادلة 92 |
يجب أن يحقق الحل المطلوب المتباينات
| المعادلة 93 |
| المعادلة 94 |
(تحدد هذه الشروط أن إحدى الدوال a 2 ، b 2 ، c 2 صغيرة مقارنة بالدالتين الأخريين، وهو ما كان عليه الحال أيضًا مع النماذج المتجانسة).
تعني المتباينة 94 أن المكونات γ a 3 صغيرة بمعنى أنه عند أي نسبة بين الإزاحتين dx a و dz ، يمكن حذف الحدود التي تحتوي على حاصل ضرب dx a dz في مربع عنصر الطول المكاني dl 2. لذلك، فإن التقريب الأول للحل هو مقياس معادلة 92 حيث γ a 3 = 0: [ ملاحظة 20 ]
| المعادلة 95 |
يمكن للمرء أن يقتنع بسهولة عن طريق حساب مكونات موتر ريتشي،،،باستخدام المعادلة المترية 95 والشرط الوارد في المعادلة 93، والذي ينص على أن جميع الحدود التي تحتوي على مشتقات بالنسبة للإحداثيات x a صغيرة مقارنةً بالحدود التي تحتوي على مشتقات بالنسبة للإحداثيات ξ و z (نسبتها ~ γ 33 / γ ab ). بعبارة أخرى، للحصول على معادلات التقريب الرئيسي، يجب اشتقاق γ 33 و γ ab في المعادلة 95 كما لو أنهما لا يعتمدان على x a .
| المعادلة 96 |
نحصل على المعادلات التالية: [ ملاحظة 21 ]
| المعادلة 97 |
| المعادلة 98 |
| المعادلة 99 |
يتم رفع وخفض المؤشر هنا بمساعدة γ ab . الكمياتوλ هما الانقباضاتووبذلك
| المعادلة 100 |
أما بالنسبة لمكونات موتر ريتشي،وبناءً على هذه الحسابات، فإن قيمتيهما تساويان صفرًا. وفي التقريب التالي (أي مع مراعاة صغر قيمة γ a 3 والمشتقات بالنسبة إلى x و y )، يتم تحديد الكميتين γ a 3 من خلال القيمتين γ 33 و γ ab المعروفتين مسبقًا .
يؤدي اختصار المعادلة 97 إلىوبالتالي،
| المعادلة 101 |
تختلف الحالات الممكنة باختلاف المتغير G. في الحالة المذكورة أعلاه ، g 00 = γ 33γ ab والحالة N > 0 (الكمية N متغيرة مع الزمن) تؤدي إلى حالات تفرد زمنية مهمة. بالتعويض في المعادلة 101، f₁ = ½ (ξ + z ) sin y و f₂ = ½ (ξ − z ) sin y ، نحصل على G من النوع
| المعادلة 102 |
لا يُقلل هذا الاختيار من عمومية الاستنتاجات؛ إذ يُمكن إثبات إمكانية العمومية (في التقريب الأول) بمجرد النظر إلى التحويلات المسموح بها المتبقية للمتغيرات. عند N < 0 (حيث N كمية مكانية)، يُمكن استبدال G = z، مما يُعمم مقياس أينشتاين-روزن المعروف . [ 53 ] عند N = 0، نصل إلى مقياس موجة روبنسون-بوندي [ 54 ] الذي يعتمد فقط على ξ + z أو فقط على ξ − z (انظر [ 55 ] ). تم وضع العامل sin y في المعادلة 102 للمقارنة المريحة مع النماذج المتجانسة. مع الأخذ في الاعتبار المعادلة 102 ، تُصبح المعادلات من 97 إلى 99 كما يلي:
| المعادلة 103 |
| المعادلة 104 |
| المعادلة 105 |
المعادلات الرئيسية هي المعادلة 103 التي تحدد مكونات γ ab ؛ ثم يتم إيجاد الدالة ψ من خلال تكامل بسيط للمعادلة 104 - المعادلة 105 .
تتراوح قيم المتغير ξ من 0 إلى ∞. يُنظر في حل المعادلة 103 عند حدين، ξ1 و1. عند قيم ξ الكبيرة، يمكن البحث عن حل يأخذ شكل تحليل 1 / √ ξ :
| المعادلة 106 |
وبذلك
| المعادلة 107 |
( تتطلب المعادلة 107 تحقق الشرط 102 ). بالتعويض بالمعادلة 103 في المعادلة 106 ، نحصل في الرتبة الأولى على
| المعادلة 108 |
حيث تشكل الكميات a و ac مصفوفة معكوسة للمصفوفة a و ac . حل المعادلة 108 له الشكل التالي :
| المعادلة 109 |
| المعادلة 110 |
حيث أن l a و m a و ρ هي دوال اختيارية للإحداثيات x و y مقيدة بالشرط eq. 110 المشتق من eq. 107 .
لإيجاد الحدود العليا لهذا التحليل، من الملائم كتابة مصفوفة الكميات المطلوبة γ ab بالشكل التالي:
| المعادلة 111 |
| المعادلة 112 |
حيث يرمز الرمز ~ إلى منقولة المصفوفة. المصفوفة H متناظرة وأثرها يساوي صفرًا. تضمن المعادلة 111 تناظر γ ab وتحقيق الشرط الوارد في المعادلة 102. إذا استُبدل exp H بالعدد 1، نحصل من المعادلة 111 على γ ab = ξ a ab ، حيث a ab مأخوذة من المعادلة 109. بعبارة أخرى، يتوافق الحد الأول من تحليل γ ab مع H = 0؛ أما الحدود الأعلى فتُحسب بتحليل القوى للمصفوفة H التي تُعتبر مكوناتها صغيرة.
تُكتب المكونات المستقلة للمصفوفة H على النحو σ و φ بحيث
| المعادلة 113 |
بإدخال المعادلة 111 في المعادلة 103 مع الإبقاء فقط على الحدود الخطية بواسطة H ، نحصل على σ و φ
| المعادلة 114 |
إذا حاولنا إيجاد حل لهذه المعادلات على شكل متسلسلة فورييه بدلالة الإحداثي z ، فسنحصل على معادلات بيسل كدوال لمعاملات المتسلسلة في ξ. أما الحدود التقاربية الرئيسية للحل عند قيم ξ الكبيرة فهي [ ملاحظة 22 ].
| المعادلة 115 |
المعاملان A و B دالتان عقديتان اختياريتان للإحداثيات x و y ، وتحققان الشروط اللازمة لقيمتي σ و φ الحقيقيتين؛ التردد الأساسي ω دالة حقيقية اختيارية للإحداثيات x و y . الآن، من المعادلتين 104 و 105، من السهل الحصول على الحد الأول من الدالة ψ:
| المعادلة 116 |
(هذا الحد يختفي إذا كانت ρ = 0؛ في هذه الحالة يكون الحد الرئيسي هو الحد الخطي لـ ξ من التفكيك: ψ = ξ q ( x , y ) حيث q دالة موجبة [ 56 ] ).
لذلك، عند قيم ξ الكبيرة، تتذبذب مركبات موتر القياس γ ab مع انخفاض ξ على خلفية انخفاض بطيء ناتج عن انخفاض عامل ξ في المعادلة 111. ينخفض المركب γ 33 = e ψ بسرعة وفقًا لقانون قريب من exp (ρ 2 ξ 2 )؛ وهذا ما يجعل الشرط في المعادلة 93 ممكنًا . [ ملاحظة 23 ]
بعد ذلك، ستنظر BKL في الحالة ξ1. يتم إيجاد التقريب الأول لحل المعادلة 103 من خلال الافتراض (الذي أكدته النتيجة) أنه في هذه المعادلات يمكن حذف الحدود التي تحتوي على مشتقات بالنسبة للإحداثيات:
| المعادلة 117 |
هذه المعادلة بالإضافة إلى الشرط الوارد في المعادلة 102 تعطي
| المعادلة 118 |
حيث أن λ a و μ a و s 1 و s 2 هي دوال اختيارية لجميع الإحداثيات الثلاثة x و y و z ، والتي ترتبط بشروط أخرى
| المعادلة 119 |
تعطينا المعادلتان 104 و 105 الآن
| المعادلة 120 |
المشتقاتتحتوي القيم المحسوبة باستخدام المعادلة 118 على الحدود ~ ξ 4 s 1 − 2 و ~ ξ 4 s 2 − 2 ، بينما الحدود المتبقية في المعادلة 117 هي ~ ξ −2 . لذلك، يُسمح بتطبيق المعادلة 103 بدلاً من المعادلة 117 في حالة الشرطين s 1 > 0 و s 2 > 0؛ ومن ثم 1 −> 0.
وبالتالي، عند قيم ξ الصغيرة، تتوقف تذبذبات الدالة γ ab بينما تبدأ الدالة γ 33 في الازدياد عند قيم ξ المتناقصة. هذا نمط كاسنر، وعند مقارنة γ 33 بـ γ ab ، لا ينطبق التقريب المذكور أعلاه.
وللتحقق من مدى توافق هذا التحليل، قامت شركة BKL بدراسة المعادلات= 0,= 0، وبحساب المكونات γ a 3 منها ، تأكدت صحة المتباينة 94. وقد أظهرت هذه الدراسة [ 52 ] أن المكونات γ a 3 في كلتا المنطقتين التقاربيتين كانت ≈ γ 33. لذا، فإن صحة المتباينة 93 تستلزم مباشرةً صحة المتباينة 94 .
يحتوي هذا الحل، كما ينبغي أن يكون عليه الحال في الحالة العامة لحقل في الفراغ، على أربع دوال اختيارية للإحداثيات المكانية الثلاثة x و y و z . في المنطقة ξ1 هذه الدوال هي، على سبيل المثال، lect 1 , lect 2 , μ 1 , s 1 . في المنطقة ξ1 يتم تعريف الدوال الأربع بواسطة متسلسلة فورييه بواسطة الإحداثي z من المعادلة 115 بمعاملات هي دوال لـ x و y ؛ على الرغم من أن تحليل متسلسلة فورييه (أو التكامل؟) يميز فئة خاصة من الدوال، إلا أن هذه الفئة كبيرة بما يكفي لتشمل أي مجموعة فرعية محدودة من مجموعة جميع الشروط الأولية الممكنة.
يحتوي الحل أيضًا على عدد من الدوال الاختيارية الأخرى للإحداثيات x و y . تظهر هذه الدوال الاختيارية ثنائية الأبعاد ، بشكل عام، لأن العلاقات بين الدوال ثلاثية الأبعاد في حلول معادلات أينشتاين تفاضلية (وليست جبرية)، بغض النظر عن المشكلة الأعمق المتعلقة بالمعنى الهندسي لهذه الدوال. لم يحسب BKL عدد الدوال ثنائية الأبعاد المستقلة لأنه في هذه الحالة يصعب استخلاص استنتاجات قاطعة نظرًا لأن الدوال ثلاثية الأبعاد تُعرَّف بواسطة مجموعة من الدوال ثنائية الأبعاد (انظر [ 52 ] لمزيد من التفاصيل). [ ملاحظة 24 ]
وأخيراً، يواصل BKL إظهار أن الحل العام يحتوي على الحل الخاص الذي تم الحصول عليه أعلاه للنماذج المتجانسة.
باستبدال متجهات الأساس للفضاء المتجانس من النوع التاسع لبيانكي في المعادلة 7، يأخذ مقياس الزمكان لهذا النموذج الشكل التالي
| المعادلة 121 |
عندما c 2في المعادلة a² و b² ، يمكن تجاهل c² في كل مكان باستثناء الحد c² dz² . للانتقال من الإطار المتزامن المستخدم في المعادلة 121 إلى إطار بشروط المعادلة 91 ، يتم إجراء التحويل dt = cd ξ/2 والتعويض z → z /2. بافتراض أن χ ≡ ln( a / b )1، نحصل من المعادلة 121 في التقريب الأول:
| المعادلة 122 |
وبالمثل، باستخدام متجهات الأساس للفضاء المتجانس من النوع الثامن لبيانكي، نحصل على
| المعادلة 123 |
وفقًا لتحليل الفضاءات المتجانسة أعلاه، في كلتا الحالتين ab = ξ (تبسيطًا)= ξ 0 ) و χ من المعادلة 51 ؛ يتم إعطاء الدالة c (ξ) بواسطة الصيغتين 53 و 61 ، على التوالي ، لنماذج النوعين التاسع والثامن.
يتم الحصول على مقياس مماثل للنوع الثامن من المعادلات 112 و 115 و 116 باختيار متجهين ثنائيي الأبعاد l و a و ma على النحو التالي :
| المعادلة 124 |
واستبدال
| المعادلة 125 |
للحصول على مقياس النوع التاسع، ينبغي استبدال
| المعادلة 126 |
(لحساب c (ξ) فإن التقريب في المعادلة 116 غير كافٍ ويتم حساب الحد في ψ الخطي بواسطة ξ [ 56 ] )
أُجري هذا التحليل للفراغ. إن تضمين المادة لا يُقلل من عمومية الحل ولا يُغير من خصائصه النوعية. [ 56 ] [ 52 ]
من أهم القيود التي تعيق الحل العام أن جميع الدوال ثلاثية الأبعاد الواردة في المعادلتين 122 و 123 يجب أن يكون لها فاصل تغير مميز واحد ومشترك. هذا الشرط وحده يسمح بتقريب جميع مشتقات المكونات المكانية للمقياس في معادلات أينشتاين بضرب هذه المكونات ببساطة في أعداد موجية مميزة، مما ينتج عنه معادلات تفاضلية عادية من النوع المُستنتج للنموذج المتجانس من النوع التاسع. وهذا هو سبب التطابق بين الحلول المتجانسة والعامة.
يستنتج من ذلك أن كلاً من نموذج النوع التاسع وتعميمه يحتويان على نمط تذبذبي ذي مقياس مكاني واحد ذي مقدار عشوائي، لا تحدده أي شروط فيزيائية. ومع ذلك، من المعروف أنه في الأنظمة غير الخطية ذات درجات الحرية اللانهائية، يكون هذا النمط غير مستقر ويتبدد جزئيًا إلى تذبذبات أصغر. في الحالة العامة للاضطرابات الصغيرة ذات الطيف العشوائي، سيكون هناك دائمًا بعض الاضطرابات التي ستزداد سعتها، مما يؤثر على طاقة العملية الكلية. ونتيجة لذلك، تظهر صورة معقدة لحركات متعددة المقاييس مع توزيع معين للطاقة وتبادل للطاقة بين التذبذبات ذات المقاييس المختلفة. ولا يقتصر هذا على الحالة التي يكون فيها تطور التذبذبات صغيرة النطاق مستحيلاً بسبب الشروط الفيزيائية. ففي هذه الحالة الأخيرة، يجب أن يوجد طول فيزيائي طبيعي يحدد الحد الأدنى للمقياس الذي تخرج عنده الطاقة من نظام ذي درجات حرية ديناميكية (كما يحدث، على سبيل المثال، في سائل ذي لزوجة معينة). ومع ذلك، لا يوجد مقياس فيزيائي متأصل لحقل الجاذبية في الفراغ، وبالتالي، لا يوجد عائق أمام تطور التذبذبات ذات المقاييس الصغيرة بشكل تعسفي. [ 57 ]
الاستنتاجات
تصف نظرية BKL حالات التفرد في الحل الكوني لمعادلات أينشتاين، والتي تتسم بطابع تذبذبي معقد. ورغم أن هذه الحالات قد دُرست بشكل أساسي على نماذج متجانسة مكانيًا، إلا أن هناك أسبابًا وجيهة للافتراض بأن حالات التفرد في الحل العام لمعادلات أينشتاين لها نفس الخصائص؛ وهذا ما يجعل نموذج BKL ذا أهمية بالغة في علم الكونيات.
يستند هذا البيان إلى حقيقة أن النمط التذبذبي في الاقتراب من التفرد ناتج عن اضطراب منفرد يُسبب أيضًا عدم استقرار في حل كاسنر المعمم. ويؤكد بناء النموذج التحليلي لفترات طويلة مع تذبذبات صغيرة عمومية هذا النموذج. ورغم أن هذا السلوك الأخير ليس عنصرًا ضروريًا لتطور المقياس بالقرب من التفرد، إلا أنه يمتلك جميع الخصائص النوعية الرئيسية: تذبذب المقياس في بُعدين مكانيين وتغير رتيب في البُعد الثالث مع اضطراب معين لهذا النمط في نهاية فترة زمنية محددة. مع ذلك، لم تُشرح الانتقالات بين حقب كاسنر في الحالة العامة للمقياس المكاني غير المتجانس بالتفصيل.
تم تأجيل دراسة المشكلة المتعلقة بالقيود المحتملة على هندسة الفضاء الناتجة عن التفرد. ومع ذلك، من الواضح منذ البداية أن نموذج BKL الأصلي قابل للتطبيق على كل من الفضاء المحدود وغير المحدود؛ ويتجلى ذلك في وجود نماذج تفرد تذبذبية لكل من الزمكان المغلق والمفتوح.
يُضفي النمط التذبذبي للاقتراب من التفرد بُعدًا جديدًا على مصطلح "محدودية الزمن". فبين أي لحظة زمنية محدودة في العالم t واللحظة t = 0، يوجد عدد لا نهائي من التذبذبات. وبهذا المعنى، تكتسب العملية طابعًا لانهائيًا. وبدلًا من الزمن t ، يُعدّ ln t متغيرًا أكثر ملاءمة لوصفها، حيث تُوسّع العملية لتشمل.
تُراعي نظرية BKL تطور المقياس في اتجاه تناقص الزمن. معادلات أينشتاين متناظرة بالنسبة لإشارة الزمن، لذا فإن تطور المقياس في اتجاه تزايد الزمن ممكنٌ بنفس القدر. مع ذلك، يختلف هذان الحالتان اختلافًا جوهريًا، لأن الماضي والمستقبل ليسا متكافئين بالمعنى الفيزيائي. لا يمكن أن يكون للتفرد المستقبلي معنى فيزيائي إلا إذا كان ممكنًا في ظل ظروف ابتدائية عشوائية موجودة في لحظة سابقة. لا يتطابق توزيع المادة والحقول في لحظة ما من تطور الكون بالضرورة مع الشروط المحددة اللازمة لوجود حل خاص معين لمعادلات أينشتاين.
يرتبط اختيار الحلول التي تتوافق مع الواقع بمتطلبات فيزيائية عميقة يستحيل إيجادها باستخدام نظرية النسبية الحالية وحدها، ويمكن التوصل إليها من خلال توليف النظريات الفيزيائية في المستقبل. لذا، قد يتضح أن هذا الاختيار يُبرز نوعًا خاصًا من التفرد (مثل التفرد المتناحي). ومع ذلك، من الطبيعي أكثر افتراض أن النمط التذبذبي، نظرًا لطبيعته العامة، هو السمة الرئيسية للمراحل التطورية الأولية.
في هذا السياق، تبرز أهمية خاصية نموذج "ميكس ماستر" الذي عرضه ميسنر [ 58 ] والمتعلقة بانتشار الإشارات الضوئية. في النموذج المتناحي، يوجد "أفق ضوئي"، أي أنه في كل لحظة زمنية، توجد مسافة قصوى يستحيل عندها تبادل الإشارات الضوئية، وبالتالي، وجود علاقة سببية: لا يمكن للإشارة أن تصل إلى هذه المسافات منذ اللحظة t = 0.
يتم تحديد انتشار الإشارة بواسطة المعادلة ds = 0. في النموذج المتناحي بالقرب من نقطة التفرد t = 0، يكون عنصر الفاصل الزمني هو، أينهي صيغة تفاضلية مكانية مستقلة عن الزمن. [ 59 ] بالاستبدالالعائد
| المعادلة 127 |
"المسافة"تصل الإشارة إلى
| المعادلة 128 |
بما أن η، مثل t ، تتغير قيمها بدءًا من 0، فإن الإشارات حتى "اللحظة" η لا يمكنها الانتشار إلا على مسافةوالذي يحدد أبعد مسافة إلى الأفق.
يُشكّل وجود أفق ضوئي في النموذج المتناحي إشكاليةً في فهم أصل التناحي الملحوظ حاليًا في الإشعاع المتبقي. فبحسب هذا النموذج، يعني التناحي الملحوظ خصائص متناحية للإشعاع القادم إلى الراصد من مناطق في الفضاء لا يمكن ربطها ببعضها البعض بشكل سببي. أما في نموذج التطور التذبذبي بالقرب من التفرد، فقد يختلف الوضع.
على سبيل المثال، في النموذج المتجانس للفضاء من النوع التاسع، تنتشر الإشارة في اتجاه تتغير فيه المقاييس لفترة طويلة وفقًا لقانون قريب من t . مربع عنصر المسافة في هذا الاتجاه هو dl² = t²، والعنصر المقابل في الفترة رباعية الأبعاد هوالاستبداليضع هذا في النموذج
| المعادلة 129 |
وبالنسبة لانتشار الإشارة، لدينا معادلة من النوع 128 مرة أخرى. والفرق المهم هو أن المتغير η يمر الآن بقيم تبدأ من(إذا كانت المعادلة المترية 129 صالحة لجميع قيم t بدءًا من t = 0).
لذلك، بالنسبة لكل "لحظة" معينة η، يتم إيجاد فترات وسيطة Δη كافية لتغطية الإشارة لكل مسافة محدودة.
بهذه الطريقة، خلال حقبة طويلة، ينفتح أفق ضوئي في اتجاه فضائي محدد. ورغم أن مدة كل حقبة طويلة تبقى محدودة، إلا أن الحقب تتغير خلال مسار تطور العالم عددًا لا نهائيًا من المرات في اتجاهات فضائية مختلفة. هذا الظرف يجعلنا نتوقع إمكانية وجود علاقة سببية بين الأحداث في الفضاء بأكمله في هذا النموذج. وبسبب هذه الخاصية، أطلق ميسنر على هذا النموذج اسم "عالم ميكس ماستر" نسبةً إلى اسم علامة تجارية لآلة خلط العجين.
مع مرور الوقت والابتعاد عن نقطة التفرد، يتزايد تأثير المادة على تطور الأبعاد، والذي كان ضئيلاً في المراحل الأولى من التطور، تدريجياً حتى يصبح هو المهيمن. ومن المتوقع أن يؤدي هذا التأثير إلى "تجانس" تدريجي للفضاء، ما يجعل خصائصه أقرب إلى نموذج فريدمان الذي يصف بدقة الحالة الراهنة للكون.
أخيرًا، يطرح BKL إشكالية جدوى اعتبار "حالة شاذة" لعالم ذي مادة لا نهائية الكثافة استنادًا إلى نظرية النسبية الحالية. ولا يمكن توضيح التطبيق الفيزيائي لمعادلات أينشتاين بصيغتها الحالية في ظل هذه الظروف إلا من خلال عملية توليف مستقبلية للنظريات الفيزيائية، ومن هذا المنطلق، لا يمكن حل هذه الإشكالية في الوقت الراهن.
من المهم ألا تفقد نظرية الجاذبية تماسكها المنطقي (أي ألا تؤدي إلى جدل داخلي) مهما كانت كثافة المادة. بعبارة أخرى، لا تتقيد هذه النظرية بالشروط التي تفرضها، والتي قد تجعل تطبيقها غير مقبول منطقيًا ومثيرًا للجدل عند الكثافات العالية جدًا؛ إذ لا يمكن أن تظهر القيود، من حيث المبدأ، إلا نتيجة لعوامل "خارجية" عن نظرية الجاذبية. هذا الظرف يجعل دراسة التفردات في النماذج الكونية مقبولة رسميًا وضرورية في إطار النظرية القائمة.
ملحوظات
- ↑ يمكن العثور على محاكاة متحركة مماثلة من تصميم ديفيد غارفينكل في [ 1 ] .
- ١ ٢ ٣ يستخدم نظام BKL نفس الاصطلاح الوارد في كتاب لاندو وليفشيتز (١٩٨٨) . تُشير المؤشرات اللاتينية إلى القيم ٠، ١، ٢، ٣؛ بينما تُشير المؤشرات اليونانية إلى القيم المكانية ١، ٢، ٣. يأخذ المقياس g ik الإشارة (+ − − −)؛ γ αβ = − g αβ هو موتر المقياس المكاني ثلاثي الأبعاد. يستخدم نظام BKL نظام وحدات تكون فيه سرعة الضوء وثابت أينشتاين للجاذبية مساويين لـ ١.
- ↑ يتم اشتقاق التعبير عن r عن طريق لوغاريتم معاملات القوة في المقياس: ln [ t 2 p α (1/ u ) ] = 2 p α (1/ u ) ln t .
- ↑ عندما ( p 1 , p 2 , p 3 ) = (0, 0, 1) فإن مقياس الزمكان المعادلة 1 مع dl 2 من المعادلة 2 يتحول إلى مقياس غاليلي مع الاستبدال t sh z = ζ, t ch z = τ، أي أن التفرد وهمي والزمكان مسطح.
- ↑ هنا وفيما يلي، يجب فهم جميع رموز عمليات المتجهات (الضرب المتجهي، وعمليات الدوران، والتدرج، وما إلى ذلك) بطريقة رسمية للغاية على أنها عمليات على المكونات المتغيرة للمتجهات l و m و n بحيث يتم تنفيذها في الإحداثيات الديكارتية x 1 و x 2 و x 3 .
- ↑ باستثناء الحالة ( p 1 , p 2 , p 3 ) = (0, 0, 1)، حيث تكون التفردية المترية وهمية.
- ↑ الثوابت λ و μ و ν هي ما يسمى بالثوابت الهيكلية لمجموعة الحركة الفضائية.
- ↑ تحتوي معادلات أينشتاين للفضاء المتجانس، في صيغتها الدقيقة، عمومًا، على ست دوال زمنية مختلفة γ ab ( t ) في المقياس. إن الحصول على نظام متسق من المعادلات الدقيقة في هذه الحالة للمقياس الذي يحتوي على ثلاث دوال زمنية فقط (γ 11 = а 2 ، γ 22 = b 2 ، γ 33 = c 2 ) يرتبط بتناظر يؤدي إلى اختفاء ستة مكونات لموتر ريتشي.
- ↑ يمكن إيجادالقيم التقاربية لـ ατ و βτ و γτ عند τ → −∞ دون حل المعادلة 29 بالكامل . يكفي ملاحظة أن المعادلة الأولى تأخذ شكل "جسيم" يتحرك في بُعد واحد في مجال جدار جهد أُسّي، حيث يلعب α دور الثابت. في هذا التشبيه، يشير نمط كاسنر إلى حركة حرة بسرعة ثابتة ατ = Λp / 1 . بعد الانعكاس عن الجدار، يتحرك الجسيم بحرية بسرعة ατ = −Λp / 1 . مع ملاحظة أنه من المعادلة... 29 α τ + β τ = const، و α τ + γ τ = const، يمكن للمرء أن يرى أن β τ و γ τ يأخذان القيم β τ = Λ( p 2 − 2 p 1 )، γ τ = Λ( p 3 − 2 p 1 ).
- ↑ إن إدخال المكونات غير القطرية لـ γ ab ( t ) يضفي بعض الميزات الجديدة على نموذج BKL: دوران المحاور المقابلة لقوى حقبة كاسنر؛ تمت دراسة هذه المشكلة في Belinsky و Khalatnikov و Lifshitz (1971).
- ↑ إن الثابت في وسيط الجيب، بالطبع، ليس بالضرورة هو نفسه ξ 0 في المعادلة 47 والمعادلة 48 ؛ ومع ذلك، فإن جعلهما متطابقين لا يغير بأي حال من الأحوال طبيعة الحل.
- ↑ في حساب أكثر دقة، يظهر مصطلح لوغاريتمي يتغير ببطء في وسيط الجيب، ويظهر مضاعف أمام الأس في التعبير عن с (ξ)، انظر Belinsky، Khalatnikov & Lifshitz 1970 ، الملحق ب .
- ↑ إذا استبدلنا sh 2χ بـ 2χ في المعادلة 49 وحلناها لجميع قيم ξ، نحصل على χ = c 1 J 0 (ξ) + c 2 N 0 (ξ) حيث J 0 و N 0 دالتان من دوال بيسل من النوع الأول والثاني على التوالي. هذا الحل يربط بين الحالتين الحدّيتين، ويسمح بربط المعاملات الثابتة في المعادلتين 52 و 55 بترتيب تقريبي.
- ↑ بما أن a و b و c لها بُعد الطول، فإن لوغاريتماتها تُعرَّف فقط حتى ثابت إضافي يعتمد على اختيار وحدات الطول؛ وبهذا المعنى ، فإن المعادلة 63 لها معنى مشروط يتوافق مع اختيار معين للقيمة الصفرية لـ α و β و γ.
- ↑ وفقًا للمعادلة 32 ، تكون الانتقالات كبيرة مع قيمة صغيرة لـ | p1 | ( أي قيمة كبيرة لـ u ) وتكون ≈1/| p1 | ~ u . ولكن حتى في هذه الحالة Δn ~ u n | αn |غير
- ↑ يُعدّ تحديد حدود الحقبة وفقًا للمعادلة 64 أمرًا ذا مغزى، لأنه في هذه الحالة، تشمل الحقبة جميع الفترات التي تتناقص فيها الدالة الثالثة، γ( t )، بشكل رتيب. إذا عُرّفت الحقبة بتسلسل قيم u من k + x إلى 1 + x ، فإن التناقص الرتيب لـ γ( t ) سيستمر خلال الفترة الأولى من الحقبة التالية.
- ↑ أطوال الحقب الزمنية كبيرة مقارنةً بالانتقالات بينها. وفقًا للمعادلة 33، تكون أطوال الانتقالات كبيرة عند قيم صغيرة لـ | p1 | (أي عند قيم كبيرة لـ u ) وتتناسب طرديًا مع 1/| p1 | ∝ u . ولكن حتى في هذه الحالة، فإن Δn ∝ u n | αn |يو إن .
- ↑ كانت المعادلة 74 معروفةً بالفعل لغوس ، وقد تناول روديون كوزمين معادلةً من النوع 73c في هذا السياق (انظر توزيع غاوس-كوزمين ). للمزيد من المعلومات حول السلوك الفوضوي والإنتروبيا للكسور المستمرة، انظر: ليناس فيبستاس، 2008، إنتروبيا الكسور المستمرة (إنتروبيا غاوس-كوزمين).
- ↑ إن رسم الدالة 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 .
- ↑ لاحظ أن هذا المقياس يسمح بتحويلات عشوائية من النوع ξ′ + z ″ = f 1 (ξ + z ), ξ′ − z ′ = f 2 (ξ − z ), x ′ a = f a ( x 1 , x 2 ).
- ↑ المعادلةهي نتيجة مباشرة للمعادلة 97 - المعادلة 99 إذاأوالقضيةلا يتطلب الأمر معالجة خاصة: يمكن إثبات أن مقياس الزمكان في هذه الحالة يتقارب (في التقريب الأول) إلى مقياس غاليليو.
- ↑ من الممكن البحث عن حل في صورة تكاملات فورييه؛ لم تُدرس هذه المسألة بالتفصيل. لذلك، لا تشترط نظرية BKL تحليل متسلسلة فورييه كشرط إلزامي لاعتماد الدالتين σ و φ على الإحداثيات.
- تؤدي الحدود المربعة H في المعادلة 103 إلى تصحيحات صغيرة فقط (≈1/ξ) في σ و φ. أما الحساب باستخدام الحدود التكعيبية فيؤدي إلى ظهور اعتماد ضعيف لـ A و B على ξ، والذي يمكن تمثيله بظهور أطوار لوغاريتمية في العوامل المتذبذبة في المعادلة 115. هذه الحسابات لحالة ρ = 0 واردة في Belinsky & Khalatnikov (1970 ، الملحق ب) (انظر الحالة المماثلة للنماذج المتجانسة، Belinsky, Khalatnikov & Lifshitz (1970 ، الملحق ب) ).
- ↑ يحتوي التفكيك المنتظم للحل العام لمعادلات أينشتاين (بالإضافة إلى الدوال الأربع ثلاثية الأبعاد) على ثلاث دوال مستقلة لإحداثيين (انظر بتروف 1969 ، الفصل 40 ؛ ليفشيتز وخالاتنيكوف (1963 ، الملحق أ) )
مراجع
- ↑ غارفينكل، ديفيد (2007). "عن التفردات وصناعة الخبز" . أينشتاين أونلاين . المجلد 3. معهد ماكس بلانك لفيزياء الجاذبية . 03-1014 . تاريخ الاسترجاع: 15 أكتوبر 2020 .
- 1 2 3 بيلينسكي، خالاتنيكوف وليفشيتز 1970
- ↑ ديماريه، هينو وسبينديل 1985 .
- ↑ ديماريه وآخرون 1986 .
- ^ ديماريت، دي روب وهينو 1989 .
- ↑ دامور وهينو 2000 .
- ↑ دامور وآخرون 2001 .
- ^ دامور وهينو ونيكولاي 2003 .
- ↑ Kac 1983 .
- ↑ دامور 2015 .
- ^ هينو وبيرسون وسبيندل 2008 .
- 1 2 3 4 5 6 ليفشيتز وخالاتنيكوف 1963
- 1 2 لاندو وليفشيتز 1988 ، الفصل 97
- ^ ليفشيتز وخالاتنيكوف 1961 أ .
- ^ ليفشيتز و خالاتنيكوف 1961 ب .
- 1 2 ليفشيتز، سوداكوف وخالاتنيكوف 1961
- ↑ هوكينج 1965 .
- ↑ هوكينج وإيليس 1968 .
- ↑ جيروش 1966 .
- 1 2 أشتيكار، هندرسون وسلون 2011
- ↑ بارو وتيبلر 1979 .
- ↑ بارو وتيبلر 1981 .
- 1 2 بيرغر 2002
- ↑ غارفينكل 2004 .
- ↑ بيرغر ومونكريف 1993 .
- ↑ بيرغر وآخرون 1998 .
- ↑ ويفر، إيزنبرغ وبيرغر 1998 .
- ^ أندرسون وريندال 2001 .
- ↑ دامور وآخرون 2002 .
- ↑ بيرغر ومونكريف 1998 .
- ↑ بيرغر ومونكريف 2000 .
- ↑ غارفينكل 2007 .
- ↑ ساوتومي، أخوري وغارفينكل 2010 .
- ↑ كاسنر 1921 .
- ↑ روغ 1994 .
- ^ بيني وشيروبيني وجانتزين 2007 .
- ↑ Landau & Lifshitz 1988 ، الفصل 117 ، المسألة 3.
- ↑ Landau & Lifshitz 1987 ، الفصل 134 ، المعادلة 134.15.
- ↑ ميسنر، ثورن وويلر 1973 ، ص 564 .
- ↑ نيلسون 1981 .
- ↑ بيلينسكي وخالاتنيكوف 1966 .
- ↑ خالاتنيكوف وليفشيتز 1970 .
- 1 2 بيلينسكي وخالاتنيكوف 1969أ
- 1 2 ليفشيتز وخالاتنيكوف 1970
- ↑ بيلينسكي، خالاتنيكوف وليفشيتز 1970 ، الملحق ج.
- ↑ ليفشيتز وخالاتنيكوف 1963 ، الملحق ج.
- ↑ تاوب 1951 .
- 1 2 3 4 ليفشيتز، ليفشيتز وخالاتنيكوف 1970
- 1 2 خالاتنيكوف وآخرون 1985
- ↑ تشيرنوف وبارو 1983 .
- ↑ بيلينسكي، خالاتنيكوف وليفشيتز 1970 ، الملحق أ.
- 1 2 3 4 بيلينسكي وخالاتنيكوف 1970
- ↑ أينشتاين وروزن 1937 .
- ^ بوندي وبيراني وروبنسون 1959 .
- ↑ Landau & Lifshitz 1988 ، الفصل 109.
- 1 2 3 بيلينسكي وخالاتنيكوف 1969ب
- ↑ بيلينسكي 1992 .
- ↑ ميسنر 1969 .
- ↑ Landau & Lifshitz 1988 ، الفصل 103–105.
فهرس
- أندرسون، لارس؛ ريندال، آلان د. (2001). "التفردات الكونية الهادئة". مجلة الاتصالات في الفيزياء الرياضية . 218 (3): 479-511 . arXiv : gr-qc/0001047 . Bibcode : 2001CMaPh.218..479A . doi : 10.1007/s002200100406 . S2CID 16167683 .
- أشتيكار، أبهاي ؛ هندرسون، آدم؛ سلون، ديفيد (2011). "صياغة هاميلتونية لتخمين بيلينسكي-خالاتنيكوف-ليفشيتز". مجلة Physical Review D. 83 ( 8) 084024. arXiv : 1102.3474 . Bibcode : 2011PhRvD..83h4024A . doi : 10.1103/PhysRevD.83.084024 . S2CID 119888235 .
- بارو، جون د .؛ تيبلر، فرانك ج. (1979). "تحليل دراسات التفرد العامة التي أجراها بيلينسكي، وخالاتنيكوف، وليفشيتز". تقارير الفيزياء . 56 (7): 371-402 . Bibcode : 1979PhR....56..371B . doi : 10.1016/0370-1573(79)90097-8 .
- بارو، جون د .؛ تيبلر، فرانك ج. (1981). "إعادة النظر في دراسات التفرد العامة". رسائل الفيزياء . 82أ (9): 441-445 . doi : 10.1016/B978-0-08-036364-6.50050-8 .
- بيلينسكي، فلاديمير أ .؛ خالاتنيكوف، إم (1966). "حل عام لمعادلات الجاذبية مع تفرد وهمي متزامن". مجلة الفيزياء التجريبية والنظرية . 49 (3): 1000. رمز Bibcode : 1966JETP...22..694B .
- بيلينسكي، فلاديمير أ .؛ خالاتنيكوف، إم (1969أ). "حول طبيعة التفردات في الحل العام لمعادلات الجاذبية". مجلة الفيزياء التجريبية والنظرية . 56 : 1700.
- بيلينسكي، فلاديمير أ .؛ خالاتنيكوف، إم (1969ب). "الحل العام لمعادلات الجاذبية مع وجود تفرد فيزيائي". مجلة الفيزياء التجريبية والنظرية . 57 : 2163.
- بيلينسكي، فلاديمير أ . خالاتنيكوف، إسحاق م . ليفشيتز، يفغيني م. (1970). "النظام الكلي مقترح لدرجة معينة في علم الكونيات النسبي" . Uspekhi Fizicheskikh Nauk . 102 (11): 463–500 . بيب كود : 1970UsFiN.102..463B . دوى : 10.3367/ufnr.0102.197011d.0463 .الترجمة الإنجليزية في: Belinskii, VA; Khalatnikov, IM; Lifshitz, EM (1970). "مقاربة تذبذبية لنقطة شاذة في علم الكونيات النسبي". Advances in Physics . 19 (80): 525–573 . Bibcode : 1970AdPhy..19..525B . doi : 10.1080/00018737000101171 .
- بيلينسكي، فلاديمير أ .؛ خالاتنيكوف، إسحاق م. (1970). "الحل العام لمعادلات الجاذبية مع تفرد تذبذبي فيزيائي". JETP . 59 : 314. Bibcode : 1971JETP...32..169B .
- بيلينسكي، فلاديمير أ .؛ خالاتنيكوف، إم. إم .؛ ليفشيتز، إي. إم. (1971). "النمط التذبذبي للاقتراب من نقطة تفرد في النماذج الكونية المتجانسة ذات المحاور الدوارة". مجلة الفيزياء التجريبية والنظرية . 60 (6): 1969-1979 .
- بيلينسكي، ف.أ. (1992). "اضطراب مجال الجاذبية بالقرب من نقطة تفرد كونية". رسائل JETP . 56 (9): 437-440 .
- بيلينسكي، فلاديمير ؛ هينو، مارك (2018). التفرد الكوني . مطبعة جامعة كامبريدج . arXiv : 1404.3864 . Bibcode : 2018cosm.book.....B . doi : 10.1017/9781107239333 . ISBN 978-1-107-04747-1.
- بيرغر، بيفرلي ؛ مونكريف، فنسنت (1993). "دراسة عددية للتفردات الكونية". مجلة Physical Review D. 48 ( 10): 4676–4687 . arXiv : gr-qc/9307032 . Bibcode : 1993PhRvD..48.4676B . doi : 10.1103 /PhysRevD.48.4676 . PMID 10016121. S2CID 450414 .
- بيرغر، بيفرلي ك .؛ مونكريف ، فنسنت (1998). "دليل عددي على أن التفرد في النماذج الكونية المتناظرة U (1) المستقطبة على T3 × R يهيمن عليه السرعة". مجلة Physical Review D. 57 ( 12): 7235–7240 . arXiv : gr-qc/9801078 . Bibcode : 1998PhRvD..57.7235B . doi : 10.1103/PhysRevD.57.7235 . S2CID 119539705 .
- بيرغر، بيفرلي ؛ غارفينكل، ديفيد؛ إيزنبرغ، جيمس ؛ مونكريف، فنسنت؛ ويفر، مارشا (1998). "التفرد في الانهيار الجذبي العام مكاني، محلي، وتذبذبي". رسائل الفيزياء الحديثة أ . 13 (19): 1565-1573 . arXiv : gr-qc/9805063 . Bibcode : 1998MPLA...13.1565B . doi : 10.1142/S0217732398001649 . S2CID 118889710 .
- بيرغر، بيفرلي ك .؛ مونكريف، فنسنت (2000). "بصمة ديناميكيات ميكس ماستر المحلية في نماذج الكون المتناظرة U (1)". مجلة Physical Review D ، 62، 123501. arXiv : gr-qc/0006071 . Bibcode : 2000PhRvD..62l3501B . doi : 10.1103/PhysRevD.62.123501 . S2CID 10389665 .
- بيرغر، بيفرلي ك. (2002). "المناهج العددية لدراسة تفردات الزمكان" . مراجعات حية في النسبية . 5 (1) 1. arXiv : gr -qc/0201056 . Bibcode : 2002LRR.....5....1B . doi : 10.12942/lrr-2002-1 . PMC 5256073. PMID 28179859 .
- بيني، دوناتو؛ تشيروبيني، مسيحي؛ جانتزن، روبرت (2007). “معاملة مؤشر ليفشيتز-خلاتنيكوف كاسنر وموتر فايل”. إل نوفو سيمينتو ب . 122 (5): 521– 536. أرخايف : 0710.4902 . بيب كود : 2007NCimB.122..521B . دوى : 10.1393/ncb/i2007-10396-4 . S2CID 15854390 .
- بيني، دوناتو؛ تشيروبيني، مسيحي؛ جيراليكو، أندريا؛ جانتزن، روبرت ت. (2009). “مخطط كهربية القلب في عالم Mixmaster”. الجاذبية الكلاسيكية والكمية . 26 (2) 025012. أرخايف : 0808.0828 . بيب كود : 2009CQGra..26b5012B . دوى : 10.1088/0264-9381/26/2/025012 . S2CID 14552653 .
- بوندي، هيرمان ؛ بيراني، إف إيه إي؛ روبنسون، آي. (1959). "موجات الجاذبية في النسبية العامة. الجزء الثالث: الموجات المستوية الدقيقة". وقائع الجمعية الملكية في لندن. السلسلة أ، العلوم الرياضية والفيزيائية . 251 (1267): 519-533 . Bibcode : 1959RSPSA.251..519B . doi : 10.1098/rspa.1959.0124 . JSTOR 100727. S2CID 122766998 .
- تشيرنوف، ديفيد ف.؛ بارو، جيه دي (1983). "الفوضى في عالم ميكس ماستر". رسائل المراجعة الفيزيائية . 50 (2): 134-137 . رمز Bibcode : 1983PhRvL..50..134C . doi : 10.1103/PhysRevLett.50.134 .
- دامور، تيبو ؛ هينو، مارك (2000). "الفوضى في علم الكونيات للأوتار الفائقة". رسائل المراجعة الفيزيائية . 85 (5): 920-923 . arXiv : hep-th/0003139 . Bibcode : 2000PhRvL..85..920D . doi : 10.1103/ PhysRevLett.85.920 . PMID 10991439. S2CID 10752273 .
- دامور، تيبو ؛ هينو، مارك ؛ جوليا، برنارد؛ نيكولاي، هيرمان (2001). "جبر كاك-مودي الزائدي والفوضى في نماذج كالوزا-كلاين". رسائل الفيزياء ب . 509 ( 3-4 ): 323-330 . arXiv : hep-th/0103094 . Bibcode : 2001PhLB..509..323D . doi : 10.1016/S0370-2693(01)00498-1 . S2CID 16964946 .
- دامور، تيبو ؛ هينو، مارك ؛ نيكولاي، هيرمان (2003). "البلياردو الكوني". الجاذبية الكلاسيكية والكمية . 20 (9): R145–R200. CiteSeerX 10.1.1.262.7733 . doi : 10.1088/0264-9381/20/9/201 . S2CID 250877925 .
- دامور، تيبو ؛ هينو، مارك ؛ ريندال، آلان د.؛ ويفر، مارشا (2002). "سلوك شبيه بكاسنر لأنظمة أينشتاين-المادة دون الحرجة". حوليات هنري بوانكاريه . 3 (6): 1049-1111 . arXiv : gr-qc/0202069 . Bibcode : 2002AnHP....3.1049D . doi : 10.1007/s000230200000 . S2CID 17092394 .
- ديماريت، جاك. هينو, مارك ; سبيندل ، فيليب (1985). “السلوك غير التذبذبي في علم الكونيات الفراغي كالوزا كلاين”. رسائل الفيزياء . 164ب ( 1– 3): 27– 29. بيب كود : 1985PhLB..164...27D . دوى : 10.1016/0370-2693(85)90024-3 .
- ديماريه، جاك؛ هانكين، جان لوك؛ هينو، مارك ؛ سبينديل، فيليب؛ تاورمينا، آنا (1986). "مصير سلوك ميكس ماستر في نماذج كالوزا-كلاين الكونية غير المتجانسة في الفراغ". رسائل الفيزياء ب . 175 (2): 129-132 . Bibcode : 1986PhLB..175..129D . doi : 10.1016/0370-2693(86)90701-X .
- ديماريه، جاك؛ دي روب، إيف؛ هينو، مارك (1989). "هل نماذج كالوزا-كلاين للكون فوضوية؟". المجلة الدولية للفيزياء النظرية . 28 (9): 1067-1079 . Bibcode : 1989IJTP...28.1067D . doi : 10.1007/BF00670349 . S2CID 121065247 .
- أينشتاين، ألبرت ؛ روزن، ناثان (1937). "حول الموجات الثقالية". مجلة معهد فرانكلين . 223 (1): 43-54 . Bibcode : 1937FrInJ.223...43E . doi : 10.1016/S0016-0032(37)90583-0 .
- فراسكا، ماركو (2006). "توسيع الاقتران القوي للنسبية العامة". المجلة الدولية للفيزياء الحديثة D. 15 ( 9): 1373-1386 . arXiv : hep-th/0508246 . Bibcode : 2006IJMPD..15.1373F . doi : 10.1142/S0218271806009091 . S2CID 15369171 .
- غارفينكل، ديفيد (2004). "محاكاة عددية للحالات الشاذة العامة". مجلة Physical Review Letters، المجلد 93 ، العدد 16، الصفحة 161101. arXiv : gr-qc/0312117 . Bibcode : 2004PhRvL..93p1101G . doi : 10.1103/PhysRevLett.93.161101 . PMID: 15524970. S2CID : 37917240 .
- غارفينكل، ديفيد (2007). "محاكاة عددية للحالات الشاذة الجاذبية العامة". الجاذبية الكلاسيكية والكمية . 24 (12): 295-306 . arXiv : 0808.0160 . Bibcode : 2007CQGra..24S.295G . doi : 10.1088/0264-9381/24/12/S19 . S2CID 16899122 .
- جيروش، روبرت ب. (1966). "التفردات في الأكوان المغلقة". رسائل المراجعة الفيزيائية . 17 (8): 445-447 . Bibcode : 1966PhRvL..17..445G . doi : 10.1103/PhysRevLett.17.445 .
- غودز، أندريه؛ بيدراك، ألكسندرا؛ بيتشوكي، Włodzimierz (2023). "ديناميات الكم المقابلة لسيناريو BKL الفوضوي" . المجلة الفيزيائية الأوروبية ج . 83 (2): 150. أرخايف : 2204.11274 . بيب كود : 2023EPJC...83..150G . دوى : 10.1140/epjc/s10052-023-11284-6 . S2CID 254044770 .
- هوكينغ، ستيفن و. (1965). "حدوث حالات شاذة في الأكوان المفتوحة" . رسائل المراجعة الفيزيائية . 15 (17): 689-690 . Bibcode : 1965PhRvL..15..689H . doi : 10.1103/PhysRevLett.15.689 .
- هوكينغ، ستيفن دبليو ؛ إليس، جي إف آر (1968). "إشعاع الجسم الأسود الكوني ووجود التفردات في كوننا". المجلة الفيزيائية الفلكية . 152 : 25. Bibcode : 1968ApJ...152...25H . doi : 10.1086/149520 .
- هاينزل، ج. مارك؛ أوغلا، كلايس؛ روهر، نيكلاس (2009). "جاذب البلياردو الكوني". مجلة الرياضيات الفيزيائية النظرية المتقدمة 13 ( 2): 293-407 . arXiv : gr-qc/0702141 . Bibcode : 2007gr.qc.....2141H . doi : 10.4310/ATMP.2009.v13.n2.a1 . S2CID 119470976 .
- هينو، مارك ؛ بيرسون، دانيال؛ سبينديل، فيليب (2008). "التفردات الفضائية والتناظرات الخفية للجاذبية" . مراجعات حية في النسبية . 11 (1): 1. arXiv : 0710.1818 . Bibcode : 2008LRR....11....1H . doi : 10.12942/ lrr -2008-1 . PMC 5255974. PMID 28179821 .
- كاك، فيكتور ج. (1983). جبر لي اللانهائي الأبعاد: مقدمة . سلسلة التقدم في الرياضيات. المجلد 44. بازل: بيركهاوزر. ص 245. doi : 10.1007/978-1-4757-1382-4 . ISBN 978-0-8176-3118-5.
- كاسنر، إدوارد (1921). "نظريات هندسية حول معادلات أينشتاين الكونية". المجلة الأمريكية للرياضيات . 43 (4): 217-221 . Bibcode : 1921AmJM...43..217K . doi : 10.2307/2370192 . JSTOR 2370192 .
- خالاتنيكوف، إسحاق م .؛ ليفشيتز، إي إم ؛ خانين، ك م؛ ششور، ل ن؛ سيناء، يا ج. (1985). "حول العشوائية في علم الكونيات النسبي". مجلة الفيزياء الإحصائية . 38 ( 1-2 ): 97-114 . Bibcode : 1985JSP....38...97K . doi : 10.1007/BF01017851 . S2CID 120521834 .
- خالاتنيكوف، إسحاق م . كامنشيك، ألكسندر يو. (2008). “ليف لانداو ومشكلة التفردات في علم الكونيات”. الفيزياء-Uspekhi . 51 (6): 609-616 . أرخايف : 0803.2684 . بيب كود : 2008PhyU...51..609K . دوى : 10.1070/PU2008v051n06ABEH006550 . S2CID 14878079 .
- خالاتنيكوف، آي إم ؛ ليفشيتز، إي إم (1970). "الحل الكوني العام لمعادلات الجاذبية مع وجود تفرد زمني". رسائل المراجعة الفيزيائية . 24 (2): 76. Bibcode : 1970PhRvL..24...76K . doi : 10.1103/PhysRevLett.24.76 .
- لاندو، ليف د .؛ ليفشيتز، يفغيني م. (1988). النظرية الكلاسيكية للحقول . سلسلة الفيزياء النظرية . المجلد 2 ( الطبعة السابعة). موسكو: ناوكا . ISBN 978-5-02-014420-0.
- لاندو، ليف د .؛ ليفشيتز، يفغيني م. (1987). ميكانيكا الموائع . سلسلة محاضرات في الفيزياء النظرية . المجلد 6 ( الطبعة الثانية). باتروورث-هاينمان . ISBN 978-0-08-033933-7.
- ليفشيتز، يفغيني م .؛ خالاتنيكوف، آي إم (1961). "حول حالات التفرد في الحلول الكونية لمعادلات الجاذبية. الجزء الأول". مجلة الفيزياء التجريبية والنظرية . 39 : 149.
- ليفشيتز، يفغيني م .؛ خالاتنيكوف، آي إم (1961). "حول حالات التفرد في الحلول الكونية لمعادلات الجاذبية. الجزء الثاني". مجلة الفيزياء التجريبية والنظرية . 39 : 800.
- ليفشيتز، يفغيني م .؛ سوداكوف، ف. ف.؛ خالاتنيكوف، إ. م. (1961). "حالات شاذة في الحلول الكونية لمعادلات الجاذبية. الجزء الثالث". مجلة الفيزياء التجريبية والنظرية . 40 : 1847.; Physical Review Letters , 6 , 311 (1961)
- ليفشيتز، إيفجيني م . خالاتنيكوف، إسحاق م. (1963). "مشاكل الكونيات النسبية" . Uspekhi Fizicheskikh Nauk . 80 (7): 391-438 . دوى : 10.3367 / UFNr.0080.196307d.0391 .الترجمة الإنجليزية في: Lifshitz, EM; Khalatnikov, IM (1963). "مشكلات في علم الكونيات النسبي". Advances in Physics . 12 (46): 185. Bibcode : 1963AdPhy..12..185L . doi : 10.1080/00018736300101283 .
- ليفشيتز، يفغيني م .؛ خالاتنيكوف، إسحاق م. (1970). "النمط التذبذبي للاقتراب من نقطة تفرد في النموذج الكوني المفتوح". رسائل JETP . 11 : 200-203 .
- ليفشيتز، إي إم ؛ ليفشيتز، آي إم ؛ خالاتنيكوف، آي إم (1970). "التحليل التقاربي للنمط التذبذبي للاقتراب من نقطة تفرد في النماذج الكونية المتجانسة". مجلة الفيزياء التجريبية والنظرية . 59 : 322. رمز Bibcode : 1970ZhETF..59..322L .
- ميسنر، تش. و. (1969). "عالم ميكس ماستر". رسائل المراجعة الفيزيائية . 22 (20): 1071-1074 . رمز Bibcode : 1969PhRvL..22.1071M . doi : 10.1103/PhysRevLett.22.1071 .
- ميسنر، تشارلز دبليو .؛ ثورن، كيب إس.؛ ويلر ، جون أرشيبالد (1973). الجاذبية . سان فرانسيسكو: دبليو إتش فريمان . ISBN 978-0-7167-0334-1.
- نيلسون، روبرت أ. (1981). "صياغة 3+1 لمعادلات النسبية العامة لسائل مثالي تعتمد على الإحداثيات". النسبية العامة والجاذبية . 13 (6): 569-580 . Bibcode : 1981GReGr..13..569N . doi : 10.1007/BF00757243 . S2CID 122196359 .
- بينروز، روجر (1965). "الانهيار الجاذبي وتفردات الزمكان" . رسائل المراجعة الفيزيائية . 14 (3): 57-59 . Bibcode : 1965PhRvL..14...57P . doi : 10.1103/PhysRevLett.14.57 .
- بيتروف، أليكسي ز. (1969). فضاءات أينشتاين . أكسفورد: دار بيرغامون للنشر . ISBN 978-0-08-012315-8.
- روغ، سفيند إي. (1994). "الفوضى في معادلات أينشتاين - توصيفها وأهميتها؟". الفوضى الحتمية في النسبية العامة . كاناناسكيس، ألبرتا، كندا: سبرينغر ساينس + بيزنس ميديا نيويورك. ص 359-422 . doi : 10.1007/978-1-4757-9993-4 . ISBN 978-1-4757-9995-8.
- ساوتومي، ريو؛ أخوري، راتيندراناث؛ غارفينكل، ديفيد (2010). "دراسة الانهيار الجذبي باستخدام حقول قياسية اختبارية". الجاذبية الكلاسيكية والكمية . 27 165019. arXiv : 1004.3569 . Bibcode : 2010CQGra..27p5019S . doi : 10.1088/0264-9381/27/16/165019 . S2CID 55229333 .
- تاوب، أ.هـ. (1951). "زمكانات فارغة تقبل مجموعة حركات بثلاثة معلمات". حوليات الرياضيات . السلسلة الثانية. 53 (3): 472-490 . Bibcode : 1951AnMat..53..472T . doi : 10.2307 / 1969567 . ISSN 0003-486X . JSTOR 1969567. MR 0041565 .
- ثورن، كيب (1994). الثقوب السوداء والتشوهات الزمنية: إرث أينشتاين المثير للجدل . دبليو دبليو نورتون . رقم ISBN 978-0-393-31276-8.]; كيب ثورن (2015). الديناميكا الهندسية: الديناميكا غير الخطية للزمكان المنحني (محاضرة فيديو). معهد فان لير القدس.
- أوغلا، كلايس؛ فان إلست، هينك؛ واينرايت، جون؛ إليس، جورج إف آر (2003). "الجاذب الماضي في علم الكونيات غير المتجانس". مجلة Physical Review D. 68 ( 10) 103502. arXiv : gr-qc/0304002 . Bibcode : 2003PhRvD..68j3502U . doi : 10.1103/PhysRevD.68.103502 . S2CID 53635765 .
- ويفر، مارشا؛ إيزنبرغ، جيمس ؛ بيرغر، بيفرلي ك. (1998). "سلوك ميكس ماستر في الزمكانات الكونية غير المتجانسة". رسائل المراجعة الفيزيائية . 80 (14): 2984-2987 . arXiv : gr-qc/9712055 . Bibcode : 1998PhRvL..80.2984W . doi : 10.1103/PhysRevLett.80.2984 . S2CID 119467400 .
- دامور، تيبو (19 نوفمبر 2015). بنية التفردات الكونية (محاضرة فيديو). معهد هنري بوانكاريه.
- النسبية العامة
- علم الكونيات الفيزيائي
- الحلول الدقيقة في النسبية العامة
- ثقب أسود
- التفردات الجاذبية
