راهنمای کامل شبیه‌سازی کوانتومی با ChemDoodle 3D، ORCA و Gaussian

از رسم ساختار تا تفسیر نتایج — آموزش گام‌به‌گام برای دانشجویان و پژوهشگران ۱۴۰۵

اگر در حوزه شیمی محاسباتی، شیمی آلی فیزیکی، یا مواد پیشرفته فعالیت می‌کنید، به احتمال زیاد با این چالش آشنا هستید: رسم ساختار سه‌بعدی در یک نرم‌افزار، خروجی گرفتن برای یک پکیج محاسباتی، اجرای محاسبه، و سپس برگرداندن نتایج برای تجزیه‌وتحلیل — همه در یک زنجیره یکپارچه. این راهنما دقیقاً همین مسیر را پوشش می‌دهد: از رسم مولکول در ChemDoodle 3D تا ارسال فایل ورودی به ORCA یا Gaussian، اجرای محاسبه روی سرور یا کلاستر دانشگاه، و بازگشت نتایج به ChemDoodle برای تجسم.

این مقاله هم برای دانشجویان کارشناسی ارشد و دکتری که تازه وارد دنیای شیمی کوانتومی می‌شوند مناسب است، هم برای اعضای هیئت علمی و پژوهشگران ارشد که به دنبال بهینه‌سازی گردش کار محاسباتی خود هستند.

فهرست مطالب

1. معرفی سه نرم‌افزار و نقش هر کدام در زنجیره محاسباتی
2. پیش‌نیازها و نصب
3. گام ۱ — رسم و بهینه‌سازی اولیه ساختار در ChemDoodle 3D
4. گام ۲ — خروجی فایل ورودی ORCA
5. گام ۳ — خروجی فایل ورودی Gaussian
6. گام ۴ — اجرای محاسبه در ORCA
7. گام ۵ — اجرای محاسبه در Gaussian
8. گام ۶ — بازگشت نتایج به ChemDoodle 3D برای تجسم
9. مثال کاربردی: بهینه‌سازی هندسه
10. مثال کاربردی: پیش‌بینی طیف NMR
11. دانلود و لایسنس
12. سوالات متداول

۱. معرفی سه نرم‌افزار و نقش هر کدام

پیش از شروع آموزش، مهم است بدانید هر نرم‌افزار در این زنجیره چه نقشی دارد. این سه ابزار مکمل یکدیگرند، نه جایگزین.

ChemDoodle 3D — محیط رسم و تجسم

ChemDoodle 3D ابزار رسم، ویرایش، و تجسم سه‌بعدی مولکول‌هاست. در نسخه ۷.۷ (دسامبر ۲۰۲۵)، پشتیبانی از خواندن و نوشتن فایل‌های ورودی ORCA به طور رسمی اضافه شد — این یعنی می‌توانید مستقیم از ChemDoodle 3D یک فایل .inp برای ORCA تولید کنید. همچنین خروجی Gaussian (فرمت .gjf / .com) از قبل پشتیبانی می‌شد. در کنار اینها، ChemDoodle 3D می‌تواند فایل‌های خروجی (.log, .out) هر دو نرم‌افزار را بخواند و هندسه بهینه‌شده، مد‌های ارتعاشی، و سطح مولکولی را نمایش دهد.

ORCA — پکیج شیمی کوانتومی رایگان و قدرتمند

ORCA یک پکیج شیمی کوانتومی رایگان برای استفاده آکادمیک است که توسط گروه Frank Neese توسعه یافته. ORCA طیف گسترده‌ای از روش‌ها را پوشش می‌دهد: DFT، HF، MP2، CCSD(T)، محاسبات حالت برانگیخته (TD-DFT, EOM-CCSD)، پیش‌بینی NMR و EPR، و محاسبات جفت‌شده نور-ماده. برای دانشگاه‌های ایران که به دنبال جایگزین رایگان برای Gaussian هستند، ORCA گزینه اول است.

Gaussian — استاندارد صنعتی شیمی محاسباتی

