ورود به حساب کاربری

نام کاربری *
رمز عبور *
یاداوری

آموزش کد محاسباتی سیستا (مولکول هیدروژن)

امتیاز کاربران

ستاره غیر فعالستاره غیر فعالستاره غیر فعالستاره غیر فعالستاره غیر فعال
 

در این آموزش سعی داریم انرژی کل و طول پیوند مولکول هیدروژن را محاسبه کنیم. برای اینکار فایل فشرده ضمیمه شده را دانلود کنید و آن را در مکان دلخواهتان (در سیستم عامل لینوکس) از حالت فشرده خارج کنید. پس از خارج کردن از حالت فشرده شما سه فایل را میبینید که یکی از آنها فایل h2.fdf است. فایل ورودی برنامه سیستا دارای پسوند fdf است . فایلی که شما دانلود کردید دارای محتوای زیر می باشد:

 

 

که توضیح هر خط آن به این شکل است:

خط اول: اسم سیستم یا ترکیب مورد نظر را اینجا می نویسیم. مهم نیست که چه چیزی باشد فقط چیزی مینویسیم که در آینده که به آن رجوع کردیم بدانیم برای چه سیستمی اجرا گرفتیم.

خط دوم: یک برچسب برای سیستم انتخاب می کنیم

خط سوم: تعداد اتم های ترکیب را اینجا می نویسیم. مثلا برای مثال اتم هیدروژن این تعداد 2 می باشد

خط چهارم: تعداد گونه های اتمی را اینجا می نویسیم. مثلا برای این مثال فقط 1 گونه اتم هیدروژن داریم.

خط ششم الی هشتم: در این بلاک در هر خط یک شاخص، عدد اتمی و برچسب گونه مورد نظر را می نویسم. مثلا در این مثال یک گونه هیدروژن داریم با عدد اتمی 1 که به آن شاخص (ایندکس) 1 نسبت میدهیم.

خط دهم: در این قسمت واحد مختصات اتمی را مینویسم که در اینجا آنگستروم است.

خط یازدهم تا چهاردهم: در این بلاک مختصات اتم های موجود در سیستم را وارد می کنیم. مثلا در این مثال دو اتم هیدروژن داریم که یکی در مختصات 0 0 0 و دیگری در مختصات 0 0 1 (مختصه x برابر 1 است) قرار دارد.

خط شانزدهم: برای توضیح این خط به صفحه 48 فایل user guide کد سیستا مراجعه کنید.

نکته: فایل ورودی برای ترکیب های مختلف با هم تفاوت دارد و برای مثال های بعدی الزاما از همین فایل نمیشه استفاده کرد. به مرور در مثال های بعدی گزینه های دیگه ای که در فایل ورودی میشود به کار برد توضیح داده میشود.

مولکول هیدروژن: اولین اجرا با سیستا

چک کنید همه فایل های مورد نیاز در پوشه مورد نظر وجود داشته باشند.

دقت کنید که علاوه بر فایل ورودی شما به ازای هر گونه اتمی یک فایل شبه پتانسیل نیاز دارید که دارای پسوند .vps یا .psd است.

برای اتم هیدروژن شما می توانید شبه پتانسیل مربوط به تقریب LDA رو از سایت سیستا دانلود کنید.

پس از اطمینان از وجود همه ی فایل های مورد نیاز می توانید با باز کردن ترمینال در مسیری که فایل های ورودی وجود دارند و نوشتن دستور زیر و زدن دکمه اینتر اولین اجرای خود را بگیرید و یک فایل خروجی دریافت کنید:

# siesta < h2.fdf > h2.1.00.out

پس از اجرای این دستور فایل خروجی با نام h2.1.00.out در پوشه جاری ساخته میشه که خروجی برنامه سیستا است.

