Pressure-Induced Structural Phase Transitions in MBANP: A First-Principles Study of Elastic, Phonon, and Optical Properties

ملخص تقني: التحولات الطورية الهيكلية المستحثة بالضغط في مادة MBANP

بيان المشكلة
تعد بلورات المواد العضوية ذات الخصائص البصرية غير الخطية (NLO)، وتحديداً مركب 2-(α-methylbenzylamino)-5-nitropyridine (MBANP)، بالغة الأهمية لتقنيات الضوئيات والإلكترونيات الضوئية من الجيل التالي نظراً لكفاءتها العالية في توليد التوافق الثاني (SHG) ومرونتها الهيكلية. وبينما تم توصيف MBANP بشكل موسع فيما يتعلق بسلوكه المعتمد على درجة الحرارة والتفاعلات بين الجزيئات، فإن استجابته للضغط الهيدروستاتيكي لا تزال غير مستكشفة تماماً. إن فهم هذه الاستجابة أمر جوهري لتوضيح كيف تتحكم التفاعلات غير التساهمية في استقرار البلورات ونشاطها البصري غير الخطي، وهو أمر ضروري عملياً لتقييم الاستقرار الميكانيكي لهذه المواد في الأجهزة المصغرة المعرضة للإجهاد.

المنهجية
تعتمد هذه الدراسة على نهج قائم على المبادئ الأولى باستخدام نظرية الدالة الوظيفية للكثافة (DFT) عبر كود CASTEP. تستخدم الحسابات تقريب التدرج المعمم (GGA) مع دالة (Perdew–Burke–Ernzerhof (PBE، معززة بتصحيح (Tkatchenko–Scheffler (TS لنمذجة التفاعلات الضعيفة بين الجزيئات بدقة. كما استُخدمت جهود شبه جهد (pseudopotentials) محافظة على النواة ومجموعة موجات مستوية بطاقة قطع قدرها 830 إلكترون فولت. وتتقصى الدراسة بشكل منهجي بلورات MBANP الأحادية تحت ضغط هيدروستاتيكي يتراوح من 0 إلى 20 جيجا باسكال. وتشمل التحليلات الرئيسية ما يلي:

• التطور الهيكلي: مراقبة معاملات الشبكة البلورية والتشكل الجزيئي.
• الخصائص المرنة: حساب الثوابت المرنة ($C_{ij}$) واشتقاق المعاملات الميكانيكية (الكتلة، القص، يونج، ونسبة بواسون) عبر متوسط (Voigt-Reuss-Hill)، جنباً إلى جنب مع التحقق من معايير بور للاستقرار.
• البنية الإلكترونية: تحليل بنى النطاقات، وكثافة الحالات (PDOS)، وكثافة الشحنة، ودالة التمركز الإلكتروني (ELF)، وآليات انتقال الشحنة.
• التفاعلات بين الجزيئات: استخدام تحليل سطح هيرشفيلد (Hirshfeld surface) ومخططات بصمات ثنائية الأبعاد لتحديد أنواع التلامس.
• الخصائص الاهتزازية والبصرية: محاكاة أطياف رامان وحساب المقاييس البصرية المعتمدة على التردد (الدالة العازلة، معامل الانكسار، الامتصاص، والانعكاسية).

المساهمات والنتائج الرئيسية

1. التحول الطوري الهيكلي عند ~14 جيجا باسكال:
تحدد الدراسة تحولاً هيكلياً متماثل الشكل (isosymmetric) مستحثاً بالضغط بالقرب من 14 جيجا باسكال. تحت هذا الحد، تظهر البلورة انضغاطية متباينة الخواص بشكل واضح (المحور $a$ > المحور $c$ > المحور $b$). وعند 14 جيجا باسكال، يصاحب التقلص غير الطبيعي في المحور $a$ توسع منفصل في المحورين $b$ و $c$ وقفزة في الزاوية $\beta$. يشير هذا إلى إعادة ترتيب جزيئي منسق أو تشوه شبيه بالقص، مما ينتج عنه طور عالي الضغط بخصائص ميكانيكية شبه متماثلة المناحي (quasi-isotropic). وبالرغم من هذه التغييرات، يتم الحفاظ على تناظر المجموعة الفراغية ($P2_1$).

2. الاستقرار الميكانيكي والتقوية:
تظل مادة MBANP مستقرة ميكانيكياً وديناميكياً طوال نطاق 0–20 جيجا باسكال، وهو ما أكده استيفاء معايير بور المعدلة وغياب الترددات التخيلية في أطياف تشتت الفونون. تمر المادة بعملية تقوية تدريجية، حيث يرتفع معامل الكتلة (Bulk Modulus) من 10.7 جيجا باسكال إلى 105.2 جيجا باسكال، ومعامل القص (Shear Modulus) من 3.76 جيجا باسكال إلى 27.2 جيجا باسكال. تزدوت نسبة بواسون ($\nu$) بشكل رتيب من 0.343 إلى 0.4 في نطاق 0–12 جيجا باسكال قبل أن تستقر حول 0.38، مما يشير إلى طبيعة مطيلية (ductile). تظهر نسبة $B/G$ تبايناً كبيراً، حيث ترتفع من 2.84 إلى 4.7 بين 0 و 12 جيجا باسكال، ثم تنخفض إلى 3.86 عند 20 جيجا باسكال، مع وجود نقطة انقلاب عند 14 جيجا باسكال ترتبط بالتحول الهيكلي.

3. تعزيز البنية الإلكترونية وانتقال الشحنة:
يعزز الضغط بشكل كبير انتقال الشحنة داخل الجزيء (CT) على طول محور المانح-المستقبل (من مجموعة methylbenzylamino إلى مجموعة nitro). تتقلص فجوة النطاق (bandgap) من 2.226 إلكترون فولت عند 0 جيجا باسكال إلى 1.056 إلكترون فولت عند 20 جيجا باسكال، مدفوعة باتساع نطاق التكافؤ (مدارات $\pi$ للبيريدين) ونطاق التوصيل (حالات $\pi^*$ للنيترو). يكشف تحليل كثافة الشحنة عن زيادة في نضوب الإلكترونات عند نيتروجين المانح وتراكمها عند أكسجين المستقبل، إلى جانب قصر طول رابطة N–H⋯O الهيدروجينية، مما يؤكد تعزيز تأثير "الدفع والسحب" (push-pull).

4. إعادة تنظيم التفاعلات بين الجزيئات:
يحدد تحليل سطح هيرشفيلد تحولاً في هيمنة التفاعلات مدفوعاً بالضغط. تنخفض المساهمة في تلامسات فان در فال الضعلة والمتماثلة المناحي (H⋯H و C⋯H)، بينما تزداد التفاعلات القوية والموجهة. ومن الجدير بالذكر أن تلامسات C⋯C (التي تشير إلى تكدس $\pi$–$\pi$) تتزايد بشدة، وتتحول تلامسات C–H⋯O المحددة إلى روابط هيدروجينية حقيقية (مسافات < 2.6 أنجستروم) فوق 14 جيجا باسكال. ويُعد هذا إعادة التنظيم آلية رئيسية للحفاظ على الاستقرار تحت الضغط الشديد.

5. الاستجابات الاهتزازية والبصرية:
تظهر أطياف رامان إزاحة زرقاء (blue-shift) موحدة عبر جميع الأنماط، تُعزى إلى انضغاط الروابط وزيادة ثوابت القوة. ويحدث انقطاع متميز في شدة القمم وتوسعها عند 14 جيجا باسكال، مما يؤكد التحول الطوري. من الناحية البصرية، يزداد ثابت العزل الساكن ومعامل الانكسار مع الضغط (يرتفع من ~3.56 إلى ~4.4 ومن ~1.86 إلى ~2.1، على التوالي). ومع ذلك، تزاح حافة الامتصاص وقمم الدالة العازلة الرئيسية نحو طاقات أعلى (إزاحة زرقاء)، وهو ما يتوافق مع الاستنتاج بأن حافة الامتصاص تنتقل إلى طاقة أعلى. وبناءً على ذلك، يضعف امتصاص المادة في المنطقة المرئية، مما يؤدي إلى انخفاض في النفاذية (من ~87% إلى ~73% لعينة بسمك 100 نانومتر) وزيادة في الانعكاسية.

الأهمية
يوفر هذا العمل أول تحقيق نظري شامل لمادة MBANP تحت الضغط العالي، مما يسد فجوة حرجة في فهم العلاقات بين البنية والخواص. وتثبت الدراسة أن الضغط الهيدروستاتيكي يعمل كأداة قوية لضبط الخصائص الوظيفية لبلورات الـ NLO العضوية دون تغيير تناظرها. ومن خلال توضيح آليات التقوية الميكانيكية، وتعزيز انتقال الشحنة، وإعادة التنظيم بين الجزيئات، تؤكد النتائج استقرار MBANP للتطبيقات في الأجهزة الإلكترونية الضوئية المصغرة، وتقدم توجيهات لتصميم الجيل القادم من البلورات الجزيئية ذات الخصائص الكهروميكانيكية والبصرية القابلة للضبط.