Gaussian قدیمی‌ترین و پراستنادترین پکیج شیمی کوانتومی جهان است. در بسیاری از مجلات معتبر، نتایج Gaussian به عنوان مرجع پذیرفته می‌شوند. Gaussian نرم‌افزار تجاری با لایسنس سالانه یا دائم است و برای دانشگاه‌ها و مراکز تحقیقاتی لایسنس سایتی ارائه می‌دهد. اگر مقاله‌ای با داده‌های Gaussian دارید یا استاد راهنمایتان از Gaussian استفاده می‌کند، این راهنما نحوه اتصال آن به ChemDoodle 3D را نشان می‌دهد.

۲. پیش‌نیازها و نصب

ChemDoodle 3D

ChemDoodle 3D بخشی از بسته ChemDoodle است و روی ویندوز، macOS، و لینوکس اجرا می‌شود. نسخه ۷.۷ آخرین نسخه پایدار در سال ۱۴۰۵ است. نصب آن نیاز به هیچ runtime یا dependency خارجی ندارد. فایل نصاب را از سایت رسمی iChemLabs یا از طریق کانال تلگرام ما دریافت کنید.

ORCA

ORCA از طریق سایت رسمی orcaforum.kofo.mpg.de برای استفاده آکادمیک رایگان قابل دانلود است. برای دانلود نیاز به ثبت‌نام دارید. نسخه‌های ویندوز، لینوکس (x86_64)، و macOS موجود هستند. برای اجرای موازی، نسخه OpenMPI را انتخاب کنید. روی کلاسترهای دانشگاهی معمولاً ORCA از پیش نصب است — با مدیر سیستم بررسی کنید.

Gaussian

Gaussian نرم‌افزار تجاری است و نصب آن معمولاً توسط IT دانشگاه یا مرکز تحقیقاتی انجام می‌شود. اگر به لایسنس Gaussian نیاز دارید، از طریق تلگرام با ما تماس بگیرید.

💬 برای دریافت قیمت و مشاوره رایگان → همین الان در تلگرام پیام دهید — پاسخ در کمتر از چند ساعت

گام ۱ — رسم و بهینه‌سازی اولیه ساختار در ChemDoodle 3D

اولین قدم در هر محاسبه شیمی کوانتومی، داشتن یک هندسه اولیه معقول است. ChemDoodle 3D چند روش برای این کار ارائه می‌دهد:

روش اول: رسم مستقیم در ChemDoodle 3D

1. ChemDoodle 3D را باز کنید.
2. از نوار ابزار، ابزار Add Atom را انتخاب کنید.
3. اتم‌ها را روی بوم کلیک کرده و با ابزار Add Bond آنها را به هم وصل کنید.
4. نوع اتم را از جدول تناوبی پانل کناری انتخاب کنید.
5. پس از رسم، از منوی Structure > Clean 3D Structure استفاده کنید تا هندسه اولیه با force field بهینه شود.

روش دوم: وارد کردن از ChemDoodle 2D

1. ساختار دو‌بعدی را در ChemDoodle 2D رسم کنید.
2. از منوی File > Send to ChemDoodle 3D استفاده کنید.
3. ChemDoodle 3D به صورت خودکار مختصات سه‌بعدی تولید می‌کند.
4. دوباره Clean 3D Structure را اجرا کنید.

روش سوم: وارد کردن از PubChem یا فایل موجود

1. از منوی File > Search PubChem نام یا CAS number مولکول را جستجو کنید.
2. ساختار سه‌بعدی مستقیماً بارگذاری می‌شود.
3. فرمت‌های .xyz، .pdb، .mol، و .sdf هم مستقیم قابل باز شدن هستند.

بهینه‌سازی force field پیش از محاسبه کوانتومی

قبل از ارسال ساختار به ORCA یا Gaussian، توصیه می‌شود یک بهینه‌سازی سریع با force field انجام دهید:

1. از منوی Structure > Optimize Geometry استفاده کنید.
2. ChemDoodle 3D از MMFF94 یا UFF (بسته به مولکول) استفاده می‌کند.
3. این مرحله هندسه را به یک minimum محلی نزدیک می‌کند و از diverge شدن محاسبه کوانتومی جلوگیری می‌کند.

گام ۲ — خروجی فایل ورودی ORCA

از نسخه ۷.۷، ChemDoodle 3D می‌تواند مستقیماً فایل .inp برای ORCA تولید کند:

1. ساختار بهینه‌شده را در ChemDoodle 3D باز کنید.
2. از منوی File > Save As استفاده کنید.
3. فرمت ORCA Input File (*.inp) را انتخاب کنید.
4. فایل را ذخیره کنید.