نکته: در لینوکس علامت > (کوچکتر) به معنای خواندن از ورودی است فایلی که در بعد از آن نوشته می شود به عنوان ورودی برای دستور قبل از این علامت مورد استفاده قرار میگیرد. همچنین علامت < (بزرگتر) به معنای نوشتن در فایل خروجی است. می توان با استفاده از عملگر پایپ (|) همزمان علاوه بر نوشتن در فایل خروجی محتوای آن را بر روی صفحه نمایش دید. برای اینکا به جای دستور فوق باید نوشت:

# siesta < h2.fdf | tee h2.1.00.out

نام فایل خروجی را می توان به صورت دلخواه انتخاب کرد. ولی بهتر طوری انتخاب شود که در آینده قابل فهم باشم. اینجا چون فاصله بین دو اتم هیدروژن یک آنگستروم بود ما اسم فایل خروجی رو به این شکل نوشتیم.

مولکول هیدروژن: نگاهی کلی به فایل های خروجی

حال میخواهیم یک بررسی در فایل های خروجی انجام دهیم ؛

 

مولکول هیدروژن: فاصله ی بین اتمی

هدف: پیدا کردن ساختار تعادلی مولکول

روش اول:

  • محاسبه انرژی کل رای مقادیر مختلف فاصله بین اتمی
  • برازش یک منحنی بر روی نقاط مختلف
  • بدست آوردن کمینه این تابع برازش شده

روش دوم:

  • محاسبه نیرو ها برای فاصله اهی بین اتمی مختلف
  • برازش یک منحنی بر روی نقاط مختلف
  • محاسبه صفر این تابع برازش شده

روش سوم:

  • استفاده از الگوریتمی اتوماتیک برای کمینه کردن انرژی

مولکول هیدروژن: بدست آوردن فاصله ی بین اتمی از روش اول

در فایل ورودی قبلی مقدار فاصله بین اتمی رو از مقدار 0.4 آنگستروم تا 3 آنگستروم با بازه های 0.1 تغییر داده و برای هر کدام یک اجرای جداگانه بگیرید.

change_in_at_dist

بعد از هر بار ذخیره این فایل ورودی باید دستوری مشابه دستور زیر اجرا کنید تا فایل خروجی مربوطه ذخیره شود:

# siesta < h2.fdf > h2.your_interatomic_distance.out

که در اینجا برای هر بار اجرا باید به جای your_interatomic_distance فاصله ای که در فایل ورودی وارد کرده اید بنویسید.

بعد از چندین بار اجرای این دستور شما یک سری فایل خروجی خواهید داشت. حال باید انرژی کل رو از تمام این فایل ها بردارید و در یک فایل ذخیره کنید. برای این کار در مسیری که فایل های خروجی قرار دارند می توانید دستور زیر را وارد کنید:

# grep “Total =” h2.*.out > h2.distance.dat

این دستور این عمل را انجام میدهد : داخل هر فایلی که اول اسم آن h2. و آخر اسمش .out است رفته ، عبارت “Total =” جستجو و پیدا کرده، آن خط را در فایل h2.distance.dat ذخیره می کند.

بعد از این کار فایل h2.distance.dat را باز کنید و فقط دو ستون باقی بگذارید که اولی فاصله بین اتمی و دومی انرژی کل است:

exercise-h2

این اعداد از سیستا نسخه 3 بتا که با کامپایلر g95 کامپایل شده بدست آمد. این اعداد ممکن است برای شما کمی تغییر پیدا کند.

در این مرحله شما باید نمودار برای این اعداد رسم کنید. برای رسم نمودار می توانید از برنامه gnuplot استفاده کنید برای اینکار یک ترمینال در مسیری که فایل مورد نظر شما هست باز کنید و د آن تایپ کنید gnuplot و enter بزنید و بعد دستور زیر را وارد کنید:

> plot "h2.distance.dat" using 1:2 with lines

بعد از اجرای این دستور نمودار به این شکل خواهید داشت:

h2-dist

 

 

-------------------------------------------

منبع : http://www.comphys.ir

