High-Harmonic Generation Probes Quantum Nuclear Fluctuations in Liquid Water
تُظهر هذه الدراسة أن مطيافية توليد التوافقيات العالية، عند دمجها مع ديناميكيات الجزيئات ذات التكامل المساري، يمكنها رصد وتوصيف التقلبات النووية الكمومية في الماء السائل من خلال الكشف عن كيفية قيام عدم التمركز النووي الكمومي بتثبيط عوائد التوافقية وإحداث حد أدنى من التداخل عبر اتساع نطاق حزمة التكافؤ الناجم عن الروابط الهيدروجينية.
المؤلفون الأصليون:George Trenins, Nicolas Tancogne-Dejean, Angana Mondal, Angel Rubio, Hans Jakob Wörner, Aaron Kelly, Ofer Neufeld, Mariana Rossi
يُعد الماء السائل الأكثر شيوعاً على كوكب الأرض، ومع ذلك لا يزال واحداً من أكثر المواد غموضاً. إن قدرته الفريدة على إذابة الكثير من المواد، وسلوكه الغريب عند التجمد، ودوره في كل خلية حية، كلها تنبع من الطريقة التي تترابط بها جزيئاته مع بعضها البعض. فهي تشكل شبكة واسعة ومتحركة ممسوكة بروابط هيدروجينية، تعمل مثل شرائط "الفيلكرو" (اللاصق) الصغيرة والمرنة التي تربط جزيء الماء بالآخر. لعقود من الزمن، عرف العلماء أنه لفهم هذه الشبكة حقاً، لا يمكنهم معاملة الذرات داخل جزيئات الماء ككرات صلبة بسيطة تتبع القواعد المتوقعة للفيزياء الكلاسيكية. بدلاً من ذلك، تتصرف أخف الذرات — وهي نويات الهيدروجين — كأنها سحب احتمالية ضبابية، تهتز بطاقة لا تتوقف أبداً، حتى عند الصفر المطلق. هذا الحركة المستمرة والمضطربة، المعروفة باسم "اللا-تموضع الكمي"، تعيد تشكيل الروابط الهيدروجينية ببراعة، مما يجعلها أقوى أو أضعف قليلاً اعتماداً على كيفية تذبذب الذرات. وبينما هذا التأثير راسخ من الناحية النظرية، فإن رؤيته بشكل مباشر في سائل، حيث تصطدم الجزيئات ببعضها البعض باستمرار، كان تحدياً استعصى على التجريبيين.
لقد اتخذ فريق من الباحثين الآن خطوة كبيرة نحو تصور هذه "التذبذبات الكمية" غير المرئية باستخدام نوع جديد وقوي من الضوء لسبر أغوار الماء السائل. لم يستخدموا مجهراً، بل تقنية تسمى "توليد التوافقيات العالية"، والتي تتضمن قذف عينة بنبضة ليزر شديدة الكثافة. عندما يصطدم الليزر بالماء، فإنه ينتزع الإلكترونات من ذراتها ثم يعيد صدمها، مما يتسبب في انبعاث ومضة ضوئية من الذرات بترددات أعلى بكما في الليزر نفسه. يحمل هذا الضوء المنبعث بصمة تفصيلية للبنية الإلكترونية وترتيب الذرات في لحظة الاصطدام. ومن خلال محاكاة هذه العملية باستخدام الحواسيب الفائبة، قارن الباحثون بين نسختين من الماء السائل: نسخة حيث عوملت ذرات الهيدروجين كجسيمات كلاسيكية صلبة، ونسخة أخرى حيث عوملت كسحب كمية. وكشفت المقارنة بين المحاكاتين أن الطبيعة الكمية للنويات تترك علامة مميزة وقابلة للقياس على الضوء المنبعث من الماء.
قام الباحثون ببناء نموذج رقمي للماء السائل في درجة حرحة الغرفة، حيث أنشأوا 50 لقطة لحركة الجزيئات للحالة الكلاسيكية، و50 لقطة أخرى للحالة الكمية. بالنسبة للقطات الكلاسيكية، استخدموا الفيزياء القياسية لحساب كيفية تحرك الذرات. أما بالنسبة للقطات الكمية، فقد استخدموا طريقة أكثر تعقيداً تسمى "ديناميكا الجزيئات ذات المسار المتكامل"، والتي تأخذ في الاعتبار حقيقة أن ذرات الهيدروجين ليست نقاطاً بل سحباً منتشرة. ثم أخضعوا كل لقطة لنبضة ليزر محاكية، وحسبوا كيف ستستجيب الإلكترونات ونوع الضوء الذي سيُنبعث. وعندما قارنوا النتائج، وجدوا فرقاً هائلاً في سطوع الضوء المنبعث عند ترددات محددة. في المحاكاة التي عوملت فيها ذرات الهيدروجين كسحب كمية، تم كبح الضوء المنبعث عند تردد محدد — وهو التوافق التاسع عشر لليزر — بمقدار يصل إلى عشرة أضعاف مقارنة بالمحاكاة الكلاسيكية. لم يكن هذا الكبح مجرد تقلب عشوائي؛ بل ظهر كفجوة حادة في الطيف، وهي إشارة واضحة على أن الطبيعة الكمية للنويات تداخلت مع عملية انبعاث الضوء.
لفهم سبب حدوث هذه الفجوة، فحص الفريق بدقة بنية الماء في عمليات المحاكاة الخاصة بهم. فحصوا المسافات بين الذرات وزوايا الجزيئات، بحثاً عن تغيير هيكلي بسيط يمكن أن يفسر هذا الانخفاض في شدة الضوء. ووجدوا أنه بينما أدى التعامل الكمي إلى توسيع توزيع المسافات بين ذرات الهيدروجين والأكسجين داخل الجزيء الواحد، إلا أن هذا التغيير الداخلي وحده لم يفسر التأثير. كما اختبروا ما إذا كان تغيير كثافة الماء سيؤدي إلى إزاحة موضع قطع الهضبة الأولى، وقد فعل ذلك بالفعل، لكن السبب الرئيسي للكبح كان أكثر دقة. اكتشف الباحثون أن "التذبذب الكمي" لذرات الهيدروجين تسبب في جعل مستويات الطاقة الإلكترونية في الماء تتوسع وتتشتت. في المحاكاة الكلاسيكية، كانت الإلكترونات تشغل حالات طاقة محددة وحادة، أما في المحاكاة الكمية، فقد أدت الحركة المستمرة للنويات إلى طمس هذه الحالات، مما خلق نطاقاً أوسع من الطاقات. هذا الطمس تسبب في إلغاء موجات الضوء المنبعثة من أجزاء مختلفة من الماء لبعضها البعض عند ذلك التردد المحدد، مما خلق الفجوة الملحوظة.
أكد الفريق هذا التفسير من خلال بناء نموذج رياضي مبسط للماء، حيث جردوا النموذج من التفاصيل المعقدة للتركيز فقط على الفيزياء الأساسية لكيفية تحرك الإلكترونات بين مستويات الطاقة. ومن خلال ضبط معايير هذا النموذج ليتوافق مع المحاكاة الكمية والكلاسيكية، أظهروا أن اتساع الحالات الإلكترونية كان بالفعل هو السبب المباشر لفجوة التداخل. ويعد هذا الاكتشاف مهماً لأنه يثبت أن تقنية "توليد التوافقيات العالية" حساسة بما يكفي لرصد الطبيعة الميكانيكية الكمية للنويات الذرية في السائل. كما يشير إلى إمكانية استخدام هذه التقنية في المستقبل لمراقبة التأثيرات الكمية في الوقت الفعلي، مما قد يسمح للعلماء بمراقبة كيفية تكسر الروابط الهيدروجية وإعادة تشكلها أثناء التفاعلات الكيميائية أو كيفية انتقال البروتونات عبر الماء. يوضح هذا العمل أنه لفهم سلوك السوائل مثل الماء بشكل كامل، يجب مراعاة "الضبابية الكمية" للذرات، حيث أن هذه الحركة الدقيقة تغير بشكل جوهري كيفية تفاعل المادة مع الضوء.
ملخص تقني: مجسات توليد التوافقيات العالية تستقصي التقلبات النووية الكمومية في الماء السائل
بيان المشكلة يمتلك الماء السائل خصائص فريدة تحكمها شبكة الروابط الهيدروجينية، وهي شبكة حساسة للغاية لما إذا كانت النوى الذرية تُعامل كجسيمات كلاسيكية أو كيانات كمومية. ومن المعروف أن التأثيرات النووية الكمومية (NQEs)، مثل الحركة النقطية الصفرية والانتشار، تؤثر على الخصائص الهيكلية (مثل الكثافة والانتشار) والخصصاف الإلكترونية (مثل تقليل فجوة النطاق وتوسيع حالات كثافة الحالات). ومع ذلك، يظل اكتشاف البصمات الطيفية لهذه التأثيرات النووية الكمومية في الطور السائل على مقاييس زمنية بالفيمتو ثانية تحديًا كبيرًا. وبينما برز توليد التوافقيات العالية (HHG) كأداة قوية لاستقصاء الديناميكيات الإلكترونية فائقة السرعة في الغازات والمواد الصلبة، فإن تطبيقه على السوائل لا يزال في مراحله الأولى. ويظل هناك سؤال جوهري مفتوح حول ما إذا كان بإمكان توليد التوافقيات العالية (HHG) رصد بصمات الحركة النووية الكمومية في السوائل الضخمة، وإذا كان الأمر كذلك، كيف تظهر هذه التأثيرات في الأطياف المنبعثة. تناولت الدراسات النظرية السابقة التأثيرات النووية الكمومية في الحالة الغازية والأنظمة البلورية، ولكن لم يسبق لأي عمل نظري أن تناول السوائل الضخمة حيث تُظهر التأثيرات النووية الكمومية وجودها حتى في نظام الاستجابة الخطية.
المنهجية استخدم المؤلفون نهجًا حاسوبيًا متعدد المقاييس يجمع بين ديناميكا المسار التكاملية (PIMD) ونظرية وظيفية الكثافة المعتمدة على الزمن (TDDFT) في الوقت الفعلي لمحاكاة توليد التوافقيات العالية (HHG) في الماء السائل من المبادئ الأولى.
أخذ العينات الهيكلية: تم أخذ عينات هياكل الماء السائل عند التوازن الحراري (T=300 K) باستخدام نهجين متميزين:
لقطات كلاسيكية: تم إنشاؤها عبر الديناميكا الجزيئية (MD) الكلاسيكية.
لقطات كمومية: تم إنشاؤها عبر (PIMD) باستخدام جهد عصبي لـ Behler-Parrinello مدرب على بيانات revPBE-D3، وتم التحقق من صحته بشكل مستقل باستخدام جهد ذري تعلم آلي (MLIP) مدرب على (RPA). استخدمت محاكاة (PIMD) عدد P=12 من النسخ مع ثرموستات PIGLET غير متوازن لضمان تقارب التوزيع الحراري الكمومي.
الديناميكيات الإلكترونية: لكل لقطة هيكلية مجمدة، تمت محاكاة الاستجابة الإلكترونية لنبضة ليزر قوية (طول موجي 1800 نانومتر، ذروة شدة ∼2×1013 واط سم−2) باستخدام (TDDFT) في الوقت الفعلي (دالة PBE) داخل كود Octopus. تم تثبيت النوى خلال انتشار (TDDFT)، بافتراض الفصل الأديباتي للمقاييس الزمنية حيث تكون الحركة النووية مهملة في مقياس زمن الليزر الذي يبلغ بضع فيمتو ثانية.
حساب الطيف: تم حساب التيار المعتمد على الزمن، وتحويله فوريًا للحصول على السعات الطيفية، ثم تم حساب المتوسط عبر 50 لقطة مستقلة وعدة اتجاهات لليزر لتوليد طيف (HHG) متماثل المناحي (isotropic). كما تم حساب الأطياف المتماسكة (المحددة بالطور) و"غير المتماسكة" (المتوسطة الطور) لعزل تأثيرات التداخل.
النمذجة: لتوضيح الأصول الفيزيائية للميزات الطيفية المرصودة، قام المؤلفون برسم نتائجهم المستمدة من المبادئ الأولى على نموذج "الربط القوي" (Tight-Binding - TB) أحادي البعد ومختزل الأبعاد. مثل هذا النموذج السائل كسلسلة من الذرات مع ضبط المعلمات (المسافات بين الذرات، تكاملات القفز، فجوات النطاق) لتتطابق مع كثافة الحالات الإلكترونية (EDOS) وبيانات الهيكل لكل من الحالات الكلاسيكية والكمومية.
النتائج الرئيسية
الاختلافات الهيكلية: تماشياً مع الأدبيات السابقة، أظهرت اللقطات الكمومية دالة توزيع شعاعي (RDF) داخل جزيئي لـ O-H أكثر اتساعًا بشكل كبير بسبب الحركة النقطية الصفرية. أما بالنسبة للروابط بين الجزيئات، فقد أظهرت (RDF) لـ O-O في الماء الكمومي قمة أولى أكثر حدة عند مسافات أقصر (2.79 أنجستروم مقابل 2.82 أنجستروم) والحد الأدنى الأول أكثر عمقًا، مما يشير إلى تقوية طفيفة للروابط الهيدروجينية على طول إحداثي مشاركة البروتون.
البصمات الطيفية لـ HHG: كانت النتيجة الأكثر إثارة هي التثبيط الكبير في عائد التوافقيات في الحالة الكمومية مقارنة بالحالة الكلاسيكية. وتحديدًا:
ظهر حد أدنى للتداخل (انخفاض في الشدة) حول التوافقي التاسع عشر (بالقرب من قطع الهضبة الأولى) في الطيف الكمومي، وهو ما كان غائبًا أو أضعف بكثير في الطيف الكلاسيكي.
كان الانخفاض في الشدة بعد قطع الهضبة الأولى أكثر وضوحًا في الحالة الكمومية.
كشف تحليل الأطوار المعقدة للسعات الطيفية أن اللقطات الكمومية أظهرت انتشارًا أوسع بكثير في الطور مقارنة باللقطات الكلاسيكية، مما يؤكد أن التثبيط ينشأ من التداخل الهدام بين التكوينات الهيكلية المختلفة.
الاعتماد على الكثافة: أظهرت المحاكاة عند كثافات متفاوتة (900، 996، و1100 كجم/م3) أن تثبيط التوافقي التاسع عشر وموقع القطع حساسان للكثافة. حيث عززت الكثافة الأعلى من التثبيط ونقلت القطع إلى رتب توافقية أقل.
الأصل الفيزيائي عبر نموذج TB: لم يسفر تحليل الارتباط المباشر بين واصفات هيكلية محددة (قمم RDF، فجوات النطاق، تموضع المدارات) وعائد التوافقي التاسع عشر في اللقطات الفردية عن أي ارتباطات ذات دلالة إحصائية. ومع ذلك، نجح نموذج (TB) المبسط في إعادة إنتاج حد التداخل الأدنى. ومن خلال تغيير معلمات النموذج بشكل مستقل، حدد المؤلفون أن التثبيط مدفوع بـ توسيع حالات نطاق التكافؤ والتوصيل (زيادة تكاملات القفز، γv و γc). وهذا التوسيع هو نتيجة مباشرة للتقلبات المعززة في شبكة الروابط الهيدروجينية الناجمة عن التأثيرات النووية الكمومية وزيادة الكثافة.
الأهمية والادعاءات يزعم البحث أن التأثيرات النووية الكمومية لا غنى عنها لالتقاط الديناميكيات الإلكترونية فائقة السرعة في السوائل بدقة، وتحديدًا في سياق (HHG). يوضح العمل أن مطيافية (HHG) للسوائل هي مسبار قابل للتطبيق للتقلبات النووية، حيث تعمل حدود التداخل الدنيا في الطيف كبصمة للانتشار الكمومي للنوى.
يؤكد المؤلفون أن الآلية وراء تثبيط الطيف هي توسيع حالات نطاق التكافؤ الناتج عن تقلبات الروابط الهيدروجينية، وهي آلية من المتوقع أن تعمم على الهياكل الأخرى المرتبطة بالروابط الهيدروجينية. ويؤكدون أن نتائجهم تحل المسألة المتعلقة بما إذا كان بإمكان (HHG) استقصاء الحركة النووية الكمومية في الطور السائل، موضحين أنها تفعل ذلك من خلال تأثيرات التداخل الناشئة عن الهيكل الإلكتروني المتوسط للمجموعة. توفر الدراسة إطارًا حاسوبيًا قويًا (يجمع بين PIMD و TDDFT) ونموذجًا فيزيائيًا مبسطًا لتفسير هذه التأثيرات، مما يطور فهم (HHG) للسوائل بما يتجاوز نظام الاستجابة الخطية. ويشير المؤلفون إلى أنه بينما تفرض محاكاتهم فصلًا أديباتيًا، فإن اعتماد الكثافة والارتباطات الهيكلية المرصودة يؤكدان أن الميزات الطيفية متجذرة في الخصائص الهيكلية للسائل المعدلة بواسطة التأثيرات النووية الكمومية.