فایل خروجی ساختار زیر را دارد:

# ORCA Input generated by ChemDoodle 3D 7.7
! B3LYP def2-SVP Opt Freq

* xyz 0 1
C   0.000000   0.000000   0.000000
H   0.629118   0.629118   0.629118
...
*

ویرایش فایل ورودی ORCA

فایل تولید‌شده را با یک ویرایشگر متنی (Notepad++، VS Code، یا vim) باز کنید و موارد زیر را تنظیم کنید:

• روش محاسباتی: به جای B3LYP می‌توانید PBE0، wB97X-D3، یا DLPNO-CCSD(T) قرار دهید.
• basis set: def2-SVP برای بهینه‌سازی سریع، def2-TZVP برای انرژی دقیق‌تر.
• بار و تعداد الکترون‌های جفت‌نشده: عدد اول در * xyz 0 1 بار مولکول و عدد دوم multiplicity است.
• حافظه و پردازنده: اضافه کنید: %maxcore 4000 و %pal nprocs 8 end

# مثال کامل برای بهینه‌سازی هندسه و محاسبه فرکانس
! B3LYP def2-TZVP Opt Freq RIJCOSX def2/J

%maxcore 4000
%pal nprocs 8 end

* xyz 0 1
C   0.000000   0.000000   0.000000
...
*

گام ۳ — خروجی فایل ورودی Gaussian

برای Gaussian، فرآیند مشابه است:

1. از منوی File > Save As استفاده کنید.
2. فرمت Gaussian Input File (*.gjf) یا (*.com) را انتخاب کنید.
3. فایل را ذخیره کنید.

ساختار فایل خروجی:

%chk=molecule.chk
%mem=8GB
%nprocshared=8
#p B3LYP/6-31G(d) Opt Freq

ChemDoodle 3D Generated Input

0 1
C   0.000000   0.000000   0.000000
H   0.629118   0.629118   0.629118
...

تنظیمات مهم فایل Gaussian

• route section: خط #p B3LYP/6-31G(d) Opt Freq را بر اساس نیاز تغییر دهید. برای محاسبه NMR اضافه کنید: NMR=GIAO
• حافظه: %mem=8GB را با RAM موجود تنظیم کنید (معمولاً ۷۵٪ RAM کل).
• پردازنده: %nprocshared=8 را با تعداد core موجود تنظیم کنید.
• خط خالی انتها: فایل Gaussian حتماً باید با یک خط خالی تمام شود، وگرنه خطا می‌دهد.

گام ۴ — اجرای محاسبه در ORCA

روی سیستم محلی (ویندوز یا لینوکس)

# لینوکس/macOS
orca molecule.inp > molecule.out 2>&1 &

# ویندوز (Command Prompt)
orca.exe molecule.inp > molecule.out

روی کلاستر دانشگاه (SLURM)

#!/bin/bash
#SBATCH --job-name=orca_opt
#SBATCH --nodes=1
#SBATCH --ntasks-per-node=8
#SBATCH --mem=32GB
#SBATCH --time=24:00:00
#SBATCH --output=orca_%j.out

module load orca/6.0
$ORCA_DIR/orca molecule.inp > molecule.out

پایش پیشرفت محاسبه ORCA

ORCA در حین اجرا اطلاعات را به فایل .out می‌نویسد. می‌توانید با دستور زیر پیشرفت را دنبال کنید:

tail -f molecule.out

جملاتی مثل GEOMETRY OPTIMIZATION CYCLE نشان‌دهنده پیشرفت بهینه‌سازی و HURRAY - OPTIMIZATION CONVERGED نشان‌دهنده موفقیت است.

گام ۵ — اجرای محاسبه در Gaussian

روی سیستم محلی

# لینوکس
g16 molecule.gjf molecule.log

# ویندوز (از طریق GaussView یا خط فرمان)
g16.exe molecule.gjf

روی کلاستر دانشگاه (SLURM)

<code">#!/bin/bash
#SBATCH --job-name=gaussian_opt
#SBATCH --nodes=1
#SBATCH --ntasks-per-node=8
#SBATCH --mem=32GB
#SBATCH --time=48:00:00

module load gaussian/16
g16 molecule.gjf > molecule.log

بررسی موفقیت محاسبه Gaussian

