Computational Methods of Wave Propagation for Semiclassical Models of High Harmonic Generation in Bulk Solids

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

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

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

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

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

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

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