شرکت دانش بنیان رایانش سریع هزاره ایرانیان به منظور طراحی، ساخت و گسترش کاربردهای رایانش سریع در علوم مختلف پایه گذاری شده است. این شرکت به دنبال ایجاد بستری با ظرفیت محاسباتی و ذخیره سازی بالا مبتنی بر CPU و GPU می باشد تا بتواند از آن طریق سرویسهای متعددی را به محققین، اساتید، دانشجویان تحصیلات تکمیلی و سایر علاقه مندان این حوزه ارائه نماید.

 

سایت مرتبط : " بنیاد توسعه رایانش سریع و ابری "

 

 

آخرین مقالات

کامپایل و نصب mfix-2016.1

درک عملکرد دستگاه های انرژی، محیط زیست و فرایندها...

ماشین هوشمند بدون کمک انس...

یک ماشین هوشمند با حل معمای مکعب روبیک بدون هیچ آم...

معرفی lammps

اصطلاح شناسی شبیه ساز: فراهم کننده ی بستری برای شب...

آموزش نرم افزار amber (زب...

در اینجا تعدادی از آموزش های آماده شده توسط توسعه...

HPC Moves to the Cloud –...

While some domains that rely on computing systems...

یکسان سازی یا نرمال سازی...

در متن کاوی منظور از پردازش زبان طبیعی (NLP)، قابل...

پردازش متن با Jhazm نسخه...

jHazm نسخه جاوایی هضم از کتابخانه پایتون برای پردا...

Why Demand for In-Memory...

The demand for in-memory computing (IMC) and new...

آموزش عملی گراف کاوی و تح...

networkx یک بسته پایتون برای ایجاد، دستکاری و مطا...

معرفی مختصر مدل اقلیمی CW...

CWRF نسخه توسعه یافته اقلیمی مدل WRF می باشد که د...

بالاخره بعد از سالها یکه...

از سال 2013، ماشین های چینی در رتبه قدرتمندترین اب...

کاربردها و خطرات کلان داد...

داده شامل اطلاعات است؛ اما این همه‌ي ماجرا نیست....

Case Study: Supercomputin...

Hyperion Research has published a new case study...

داده کاوی بزرگ مقرون به ص...

چکیده : کاوش داده ی بزرگ، اغلب نیازمند منابع محاس...

Python: Unlocking the Pow...

The programming language Python is being used to...

Opportunities Abound: HPC...

The is the first entry in a five-part insideHPC s...

آموزش کد محاسباتی سیستا (...

در این آموزش سعی داریم انرژی کل و طول پیوند مولکو...

درمورد مدل‌های آب و هوایی...

دو اصطلاح اقلیم‌شناسی و هوا‌شناسی بیشتر مواقع به‌ج...

الگوریتم تجزیه بندرز خودک...

در نسخه‌ی ۱۲٫۷ نرم‌افزار CPLEX رویه‌ی جدیدی جایگز...

نرم افزار tecplot

معرفی : Tecplot در واقع یک ابزار تجسمی و رسم نمودا...

تجزیه گر یا پارسر متون و...

آپاچی تیکا: آپاچی تیکا یکی از پروژه‌های متن باز A...

OpenFOAM v5 User Guide: I...

This guide accompanies the release of version 5.0...

OpenFOAM User Guide - Ver...

Introduction The OpenFOAM User Guide provides an...

How to Prepare Weather an...

In this video from the NVIDIA GPU Technology Conf...

مثالی از داده کاوی و تحلی...

در این پست یک تحلیل و بصری سازی خوب از گراف ۶۰،۰۰...

نرم افزار سیپلکس

معرفی : Cpelx یک بسته ی نرم افزاری شامل مجموعه ای...

تلاش AMD برای ساخت تراشه‌...

ریسک‌های کمی در دنیای فناوری به بزرگی خرج میلیارده...

معرفی، آموزش و لینک دانلو...

معرفی : ابزاری به منظور آنالیز فنی و طراحی کامپیوت...

Neurala Reduces Training...

Today Neurala announced a breakthrough update to...

با 10 باور نادرست در ارتب...

برای آنکه اطمینان حاصل کنید برنامه‌های سازمان شما...

گالری تصاویر