در انتهای فایل .log باید جمله Normal termination of Gaussian دیده شود. اگر این جمله نبود، خطایی رخ داده — دنبال کلمه Error در فایل log بگردید.

گام ۶ — بازگشت نتایج به ChemDoodle 3D برای تجسم

پس از اتمام محاسبه، فایل‌های خروجی را در ChemDoodle 3D باز کنید:

1. در ChemDoodle 3D از منوی File > Open استفاده کنید.
2. فایل molecule.out (ORCA) یا molecule.log (Gaussian) را انتخاب کنید.
3. ChemDoodle 3D به صورت خودکار هندسه بهینه‌شده آخرین cycle را بارگذاری می‌کند.

قابلیت‌های تجسم در ChemDoodle 3D

• هندسه بهینه‌شده: مستقیم نمایش داده می‌شود. می‌توانید bond length و bond angle را اندازه‌گیری کنید.
• مدهای ارتعاشی (Frequency): از منوی View > Animation می‌توانید هر mode ارتعاشی را به صورت متحرک ببینید.
• سطح مولکولی: از منوی Structure > Generate Surface سطح Van der Waals، SAS، یا SES را رندر کنید.
• MO و electron density: اگر فایل .cube از Gaussian یا ORCA داشته باشید، می‌توانید orbital‌ها را روی سطح مولکول نقشه‌برداری کنید.

ذخیره تصویر برای مقاله

از منوی File > Export Image می‌توانید تصویر با کیفیت بالا (PNG یا SVG) برای استفاده در مقاله خروجی بگیرید. رزولوشن را حداقل روی ۳۰۰ DPI تنظیم کنید.

مثال کاربردی: بهینه‌سازی هندسه اتانول

یک مثال کامل از ابتدا تا انتها برای بهینه‌سازی هندسه اتانول (CH₃CH₂OH):

1. در ChemDoodle 2D اتانول را رسم کنید و با Send to ChemDoodle 3D انتقال دهید.
2. Clean 3D Structure را اجرا کنید.
3. فایل را به فرمت ORCA Input ذخیره کنید (ethanol.inp).
4. فایل را ویرایش کنید:

! B3LYP def2-SVP Opt Freq

%maxcore 2000
%pal nprocs 4 end

* xyz 0 1
C  -0.748833  -0.021715   0.005547
C   0.748833   0.021715  -0.005547
O   1.178920   1.170000   0.598000
H  -1.165660   0.985320   0.013720
H  -1.104930  -0.560400   0.888580
H  -1.104930  -0.560400  -0.877490
H   1.165660  -0.985320  -0.013720
H   1.104930   0.560400  -0.888580
H   1.852000   1.350000   0.010000
*

5. ORCA را اجرا کنید: orca ethanol.inp > ethanol.out
6. پس از اتمام (حدود ۲-۵ دقیقه)، ethanol.out را در ChemDoodle 3D باز کنید.
7. هندسه بهینه‌شده را مشاهده کنید — طول پیوند C-O و زاویه C-O-H را با منوی Measurement اندازه بگیرید.

مثال کاربردی: پیش‌بینی طیف ¹H NMR

برای محاسبه NMR با Gaussian، از همان ساختار بهینه‌شده استفاده کنید:

1. هندسه بهینه‌شده از مثال قبل را در ChemDoodle 3D باز کنید.
2. فایل را به فرمت Gaussian Input ذخیره کنید.
3. route section را ویرایش کنید:

%chk=ethanol_nmr.chk
%mem=4GB
%nprocshared=4
#p B3LYP/6-311+G(2d,p) NMR=GIAO

NMR calculation for ethanol

0 1
C  -0.748833  -0.021715   0.005547
...

4. Gaussian را اجرا کنید.
5. در فایل log، بخش SCF GIAO Magnetic shielding tensor را پیدا کنید — مقادیر isotropic shielding برای هر اتم آمده است.
6. برای تبدیل به chemical shift (δ)، مقدار shielding را از shielding TMS (که باید جداگانه محاسبه کنید: σ_TMS ≈ 31.7 ppm با همین سطح نظری) کم کنید.
7. ChemDoodle 2D هم ابزار پیش‌بینی NMR داخلی دارد که برای بررسی سریع‌تر مفید است.

