تقع هذه الدراسة في موقع منعزل: أقرب الأبحاث المجاورة في الأرشيف تدرس الأوبئة عبر الاتصالات، والتفشيات الفيروسية، ومسار المرض في المدينة، لكن لا أحد تقريبًا يبني صورة مجهرية بوكلاء ذاتيين الحركة. تصف نماذج SIR وSI الكلاسيكية التطور الزمني للوباء بشكل جيد، لكنها صامتة حول كيفية تأثير حركة الناس وتنظيمهم المكاني على انتقال العدوى. أظهرت محاولات سابقة لإدخال نشاط الجسيمات في علم الأوبئة أن سرعة الوكلاء ونصف قطر العدوى يغيران منحنيات المرض، لكن دور الحماية الفردية، مثل التطعيم أو الكمامات، بقي غير مستكشف على المستوى المجهري.
حاكى الباحثون 3000 وكيل كجسيمات براونية نشطة في منطقة مربعة ثنائية الأبعاد بشروط حدية دورية. يتحرك كل جسيم ذاتيًا بسرعة ثابتة ويتعرض لانعطافات عشوائية، وتتنافر الجسيمات فيما بينها وفق جهد ويكس-تشاندلر-أندرسن (Weeks–Chandler–Andersen). ينقسم الوكلاء إلى قابلين للإصابة، ومصابين، ومحميين؛ يُصاب القابل إذا كان على مسافة أقل من نصف قطر العدوى من مصاب لمدة زمنية محددة. تراوحت نسبة الوكلاء المحميين بين 0 و50%، وتغيرت الكثافة السكانية بتغيير حجم المنطقة مع ثبات عدد الجسيمات. لكل مجموعة معايير، أجريت عشر محاكاة مستقلة على الحاسوبين العملاقين Miztli وLaboratorio de Supercómputo del Bajío، وقُربت منحنيات الإصابة بدالة الخطأ.
منحنيات عدد المصابين بمرور الزمن تأخذ شكلًا سينيًا، ويصل مشتقها، أي سرعة الإصابة، إلى أقصاه في اللحظة التي ينتشر فيها الوباء بأسرع ما يمكن. عندما تكون نسبة المحميين أقل من 30%، تنطبق منحنيات الإصابة المعيارية لمختلف الكثافات على منحنى عالمي واحد، وتنمو السرعة القصوى للإصابة كالجذر التربيعي للكثافة. يتوافق هذا مع تقدير بسيط: المسافة المميزة بين الجسيمات تتناسب عكسيًا مع الجذر التربيعي للكثافة، لذا فإن تواتر اللقاءات يتناسب مع الجذر التربيعي للكثافة. عندما تتجاوز نسبة المحميين 30%، تختفي العالمية: تتباعد المنحنيات، وتبدأ الحماية بالهيمنة. في التجمعات منخفضة الكثافة، يكون تثبيط الوباء أقوى منه في التجمعات الكثيفة، حيث يعوض التواتر العالي للقاءات جزئيًا تأثير الحماية.
يتميز هذا العمل عن الأعمال المجاورة بأنه يقدم لأول مرة بشكل منهجي وكلاء محميين في نموذج مجهري للمادة النشطة ويكشف عن تحول واضح يفصل بين نمطين. تدرس الأبحاث المجاورة إما ديناميكيات الاتصال البحتة دون حركة ذاتية، أو التفشيات كصدمات عشوائية، لكنها لا تُظهر كيف تتفاعل الحماية الفردية مع الكثافة على مستوى حركة الجسيمات. هنا اكتُشف قانون عالمي للحماية المنخفضة وانتهاكه عند الحماية العالية، مما يعطي تفسيرًا مجهريًا للعلاقة غير الخطية بين معدلات الاتصال والكثافة، والتي لوحظت سابقًا في النماذج المكانية.
يمكن أن يتطور الموضوع نحو تحديد أدق للتحول: يخطط الباحثون لدراسة قيم إضافية لنسبة المحميين حول 30% لفهم ما إذا كان الانتقال ظاهرة حرجة حادة أم تحولًا تدريجيًا. كما أن مقارنة النموذج بشبكات التخلل الوبائي (epidemic percolation networks) واعدة لربط التنظيم المكاني للوكلاء بظهور العدوى واسعة النطاق وتحديد الأصل المجهري للقانون الجذري.
قد تؤثر النتائج على نمذجة الأوبئة في البيئات الحضرية، حيث تختلف الكثافة السكانية ومستوى التطعيم أو ارتداء الكمامات من حي إلى آخر. كما أنها مفيدة لفهم التفشيات الحيوانية المنشأ في الأماكن ذات الكثافة الحيوانية العالية، مثل مزارع الفراء.
الخطوة التالية هي دراسة منطقة التحول حول 30% من الحماية بالتفصيل، بتغيير نسبة المحميين بخطوات صغيرة. ثم تُخطط مقارنة النموذج بشبكات التخلل الوبائي لإيجاد علاقة كمية بين البنية المكانية وعتبة الوباء.
يتناول العمل مشكلة قديمة في علم الأوبئة: كيف ترتبط معدلات الاتصال بالكثافة السكانية بشكل غير خطي، وكيف تغير إجراءات الحماية الفردية هذه العلاقة. تفترض نماذج الحجرات التقليدية اختلاطًا متجانسًا، وهو ما لا يعكس التغاير المكاني الحقيقي؛ تُظهر هذه الدراسة أن أخذ الحركة والتفاعلات المحلية في الاعتبار يؤدي إلى أنماط جديدة لا يمكن استنتاجها من معادلات SIR.
في النموذج، عندما تكون نسبة المحميين 30% بالضبط، يتأخر زمن بلوغ ذروة الوباء من 138 إلى 161 وحدة زمنية غير بعدية، أي أن الحماية لا تقلل ارتفاع الذروة فحسب، بل تؤخرها بشكل ملحوظ.