ويمكن حساب تركيزات المجموعات الجزيئية الضعيفة من الخصائص الكيميائية الحرارية لهياكل الطاقة المنخفضة التي توجد من خلال منهجية أخذ العينات التكوينية المتعددة الخطوات باستخدام خوارزمية جينية وكيمياء الكم شبه التجريبية والكيميائية اللاشرعية.
تتطلب الدراسة الحسابية لتكوين ونمو الهباء الجوي سطح طاقة حر دقيق من جيبس ، والذي يمكن الحصول عليه من الهيكل الإلكتروني لمرحلة الغاز وحسابات التردد الاهتزازي. وهذه الكميات صالحة بالنسبة للمجموعات الجوية التي تتوافق هندستها مع الحد الأدنى على أسطح الطاقة المحتملة لديها. يمكن استخدام الطاقة الحرة جيبس من بنية الطاقة الدنيا للتنبؤ تركيزات الغلاف الجوي للكتلة في ظل مجموعة متنوعة من الظروف مثل درجة الحرارة والضغط. نحن نقدم إجراء غير مكلف حسابيا مبنيا على أخذ عينات التكوينية القائمة على الخوارزميات الجينية تليها سلسلة من حسابات الفحص الدقيقة بشكل متزايد. يبدأ الإجراء بتوليد وتطوير هندسات مجموعة كبيرة من التكوينات باستخدام نماذج شبه تجريبية ثم يصقل الهياكل الفريدة الناتجة في سلسلة من مستويات initio ab عالية المستوى من النظرية. وأخيراً، يتم حساب التصحيحات الحرارية لمجموعة الناتجة من هياكل الطاقة الدنيا وتستخدم لحساب الطاقات الحرة جيبس من التكوين، وثوابت التوازن، وتركيزات الغلاف الجوي. نقدم تطبيق هذا الإجراء على دراسة مجموعات الجليسين الرطبة في ظل الظروف المحيطة.
وأكثر البارامترات غموضا في الدراسات الجوية لتغير المناخ هي المدى الدقيق الذي تعكس به جسيمات السحب الإشعاع الشمسي الوارد. الهباء الجوي، وهي جسيمات معلقة في الغاز، تشكل جزيئات سحابية تسمى نواة التكثيف السحابي (CCN) التي تبعثر الإشعاع الوارد، وبالتالي منع امتصاصه والتدفئة اللاحقة للغلاف الجوي1. ويتطلب الفهم المفصل لهذا التأثير الصافي للتبريد فهم نمو الهباء الجوي إلى شبكات الكربون المتعددة الاستراتيجية، الأمر الذي يتطلب بدوره فهم نمو المجموعات الجزيئية الصغيرة في جزيئات الهباء الجوي. وقد أشارت الأعمال الأخيرة إلى أن تشكيل الهباء الجوي يبدأ بواسطة مجموعات جزيئية قطرها 3 نانومتر أو أقل2؛ ومع ذلك ، فإن نظام الحجم هذا يصعب الوصول إليه باستخدام التقنيات التجريبية3،4. لذلك ، يتم طلب نهج النمذجة الحسابية من أجل التغلب على هذا القيد التجريبي.
باستخدام نهج النمذجة لدينا المذكورة أدناه، يمكننا تحليل نمو أي مجموعة رطبة. لأننا مهتمون بدور الماء في تكوين جزيئات بيولوجية كبيرة من مكونات أصغر في بيئات ما قبل الحيوية ، فإننا نوضح نهجنا مع الجليسين. التحديات التي تواجهها والأدوات اللازمة لمعالجة تلك الأسئلة البحثية هي مشابهة جدا لتلك التي شاركت في دراسة الهباء الجوي ومجموعات prenucleation5،6،7،8,،9،10,،11،12،13،14،15. هنا، نقوم بفحص مجموعات الجليسين الرطبة بدءًا من جزيء جليسين معزول تليه سلسلة من الإضافات التدريجية تصل إلى خمسة جزيئات مائية. الهدف النهائي هو حساب تركيزات التوازن من Gly (H2O)n = 0-5 مجموعات في الغلاف الجوي في درجة حرارة الغرفة في مستوى سطح البحر والرطوبة النسبية (RH) من 100 ٪.
وهناك عدد صغير من هذه المجموعات الجزيئية دون النانومتر تنمو لتصبح كتلة حرجة metastable (1-3 نانومتر في القطر) إما عن طريق إضافة جزيئات بخار أخرى أو تخثر على المجموعات القائمة. هذه المجموعات الحرجة لديها لمحة النمو مواتية مما يؤدي إلى تشكيل أكبر بكثير (تصل إلى 50-100 نانومتر) نواة التكثيف سحابة (CCN)، والتي تؤثر بشكل مباشر على كفاءة هطول الأمطار من الغيوم وكذلك قدرتها على عكس ضوء الحادث. ولذلك، فإن الفهم الجيد للديناميكا الحرارية للمجموعات الجزيئية وتوزيعات توازنها ينبغي أن يؤدي إلى تنبؤات أكثر دقة بأثر الهباء الجوي على المناخ العالمي.
يتطلب النموذج الوصفي لتشكيل الهباء الجوي الديناميكا الحرارية الدقيقة لتشكيل الكتلة الجزيئية. يتطلب حساب الديناميكا الحرارية الدقيقة لتكوين الكتلة الجزيئية تحديد التكوينات الأكثر استقرارًا ، والتي تنطوي على العثور على المينيما العالمية والمحلية على سطح الطاقة المحتمل للكتلة (PES)16. وتسمى هذه العملية أخذ العينات التكوينية ويمكن تحقيقها من خلال مجموعة متنوعة من التقنيات، بما في ذلك تلك القائمة على الديناميات الجزيئية (MD)17،18،,19،20، مونتي كارلو (MC)21،22، والخوارزميات الوراثية (GA)23،24،25.
وقد وضعت بروتوكولات مختلفة على مر السنين للحصول على هيكل والديناميكا الحرارية من هيدرات الغلاف الجوي على مستوى عال من النظرية. واختلفت هذه البروتوكولات في اختيار ‘1’ طريقة أخذ العينات التكوينية، و’2′ طبيعة الأسلوب المنخفض المستوى المستخدم في أخذ العينات التكوينية، و’3′ التسلسل الهرمي للأساليب الرفيعة المستوى المستخدمة لتحسين النتائج في الخطوات اللاحقة.
وشملت أساليب أخذ العينات التكوينية الحدس الكيميائي26، وأخذ العينات العشوائية27،28، الديناميات الجزيئية (MD)29،30، حوض التنقل (BH)31، والخوارزمية الوراثية (GA)24،25،32. وأكثر الطرق المنخفضة المستوى شيوعاً المستخدمة في أساليب أخذ العينات هذه هي حقول القوة أو النماذج شبه التجريبية مثل PM6 و PM7 و SCC-DFTB. وغالبا ما يتبع هذه الحسابات DFT مع مجموعات أساس أكبر على نحو متزايد والوظائف أكثر موثوقية من الدرجات العليا من سلم يعقوب33. في بعض الحالات، ويتبع هذه أساليب الموجي مستوى أعلى مثل MP2، CCSD (T)، وكفاءة التكلفة DLPNO-CCSD (T)34،35.
طوّر كلدغارد وآخرون36 طريقة منهجية يتم فيها إضافة جزيئات الماء عند نقاط على مجالات فيبوناتشي37 حول مجموعات أصغر رطبة أو غير رطبة لتوليد مرشحين لمجموعات أكبر. تتم إزالة المرشحين غير الماديين والزائدين على أساس عتبات الاتصال الوثيق والمسافة الجذرية المتوسطة المربعة بين مختلف المناومين. يتم استخدام التحسينات اللاحقة باستخدام طريقة PM6 شبه التجريبية والتسلسل الهرمي لأساليب DFT والموجات للحصول على مجموعة من أجهزة التوافق في الطاقة المنخفضة على مستوى عال ٍ من النظرية.
مستعمرة النحل الاصطناعية (ABC) خوارزمية38 هو نهج جديد أخذ العينات التكوينية التي تم تنفيذها مؤخرا من قبل تشانغ وآخرون لدراسة المجموعات الجزيئية في برنامج يسمى ABCluster39. Kubecka وآخرون40 تستخدم ABCluster لأخذ العينات التكوينية تليها عمليات إعادة تحسين منخفضة المستوى باستخدام ضيق ملزمة GFN-xTB شبه التجريبية الطريقة41. كما قاموا بصقل الهياكل والطاقات باستخدام أساليب DFT تليها الطاقات النهائية باستخدام DLPNO-CCSD (T).
بغض النظر عن الأسلوب، يبدأ أخذ العينات التكوينية بتوزيع عشوائي أو غير عشوائي للنقاط على PES. كل نقطة تتوافق مع هندسة محددة من الكتلة الجزيئية المعنية ويتم إنشاؤها بواسطة طريقة أخذ العينات. ثم تم العثور على أقرب الحد الأدنى المحلي لكل نقطة باتباع اتجاه “انحدار” على PES. مجموعة من minima وجدت بالتالي تتوافق مع تلك الهندسة من الكتلة الجزيئية التي هي مستقرة، على الأقل لبعض الوقت. هنا ، فإن شكل PES وتقييم الطاقة في كل نقطة على السطح تكون حساسة للوصف المادي للنظام حيث وصفا ماديا أكثر دقة يؤدي إلى حساب الطاقة أكثر تكلفة حسابيا. وسوف نستخدم على وجه التحديد طريقة الجمعية العامة التي تم تنفيذها في برنامج OGOLEM25 ، والتي تم تطبيقها بنجاح على مجموعة متنوعة من مشاكل التحسين والعينات التكوينية العالمية42،43،44،4544،لتوليد المجموعة الأولية من نقاط أخذ العينات. سيتم وصف PES من قبل نموذج PM746 الذي تم تنفيذه في برنامج MOPAC201647. يتم استخدام هذا المزيج لأنه يولد مجموعة أكبر من النقاط مقارنة بأساليب MD و MC ويجد minima المحلية أسرع من الأوصاف الأكثر تفصيلاً للPES.
يتم أخذ مجموعة من minima المحلية المحسنة من قبل الجمعية العامة كهندسة البداية لسلسلة من خطوات الفحص ، والتي تؤدي إلى مجموعة من الحد الأدنى من الطاقة المنخفضة. يبدأ هذا الجزء من البروتوكول بتحسين مجموعة من الهياكل المثلى لـ GA الفريدة باستخدام نظرية الكثافة الوظيفية (DFT) مع مجموعة أساس صغيرة. هذه المجموعة من التحسينات سوف تعطي عموما مجموعة أصغر من الهياكل الدنيا المحلية الفريدة التي هي على غرار بمزيد من التفصيل بالمقارنة مع الهياكل شبه التجريبية الأمثل من قبل الجمعية العامة. ثم يتم تنفيذ جولة أخرى من تحسينات DFT على هذه المجموعة الأصغر من الهياكل باستخدام مجموعة أساس أكبر. مرة أخرى، هذه الخطوة عموما ً ستعطي مجموعة أصغر من الهياكل الفريدة التي تم تصميمها بمزيد من التفصيل مقارنة بخطوة DFT الصغيرة. ثم يتم تحسين المجموعة النهائية من الهياكل الفريدة إلى تقارب أكثر إحكاماويتم حساب ترددات الاهتزازات التوافقية. بعد هذه الخطوة لدينا كل ما نحتاجه لحساب تركيزات التوازن للمجموعات في الغلاف الجوي. يتم تلخيص النهج العام بشكل شكلي في الشكل 1. سنستخدم PW9148 تعميم التدرج التقريب (GGA) تبادل الارتباط الوظيفي في Gaussian0949 تنفيذ DFT جنبا إلى جنب مع اثنين من الاختلافات من مجموعة أساس Pople50 (6-31 + G * لخطوة أساس صغير و 6-311 + + G ** لخطوة أساس كبير). تم اختيار هذا المزيج الخاص من تبادل الارتباط الوظيفي ومجموعة الأساس نظرا لنجاحها السابق في الحوسبة الدقيقة جيبس الطاقات الحرة لتشكيل لمجموعات الغلاف الجوي51,52.
يفترض هذا البروتوكول أن المستخدم لديه إمكانية الوصول إلى مجموعة حوسبة عالية الأداء (HPC) مع نظام الدفعي المحمول53 (PBS) و MOPAC2016 (http://openmopac.net/MOPAC2016.html)47و OGOLEM (https://www.ogolem.org)25و Gaussian 09 (https://gaussian.com)49وبرنامج OpenBabel54 (http://openbabel.org/wiki/Main_Page) المثبت باتباع إرشادات التثبيت الخاصة بهم. تستخدم كل خطوة في هذا البروتوكول أيضاً مجموعة من البرامج النصية shell و Python 2.7 الداخلية التي يجب حفظها إلى دليل مضمن في المتغير البيئي $PATH للمستخدم. يجب أيضاً تحميل كافة الوحدات النمطية البيئية الضرورية وأذونات التنفيذ لتشغيل كافة البرامج المذكورة أعلاه في جلسة عمل المستخدم. القرص واستخدام الذاكرة من قبل رمز الجمعية العامة (OGOLEM) والرموز شبه التجريبية (MOPAC) صغيرة جدا وفقا لمعايير موارد الكمبيوتر الحديثة. يعتمد الاستخدام العام للذاكرة والقرص لـ OGOLEM/MOPAC على عدد مؤشرات الترابط التي يريد المرء استخدامها، وحتى ذلك الحين، سيكون استخدام الموارد صغيرًا مقارنة بقدرات معظم أنظمة HPC. تعتمد احتياجات الموارد لأساليب QM على حجم المجموعات ومستوى النظرية المستخدمة. ميزة استخدام هذا البروتوكول هو أن المرء يمكن أن تختلف مستوى النظرية لتكون قادرة على حساب المجموعة النهائية من هياكل الطاقة المنخفضة، مع الأخذ في الاعتبار أن الحسابات الأسرع عادة ما تؤدي إلى مزيد من عدم اليقين في دقة النتائج.
من أجل الوضوح، سيتم الإشارة إلى الكمبيوتر المحلي للمستخدم باسم“الكمبيوتر المحلي”في حين سيتم الإشارة إلى الكتلة HPC لديهم حق الوصول إلى“كتلة البعيد”.
وتتوقف دقة البيانات التي يولدها هذا البروتوكول أساساً على ثلاثة أمور: ‘1’ تنوع التكوينات التي أخذت عينات منها الخطوة 2، ‘2’ دقة الهيكل الإلكتروني للنظام، ‘3’ ودقة التصحيحات الحرارية. يمكن معالجة كل من هذه العوامل عن طريق تعديل الأسلوب عن طريق تحرير البرامج النصية المضمنة. يتم التغلب بسهولة على العامل الأول باستخدام مجموعة أولية أكبر من الهياكل التي تم إنشاؤها عشوائيًا ، وتكرارات أكثر عددًا من الجمعية العامة ، وتعريف أكثر مرونة للمعايير المشاركة في الجمعية العامة. بالإضافة إلى ذلك ، يمكن للمرء استخدام طريقة مختلفة شبه تجريبية مثل كثافة الشحنة المتسقة ذاتيًا – المحدودة -الربط المحكم (SCC-DFTB)62 ونموذج السعة الفعالة للتجزئة (EFP)63 من أجل استكشاف آثار الأوصاف الفيزيائية المختلفة. القيد الرئيسي هنا هو عدم قدرة الطريقة على تشكيل أو كسر السندات الكوفانية ، مما يعني أن يتم تجميد مونومرات. يجد الإجراء GA فقط المواقف النسبية الأكثر استقرارا من هذه مونومرات المجمدة وفقا للوصف شبه التجريبية.
ويمكن تحسين دقة الهيكل الإلكتروني للنظام بطرق متنوعة، لكل منها تكلفته ابتيالية. يمكن للمرء أن يختار أفضل كثافة وظيفية، مثل M06-2X64 وwB97X-V65،أو الكم الميكانيكية (QM) طريقة مثل Møller-Plesset66،,67،,68 (MPn) نظريات الاضطرابات والمقترنة الكتلة69 (CC) أساليب من أجل تحسين الوصف المادي للنظام. في التسلسل الهرمي للوظائف ، يتحسن الأداء بشكل عام عند الانتقال من وظائف تقريب التدرج المعمم (GGA) مثل PW91 إلى الوظائف الهجينة المنفصلة عن النطاق مثل wB97X-D والوظائف الهجينة الفوقية GGA مثل M06-2X.
وعيب أساليب DFT هو أن التقارب المنهجي نحو قيمة دقيقة غير ممكن؛ ومع ذلك ، أساليب DFT غير مكلفة حسابيا وهناك مجموعة واسعة من الوظائف لمجموعة واسعة من التطبيقات.
الطاقات المحسوبة باستخدام أساليب الموجات مثل MP2 و CCSD (T) بالاقتران مع مجموعات أساس متسقة للارتباط من العدد الكاردينالي المتزايد ([aug-] cc-pV [D, T, Q,…] Z) تتلاقى نحو حد مجموعة الأساس الكامل بشكل منهجي ، ولكن التكلفة الحسابية لكل حساب تصبح باهظة مع نمو حجم النظام. ويمكن إنجاز المزيد من التنقيح للهيكل الإلكتروني باستخدام مجموعات أساس مترابطة صراحة70 وبالاستقراء على حد الأساس الكامل (CBS)71. عملنا الأخير يشير إلى أن كثافة المجهزة بشكل صريح مترابطة من الدرجة الثانية مولر-Plesset (DF-MP2-F12) نهج مضطرب ة تسفر عن طاقات تقترب من أن من حسابات MP2 / CBS32. تعديل البروتوكول الحالي لاستخدام أساليب بنية إلكترونية مختلفة ينطوي على خطوتين: (1) إعداد ملف إدخال قالب بعد بناء الجملة التي قدمها البرنامج، (2) وتحرير تشغيل pw91-sb.csh، تشغيل pw91-lb.csh،والبرامج النصية تشغيل pw91-lb-ultrafine.csh لتوليد بناء جملة ملف الإدخال الصحيح وكذلك النصي تقديم الصحيح للبرنامج.
وأخيراً، تعتمد دقة التصحيحات الحرارية على طريقة الهيكل الإلكتروني وكذلك على وصف نظام المعلومات الجغرافية حول الحد الأدنى العالمي. يتطلب الوصف الدقيق للPES حساب المشتقات الثالثة والعليا من PES فيما يتعلق بالتشريد في الدرجات النووية للحرية ، مثل حقل القوة الكوار72،73 (QFF) ، وهي مهمة مكلفة بشكل استثنائي. يستخدم البروتوكول الحالي تقريب المذبذب التوافقي إلى ترددات الاهتزازية ، مما يؤدي إلى الحاجة إلى حساب ما يصل فقط إلى المشتقات الثانية من PES. هذا النهج يصبح إشكالية في النظم مع الانسجام عالية، مثل الجزيئات المرنة جدا ومتماثل ة مزدوجة البئر المحتملة بسبب الفرق الكبير في PES الحقيقي وPES التوافقي. وعلاوة على ذلك، فإن تكلفة وجود نظام إدارة المؤشرات الإلكترونية عالي الجودة من طريقة بنية إلكترونية تتطلب استخداماً حسابياً لا يؤدي إلا إلى تفاقم مشكلة تكلفة حسابات التردد الاهتزازي. ويتمثل أحد النهوج المتبعة في التغلب على ذلك في استخدام الطاقات الإلكترونية المستمدة من حساب هيكل إلكتروني عالي الجودة إلى جانب ترددات اهتزازية محسوبة على نظام PES أقل جودة، مما يؤدي إلى تحقيق توازن بين التكلفة والدقة. ويمكن تعديل البروتوكول الحالي لاستخدام أوصاف PES مختلفة كما هو موضح في الفقرة السابقة؛ ومع ذلك، يمكن للمرء أيضاتحرير الكلمات الرئيسية تردد اهتزازي في البرامج النصية والقوالب لحساب ترددات اهتزازية أنوينية.
وثمة مسألتان حاسمتان لأي بروتوكول لأخذ العينات التكوينية هما الطريقة الأولية لأخذ عينات من سطح الطاقة المحتمل والمعايير المستخدمة لتحديد كل مجموعة. لقد استخدمنا على نطاق واسع مجموعة متنوعة من الأساليب في عملنا السابق. بالنسبة للمسألة الأولى ، وهي الطريقة الأولية لأخذ عينات من سطح الطاقة المحتمل ، قمنا باختيار استخدام GA مع أساليب شبه تجريبية استنادًا إلى هذه العوامل. أخذ العينات التكوينية باستخدام الحدس الكيميائي26، وأخذ العينات العشوائية ، والديناميات الجزيئية (MD)29،30، تفشل في العثور على minima العالمية المفترضة بانتظام لمجموعات أكبر من 10 مونومرات ، كما لاحظنا في دراساتنا لمجموعات المياه18. لقد استخدمنا بنجاح حوض التنقل (BH) لدراسة PES معقدة من (H2O)1174، لكنه يتطلب إدراج دليل بعض ايزومرات الطاقة المنخفضة المحتملة خوارزمية BH لم تجد. وأثبتت مقارنة أداء البوسنة والهرسك والجمعية العامة في العثور على الحد الأدنى العالمي لمجموعات المياه(H2O)n=10-20 أن الجمعية العامة وجدت باستمرار الحد الأدنى العالمي أسرع من75BH . الجمعية العامة كما نفذت في OGOLEM وCLUSTER هو تنوعا جدا لأنه يمكن تطبيقها على أي كتلة الجزيئية ويمكن أن تتفاعل مع عدد كبير من الحزم مع حقل القوة الكلاسيكية، شبه التجريبية، كثافة وظيفية، وقدرات ab initio. يتم اختيار PM7 من خلال سرعته اما الدقة المعقولة. أي طريقة شبه تجريبية أخرى تقريبا سيكون لها تكلفة حسابية أعلى بكثير.
أما بالنسبة للمسألة الثانية ، فقد استكشفنا استخدام معايير مختلفة لتحديد هياكل فريدة تتراوح بين الطاقات الإلكترونية ، ولحظات ثنائي القطب ، وتداخل RMSDs والثوابت التناوبية. استخدام لحظات ثنائي القطب ثبت من الصعب لأن كل من مكونات لحظة ثنائي القطب كانت تعتمد على اتجاه الجزيء واللحظة الديبول الكلي كان حساسا جدا لاختلافات الهندسة في مثل هذه الطريقة التي كان من الصعب تحديد عتبات تحديد هو الهياكل هي نفسها أو فريدة من نوعها. وثبت أن الجمع بين الطاقات الإلكترونية والثوابت التناوبية هو الأكثر فائدة.
وتستند المعايير الحالية لاعتبار هيكلين فريدين على عتبة فرق الطاقة 0.10 سعرحراري-1 وفارق ثابت دوراني قدره 1 في المائة. لذلك ، يعتبر هيكلان مختلفان إذا كانت طاقاتهما تختلف بأكثر من 0.10 سعرحراريمول -1 (~ 0.00015 وحدة حرارية) وأي من ثوابتها الدورانية الثلاثة (A، B، C) تختلف بأكثر من 1٪. ووجدت المعايير الداخلية الهامة على مر السنين أن هذه العتبات هي خيارات معقولة. وقد تم تطبيق نهجنا أخذ العينات التكوينية ومنهجية الفحص إلى مجموعات مقيدة ضعيفة للغاية مثل الهيدروكربونات متعددة العطريات المعقدة بالماء76،77 وكذلك هيدرات كبريتات ternary ملزمة بقوة تحتوي على الأمونيا والأمينات32. بالنسبة للمجموعات التي توجد فيها حالات بروتون مختلفة يجب النظر فيها ، فإن أفضل نهج هو إجراء حسابات الجمعية العامة المختلفة ، كل منها يبدأ بمونومرات في حالات البروتون المختلفة. وهذا يضمن أن الهياكل ذات حالات البروتونات المختلفة يتم النظر فيها بعناية. ومع ذلك ، فإن حسابات DFT منخفضة المستوى غالبا ً ما تسمح لدول البروتونات بالتغيير أثناء تحسين الهندسة ، وبالتالي تسفر عن حالة البروتونالأكثر استقرارًا بغض النظر عن هندسة البداية.
يجب أن تعمل طرق أخذ العينات التكوينية الخاصة بنا بشكل جيد حتى بالنسبة للجزيئات المرنة طالما أن رموز GA مُربطة بأساليب عامة غير معلمة تسمح للمونومرات باعتماد تكوينات مختلفة أثناء تشغيل GA. على سبيل المثال، من شأن ربط الجمعية العامة مع PM7 أن يسمح بتغيّر هياكل مونومرات، ولكن إذا انكسرت روابطها كما سيحدث عندما تتغير دول البروتونات، فقد يتم التخلص من الهياكل كمرشحين غير مقبولين.
لقد نظرنا في طرق مختلفة لتصحيح أوجه القصور في التقريب التوافقي ، وخاصة تلك الناشئة عن الترددات الذبذبية المنخفضة. وليس من الصعب إدراج التقريب شبه التوافقي في المنهجية الحالية. ومع ذلك ، لا تزال هناك أسئلة حول الطريقة شبه التوافقية ، خاصة عندما يتعلق الأمر بتردد القطع الذي سيتم تطبيقه أدناه. أيضا، لا توجد أعمال قياس صارمة دراسة موثوقية تقريب شبه RRHO على الرغم من أن الحكمة التقليدية تشير إلى أنه ينبغي أن يكون تحسنا على تقريب RRHO.
ويمكن تعميم البروتوكول المعروض على هذا النحو على أي نظام من الكتل الجزيئية لمرحلة الغاز غير المنضمة. كما يمكن تعميم استخدام أي طريقة شبه تجريبية، وطريقة الهيكل الإلكتروني ة والبرمجيات، وطريقة التحليل الاهتزازي والبرمجيات عن طريق تحرير البرامج النصية والقوالب. هذا يفترض أن المستخدم مرتاح مع واجهة سطر الأوامر لينكس ، الكتابة النصية بيثون ، والحوسبة عالية الأداء. بناء الجملة غير مألوفة ونظرة من نظام التشغيل لينكس وعدم وجود خبرة البرمجة هو أكبر عثرة في هذا البروتوكول، وحيث الطلاب الجدد النضال أكثر. وقد استخدم هذا البروتوكول بنجاح في مجموعة متنوعة من التطبيقات لسنوات في مجموعتنا، مع التركيز في الغالب على آثار حمض الكبريتيك والأمونيا على تكوين الهباء الجوي. وستشمل التحسينات الإضافية التي أدخلت على هذا البروتوكول وجود واجهة أقوى لبرامج الهيكل الإلكترونية، والتطبيقات البديلة للخوارزمية الجينية، وربما استخدام أساليب أحدث لإجراء حسابات أسرع للطاقات الإلكترونية والاهتزازية. تطبيقاتنا الحالية لهذا البروتوكول هي استكشاف أهمية الأحماض الأمينية في المراحل المبكرة من تكوين الهباء الجوي في الغلاف الجوي الحالي وفي تشكيل جزيئات بيولوجية أكبر في بيئات البريبايوتيك.
The authors have nothing to disclose.
تم دعم هذا المشروع من خلال المنح CHE-1229354، CHE-1662030، CHE-1721511، وCHE-1903871 من المؤسسة الوطنية للعلوم (GCS)، وجائزة أرنولد ومابيل بيكمان للباحثين (AGG)، ومنحة باري م. جولدووتر (AGG). واستُخدمت موارد الحوسبة العالية الأداء التابعة لاتحاد الزئبق (http://www.mercuryconsortium.org).
Avogadro | https://avogadro.cc | Open-source molecular visualization program | |
Gaussian [09/16] Software | http://www.gaussian.com/ | Commercial ab initio electronic structure program | |
MOPAC 2016 | http://openmopac.net/MOPAC2016.html | Open-source semi-empirical program | |
OGOLEM Software | https://www.ogolem.org | Genetic algorithm-based global optimization program | |
OpenBabel | http://openbabel.org/wiki/Main_Page | Open-source cheminformatics library | |
calcRotConsts.py | Shields Group, Department of Chemistry, Furman University | Python script to compute rotational constants | |
calcSymmetry.csh | Shields Group, Department of Chemistry, Furman University | Shell script to calculate symmetry number of a molecule given Cartesian coordinates | |
combine-GA.csh | Shields Group, Department of Chemistry, Furman University | Shell script to combine energy and rotational constants from different GA directories | |
combine-QM.csh | Shields Group, Department of Chemistry, Furman University | Shell script to combine energy and rotational constants from different QM directories | |
gaussianE.csh | Shields Group, Department of Chemistry, Furman University | Shell script to extract Gaussian 09 energies | |
gaussianFreqs.csh | Shields Group, Department of Chemistry, Furman University | Shell script to extract Gaussian 09 vibrational frequencies | |
getrotconsts | Shields Group, Department of Chemistry, Furman University | Executable to calculate rotational constants given a molecule's Cartesian coordinates | |
getRotConsts-dft-lb.csh | Shields Group, Department of Chemistry, Furman University | Shell script to compute rotational constants for a batch of large basis DFT optimized structures | |
getRotConsts-dft-lb-ultrafine.csh | Shields Group, Department of Chemistry, Furman University | Shell script to compute rotational constants for a batch of ultrafine DFT optimized structures | |
getRotConsts-dft-sb.csh | Shields Group, Department of Chemistry, Furman University | Shell script to compute rotational constants for a batch of small basis DFT optimized structures | |
getRotConsts-GA.csh | Shields Group, Department of Chemistry, Furman University | Shell script to compute rotational constants for a batch of genetic algorithm optimized structures | |
global-minimum-coords.xyz | Shields Group, Department of Chemistry, Furman University | Cartesian coordinates of global minimum structures of gly-(h2o)n, where n=0-5 | |
make-thermo-gaussian.csh | Shields Group, Department of Chemistry, Furman University | Shell script to extract data from Gaussian output files and make input files for the thermo.pl script | |
ogolem-input-file.ogo | Shields Group, Department of Chemistry, Furman University | Ogolem sample input file | |
ogolem-submit-script.pbs | Shields Group, Department of Chemistry, Furman University | PBS batch submission file for Ogolem calculations | |
README.docx | Shields Group, Department of Chemistry, Furman University | Clarifications to help readers use the scripts effectively | |
runogolem.csh | Shields Group, Department of Chemistry, Furman University | Shell script to run OGOLEM | |
run-pw91-lb.csh | Shields Group, Department of Chemistry, Furman University | Shell script to run a batch of large basis DFT optimization calculations | |
run-pw91-lb-ultrafine.csh | Shields Group, Department of Chemistry, Furman University | Shell script to run a batch of ultrafine DFT optimization calculations | |
run-pw91-sb.csh | Shields Group, Department of Chemistry, Furman University | Shell script to run a batch of small basis DFT optimization calculations | |
run-thermo-pw91.csh | Shields Group, Department of Chemistry, Furman University | Shell script to compute the thermodynamic corrections for a batch of DFT optimized structures | |
similarityAnalysis.py | Shields Group, Department of Chemistry, Furman University | Python script to determine unique structures based on rotational constants and energies | |
symmetry | Shields Group, Department of Chemistry, Furman University | Executable to calculate molecular symmetry given Cartesian coordinates | |
symmetry.c | (C) 1996, 2003 S. Patchkovskii, Serguei.Patchkovskii@sympatico.ca | C code to determine the molecular symmstry of a molecule given Cartesian coordinates | |
template-marcy.pbs | Shields Group, Department of Chemistry, Furman University | Template for a PBS submit script which uses OGOLEM | |
template-pw91.com | Shields Group, Department of Chemistry, Furman University | Template Gaussian 09 input | |
template-pw91-HL.com | Shields Group, Department of Chemistry, Furman University | Template Gaussian 09 input for ultrafine DFT optimization | |
thermo.pl | https://www.nist.gov/mml/csd/chemical-informatics-research-group/products-and-services/program-computing-ideal-gas | Perl open-source script to compute ideal gas thermodynamic corrections | |
gly-h2o-n.xlsx | Shields Group, Department of Chemistry, Furman University | Excel spreadsheet for the complete protocol | |
table-1.xlsx | Shields Group, Department of Chemistry, Furman University | Excel spreadsheet | |
table-2.xlsx | Shields Group, Department of Chemistry, Furman University | Excel spreadsheet | |
table-3.xlsx | Shields Group, Department of Chemistry, Furman University | Excel spreadsheet | |
table-4.xlsx | Shields Group, Department of Chemistry, Furman University | Excel spreadsheet | |
water.xyz | Shields Group, Department of Chemistry, Furman University | Cartesian coordinates of water | |
glycine.xyz | Shields Group, Department of Chemistry, Furman University | Cartesian coordinates of glycine |