💬 برای دریافت قیمت و مشاوره رایگان → همین الان در تلگرام پیام دهید — پاسخ در کمتر از چند ساعت

دانلود و لایسنس

ChemDoodle 3D

ChemDoodle 3D بخشی از بسته ChemDoodle است. لایسنس فردی از ۲۹ دلار شروع می‌شود و شامل ۱۴ روز آزمایش رایگان است. لایسنس سایتی برای دانشگاه‌ها و مراکز تحقیقاتی با قیمت ویژه در دسترس است.

ORCA

ORCA برای استفاده آکادمیک (غیرتجاری) رایگان است. از طریق سایت رسمی با ثبت‌نام دانلود کنید. برای استفاده تجاری، لایسنس تجاری ORCA از طریق FAccTs GmbH قابل تهیه است.

Gaussian

Gaussian نرم‌افزار تجاری است. لایسنس آکادمیک سالانه و لایسنس دائم هر دو موجود است. برای خرید لایسنس Gaussian برای ایران، با ما تماس بگیرید.

💬 برای دریافت قیمت و مشاوره رایگان → همین الان در تلگرام پیام دهید — پاسخ در کمتر از چند ساعت

سوالات متداول

آیا ORCA واقعاً رایگان است؟

بله، ORCA برای استفاده آکادمیک و تحقیقاتی غیرتجاری کاملاً رایگان است. فقط نیاز به ثبت‌نام در سایت رسمی دارید. برای استفاده در شرکت‌های تجاری باید لایسنس تجاری تهیه شود.

تفاوت ORCA و Gaussian در چیست؟ کدام را انتخاب کنم؟

هر دو پکیج‌های قدرتمند شیمی کوانتومی هستند. ORCA در روش‌های DFT مدرن، محاسبات coupled cluster، و طیف‌سنجی EPR بسیار قوی است و رایگان بودن آن برای محیط‌های آکادمیک مزیت بزرگی است. Gaussian سابقه طولانی‌تری دارد و در برخی مجلات و گروه‌های تحقیقاتی استاندارد است. اگر استاد راهنما یا همکارانتان Gaussian استفاده می‌کنند، Gaussian بهتر است. در غیر اینصورت ORCA گزینه اول است.

آیا ChemDoodle 3D جایگزین GaussView می‌شود؟

برای بسیاری از کاربردها بله. ChemDoodle 3D می‌تواند فایل‌های ورودی و خروجی Gaussian را بخواند و بنویسد، هندسه بهینه‌شده را نمایش دهد، و مدهای ارتعاشی را متحرک کند. GaussView قابلیت‌های بیشتری برای کار اختصاصی با Gaussian دارد (مثل ویرایش route section با GUI)، اما برای اکثر کاربردهای روزمره ChemDoodle 3D کافی است و قیمت آن بسیار پایین‌تر است.

چطور از ChemDoodle 3D فایل ورودی ORCA بگیرم؟

در نسخه ۷.۷ این قابلیت اضافه شده است. از منوی File > Save As فرمت ORCA Input را انتخاب کنید. اگر این گزینه را نمی‌بینید، نسخه خود را به‌روز کنید.

محاسبه‌ام خطا می‌دهد. مشکل از کجاست؟

شایع‌ترین دلایل: (۱) هندسه اولیه بد — ابتدا با force field در ChemDoodle 3D بهینه کنید. (۲) بار یا multiplicity اشتباه — مطمئن شوید بار مولکول و تعداد الکترون‌های جفت‌نشده درست است. (۳) حافظه ناکافی — مقدار %mem یا %maxcore را کاهش دهید. (۴) basis set برای عناصر سنگین — از basis set های مناسب مثل def2-TZVP استفاده کنید.

آیا می‌توانم این نرم‌افزارها را روی کلاستر دانشگاه اجرا کنم؟

بله. ORCA و Gaussian هر دو نسخه‌های MPI برای اجرای موازی دارند. ChemDoodle 3D روی لپ‌تاپ شما کار می‌کند و فایل‌های ورودی/خروجی را با کلاستر از طریق SSH/SFTP تبادل می‌کنید.

چطور لایسنس ChemDoodle یا Gaussian تهیه کنم؟

💬 برای دریافت قیمت و مشاوره رایگان → همین الان در تلگرام پیام دهید — پاسخ در کمتر از چند ساعت