ریاضیات

طیف چرخشی رامان و فاصله خطوط طیفی (4B)

توسط علی رضا نقش نیلچی | یکشنبه سی و یکم خرداد ۱۴۰۵ | 0:2

طیف چرخشی رامان و فاصله خطوط 4B

در طیف رامان، علاوه بر گذارهای ارتعاشی، گذارهای چرخشی نیز قابل مشاهده هستند.

این طیف را Rotational Raman Spectrum می‌نامند و قوانین گزینش آن با طیف چرخشی مایکروویو متفاوت است.

۱. قانون گزینش در رامان چرخشی

در طیف چرخشی مایکروویو داشتیم:

ΔJ = ±1

اما در طیف رامان قانون گزینش به صورت زیر است:

ΔJ = ±2

بنابراین دو نوع خط مشاهده می‌شود:

S branch : ΔJ = +2

O branch : ΔJ = −2

در وسط طیف نیز یک خط قوی مربوط به پراکندگی رایلی قرار دارد.

۲. انرژی ترازهای چرخشی

انرژی چرخشی یک مولکول دو اتمی از رابطه زیر به دست می‌آید:

F(J) = B J(J+1)

که در آن B ثابت چرخشی مولکول است.

۳. محاسبه مکان خطوط رامان

برای گذار J → J+2 اختلاف انرژی چنین خواهد بود:

ΔF = F(J+2) − F(J)

F(J+2) = B (J+2)(J+3)

F(J) = B J(J+1)

ΔF = B[(J+2)(J+3) − J(J+1)]

(J+2)(J+3) = J² +5J +6

J(J+1) = J² + J

ΔF = B(4J + 6)

ΔF = 2B(2J + 3)

۴. موقعیت چند خط اول

اگر J=0 باشد:

Δν = 6B

برای J=1:

Δν = 10B

برای J=2:

Δν = 14B

بنابراین خطوط طیفی چنین خواهند بود:

6B , 10B , 14B , 18B , ...

۵. فاصله خطوط طیفی

اختلاف بین خطوط متوالی:

10B − 6B = 4B

14B − 10B = 4B

18B − 14B = 4B

پس فاصله بین خطوط در طیف چرخشی رامان برابر است با:

Δν = 4B

۶. نتیجه مهم

در طیف چرخشی مایکروویو فاصله خطوط = 2B

در طیف چرخشی رامان فاصله خطوط = 4B

بنابراین با اندازه‌گیری فاصله خطوط رامان می‌توان ثابت چرخشی B و سپس ممان اینرسی و طول پیوند مولکول را محاسبه کرد.

طیف‌سنجی رامان (Raman Spectroscopy)

توسط علی رضا نقش نیلچی | شنبه سی ام خرداد ۱۴۰۵ | 23:58

طیف‌سنجی رامان (Raman Spectroscopy)

در طیف‌سنجی رامان، نور تک‌رنگ (معمولاً لیزر) به مولکول تابیده می‌شود و بخشی از نور پس از برخورد با مولکول پراکنده می‌شود.

بیشتر نور بدون تغییر انرژی پراکنده می‌شود، اما مقدار کمی از آن با تغییر انرژی ظاهر می‌شود که اطلاعات ساختاری مهمی از مولکول می‌دهد.

۱. انواع پراکندگی نور

وقتی نور به مولکول برخورد می‌کند سه حالت ممکن است رخ دهد:

Rayleigh Scattering: انرژی نور تغییر نمی‌کند (رایج‌ترین حالت).

Stokes Raman: مولکول انرژی می‌گیرد و نور با انرژی کمتر پراکنده می‌شود.

Anti‑Stokes Raman: مولکول انرژی از دست می‌دهد و نور با انرژی بیشتر پراکنده می‌شود.

ν_scattered = ν_incident ± ν_vibration

۲. شرط فعال بودن در رامان

برای اینکه یک ارتعاش در طیف رامان دیده شود باید در طول ارتعاش، قطبش‌پذیری مولکول (polarizability) تغییر کند.

این شرط با شرط فعال بودن در IR متفاوت است.

IR active: تغییر ممان دوقطبی

Raman active: تغییر قطبش‌پذیری

۳. مثال مهم

مولکول‌های متقارن مانند:

H₂ , N₂ , O₂

در طیف IR دیده نمی‌شوند چون ممان دوقطبی ندارند.

اما در طیف رامان کاملاً قابل مشاهده‌اند زیرا قطبش‌پذیری آنها در طول ارتعاش تغییر می‌کند.

۴. تفاوت خطوط Stokes و Anti‑Stokes

در دمای اتاق، بیشتر مولکول‌ها در تراز پایه ارتعاشی هستند. بنابراین:

خطوط Stokes قوی‌تر هستند.

خطوط Anti‑Stokes ضعیف‌تر هستند.

شدت این دو نوع خط با توزیع بولتزمن مرتبط است.

۵. شیفت رامان

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

Δν = ν_incident − ν_scattered

این مقدار را Raman Shift می‌نامند و معمولاً بر حسب cm^-1 بیان می‌شود.

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

بررسی حل تحلیلی معادله شرودینگر برای پتانسیل خطی

توسط علی رضا نقش نیلچی | شنبه سی ام خرداد ۱۴۰۵ | 23:50

:::writing
حل قسمت (ج)

بررسی حل تحلیلی معادله شرودینگر برای پتانسیل خطی

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

1) پتانسیل مسئله

پتانسیل خطی به صورت زیر داده می‌شود:

V(x) = kx

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

2) معادله شرودینگر مستقل از زمان

معادله شرودینگر یک بعدی:

-(h^2 / 2m) * d^2ψ/dx^2 + V(x) ψ = E ψ

با قرار دادن پتانسیل:

-(h^2 / 2m) * d^2ψ/dx^2 + kx ψ = E ψ

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

3) مرتب کردن معادله

انتقال جمله انرژی:

-(h^2 / 2m) d^2ψ/dx^2 + (kx - E) ψ = 0

ضرب در (2m / h^2):

d^2ψ/dx^2 = (2m / h^2) (kx - E) ψ

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

4) تغییر متغیر

متغیر جدید تعریف می‌کنیم:

z = (2mk / h^2)^(1/3) ( x - E/k )

با این تغییر متغیر معادله به شکل ساده تبدیل می‌شود:

d^2ψ/dz^2 - z ψ = 0

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

5) معادله Airy

معادله بالا همان معادله Airy است.

حل عمومی آن:

ψ(z) = A Ai(z) + B Bi(z)

که در آن:

Ai(z) تابع Airy نوع اول
Bi(z) تابع Airy نوع دوم

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

6) شرط فیزیکی

تابع Bi(z) برای z بزرگ واگرا می‌شود، بنابراین برای قابل نرمال‌سازی بودن تابع موج باید ضریب آن صفر باشد.

پس:

ψ(z) = A Ai(z)

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

نتیجه

معادله شرودینگر برای پتانسیل خطی دارای حل تحلیلی است، اما جواب آن با توابع ویژه Airy بیان می‌شود:

ψ(x) = A Ai(z)

که در آن:

z = (2mk / h^2)^(1/3) ( x - E/k )
:::

بررسی حل تحلیلی معادله شرودینگر برای پتانسیل خطی

توسط علی رضا نقش نیلچی | شنبه سی ام خرداد ۱۴۰۵ | 23:13

:::writing

حل قسمت (ج)

بررسی حل تحلیلی معادله شرودینگر برای پتانسیل خطی

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

1) پتانسیل مسئله

پتانسیل خطی به صورت زیر داده می‌شود:

V(x) = kx

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

2) معادله شرودینگر مستقل از زمان

معادله شرودینگر یک بعدی:

-(h^2 / 2m) * d^2ψ/dx^2 + V(x) ψ = E ψ

با قرار دادن پتانسیل:

-(h^2 / 2m) * d^2ψ/dx^2 + kx ψ = E ψ

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

3) مرتب کردن معادله

انتقال جمله انرژی:

-(h^2 / 2m) d^2ψ/dx^2 + (kx - E) ψ = 0

ضرب در (2m / h^2):

d^2ψ/dx^2 = (2m / h^2) (kx - E) ψ

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

4) تغییر متغیر

متغیر جدید تعریف می‌کنیم:

z = (2mk / h^2)^(1/3) ( x - E/k )

با این تغییر متغیر معادله به شکل ساده تبدیل می‌شود:

d^2ψ/dz^2 - z ψ = 0

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

5) معادله Airy

معادله بالا همان معادله Airy است.

حل عمومی آن:

ψ(z) = A Ai(z) + B Bi(z)

که در آن:

Ai(z) تابع Airy نوع اول

Bi(z) تابع Airy نوع دوم

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

6) شرط فیزیکی

تابع Bi(z) برای z بزرگ واگرا می‌شود، بنابراین برای قابل نرمال‌سازی بودن تابع موج باید ضریب آن صفر باشد.

پس:

ψ(z) = A Ai(z)

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

نتیجه

معادله شرودینگر برای پتانسیل خطی دارای حل تحلیلی است، اما جواب آن با توابع ویژه Airy بیان می‌شود:

ψ(x) = A Ai(z)

که در آن:

z = (2mk / h^2)^(1/3) ( x - E/k )

:::

کاربرد نظریه اختلال در اتم هلیوم (حالت سینگلت)

توسط علی رضا نقش نیلچی | شنبه سی ام خرداد ۱۴۰۵ | 23:10

:::writing

حل قسمت (ب)

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

V = k |r1 − r2|

می‌خواهیم تصحیح انرژی حالت سینگلت را با استفاده از نظریه اختلال مرتبه اول محاسبه کنیم.

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

1) تابع موج حالت پایه

تابع موج اوربیتال 1s:

phi(r) = (1 / sqrt(pi a0^3)) * exp(-r / a0)

تابع موج فضایی دو الکترون:

psi(r1 , r2) = phi(r1) * phi(r2)

بنابراین:

|psi|^2 = (1 / (pi^2 a0^6)) * exp[-2(r1 + r2)/a0]

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

2) حالت اسپینی سینگلت

تابع اسپین سینگلت:

chi = (1 / sqrt(2)) (↑↓ − ↓↑)

چون تابع اسپین پادمتقارن است، تابع فضایی متقارن خواهد بود.

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

3) تصحیح انرژی مرتبه اول

طبق نظریه اختلال:

E^(1) =

پس:

E^(1) = k ∫ |psi(r1 , r2)|^2 |r1 − r2| d^3r1 d^3r2

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

4) جایگذاری تابع موج

E^(1) = (k / (pi^2 a0^6)) ∫ exp[-2(r1 + r2)/a0] |r1 − r2| d^3r1 d^3r2

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

5) محاسبه انتگرال

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

∫ exp[-2(r1 + r2)/a0] |r1 − r2| d^3r1 d^3r2

= (5/8) pi^2 a0^7

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

6) نتیجه

با قرار دادن نتیجه انتگرال در رابطه انرژی:

E^(1) = (k / (pi^2 a0^6)) * (5/8) pi^2 a0^7

ساده‌سازی:

E^(1) = (5/8) k a0

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

نتیجه نهایی

تصحیح انرژی حالت سینگلت برابر است با:

E^(1) = (5/8) k a0

و انرژی کل سیستم:

E = E0 + (5/8) k a0

:::

برهم‌کنش بین دو الکترون

توسط علی رضا نقش نیلچی | شنبه سی ام خرداد ۱۴۰۵ | 23:3

:::writing

حل قسمت (ب)

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

V = k |r1 - r2|

می‌خواهیم تصحیح انرژی حالت سینگلت را پیدا کنیم.

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

1) تابع موج حالت پایه

تابع موج اوربیتال 1s:

phi(r) = (1 / sqrt(pi * a0^3)) * exp(-r / a0)

تابع موج کل دو الکترون:

psi(r1 , r2) = phi(r1) * phi(r2)

بنابراین:

|psi|^2 = (1 / (pi^2 * a0^6)) * exp(-2(r1 + r2)/a0)

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

2) حالت اسپینی سینگلت

تابع اسپین:

chi = (1 / sqrt(2)) (up down - down up)

در نتیجه تابع فضایی متقارن است.

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

3) تصحیح انرژی مرتبه اول

طبق نظریه اختلال:

E1 =

پس:

E1 = k * integral |psi(r1,r2)|^2 * |r1 - r2| d^3r1 d^3r2

با جایگذاری تابع موج:

E1 = (k / (pi^2 * a0^6)) *

integral exp(-2(r1+r2)/a0) * |r1 - r2| d^3r1 d^3r2

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

4) نتیجه انتگرال

نتیجه انتگرال استاندارد:

integral exp(-2(r1+r2)/a0) |r1-r2| d^3r1 d^3r2

= (5/8) * pi^2 * a0^7

بنابراین:

E1 = (5/8) * k * a0

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

نتیجه نهایی

E1 = (5/8) k a0

و انرژی کل سیستم:

E = E0 + (5/8) k a0

:::

اگر بخواهید، می‌توانم یک **نسخه حتی سازگارتر با بلاگفا (مثل جزوه دانشگاهی)** هم بنویسم که در آن:

- کسرها با `/` نوشته شوند

- توان‌ها با `^` نوشته شوند

- و دقیقاً مثل **حل امتحانی مکانیک کوانتومی** مرتب شود.

مسئله برهم‌کنش بین دو الکترون در اتم هلیوم

توسط علی رضا نقش نیلچی | شنبه سی ام خرداد ۱۴۰۵ | 23:0

:::writing

حل قسمت (ب)

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

V = k | r1 - r2 |

می‌خواهیم تصحیح انرژی حالت سینگلت را با استفاده از نظریه اختلال مرتبه اول محاسبه کنیم.

1) تابع موج حالت پایه

تابع موج اوربیتال 1s اتم هیدروژن به صورت زیر است:

phi_1s(r) = ( 1 / sqrt(pi a0^3) ) exp( -r / a0 )

تابع موج فضایی دو الکترون:

psi(r1 , r2) = phi_1s(r1) phi_1s(r2)

بنابراین:

|psi|^2 = ( 1 / ( pi^2 a0^6 ) ) exp( -2 ( r1 + r2 ) / a0 )

2) حالت اسپینی سینگلت

تابع اسپین سینگلت به صورت زیر است:

chi = (1 / sqrt(2)) ( ↑↓ - ↓↑ )

در نتیجه تابع فضایی متقارن باقی می‌ماند.

3) تصحیح انرژی مرتبه اول

طبق نظریه اختلال:

E1 = < psi | V | psi >

پس داریم:

E1 = k ∫ |psi(r1 , r2)|^2 | r1 - r2 | d^3 r1 d^3 r2

با جایگذاری تابع موج:

E1 = ( k / ( pi^2 a0^6 ) )

∫ exp( -2 ( r1 + r2 ) / a0 ) | r1 - r2 | d^3 r1 d^3 r2

4) محاسبه انتگرال

نتیجه انتگرال استاندارد:

∫ exp( -2 ( r1 + r2 ) / a0 ) | r1 - r2 | d^3 r1 d^3 r2

= ( 5 / 8 ) pi^2 a0^7

پس انرژی تصحیحی برابر می‌شود با:

E1 = ( 5 / 8 ) k a0

نتیجه

تصحیح انرژی حالت سینگلت برابر است با:

E1 = ( 5 / 8 ) k a0

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

E = E0 + ( 5 / 8 ) k a0

:::

اتم هلیوم با برهم‌کنش خطی

توسط علی رضا نقش نیلچی | شنبه سی ام خرداد ۱۴۰۵ | 22:56

در ادامه **حل کامل قسمت (ب) با محاسبه انتگرال**

:::writing

حل قسمت (ب):اتم هلیوم با برهم‌کنش خطی

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

V = k |r₁ - r₂|

می‌خواهیم تصحیح انرژی برای حالت اسپینی سینگلت را به دست آوریم.

1) تابع موج حالت پایه

در تقریب مستقل، هر دو الکترون در اوربیتال \(1s\) قرار دارند:

φ₁s(r) = (1/√(π a₀³)) e-r/a₀

تابع موج فضایی دو الکترون:

ψ(r₁,r₂) = φ₁s(r₁) φ₁s(r₂)

بنابراین:

|ψ|² =

(1 / π² a₀⁶)

e-2(r₁+r₂)/a₀

2) حالت اسپینی سینگلت

تابع اسپینی سینگلت پادمتقارن است:

χ = (1/√2)(↑↓ − ↓↑)

بنابراین تابع فضایی متقارن باقی می‌ماند.

3) تصحیح انرژی مرتبه اول

طبق نظریه اختلال:

E¹ = 〈ψ|V|ψ〉

پس داریم:

E¹ =

k ∫ |ψ(r₁,r₂)|² |r₁ - r₂|

d³r₁ d³r₂

با جایگذاری تابع موج:

E¹ =

(k / π² a₀⁶)

∫ e-2(r₁+r₂)/a₀

|r₁ - r₂|

d³r₁ d³r₂

4) محاسبه انتگرال

با استفاده از نتایج استاندارد انتگرال‌های دو الکترونی برای توابع نمایی، مقدار انتگرال به صورت زیر به دست می‌آید:

∫ e^{-2(r₁+r₂)/a₀} |r₁ - r₂| d³r₁ d³r₂

= (5/8) π² a₀⁷

با جایگذاری در رابطه انرژی:

E¹ = (5/8) k a₀

نتیجه نهایی

تصحیح انرژی حالت سینگلت برابر است با:

E¹ = (5/8) k a₀

بنابراین انرژی حالت پایه هلیوم در حضور این برهم‌کنش خطی برابر خواهد بود با:

E = E₀ + (5/8) k a₀

:::

مثال چهارم: روش وردشی برای نوسانگر آنهارمونیک (پتانسیل x⁴)

توسط علی رضا نقش نیلچی | شنبه سی ام خرداد ۱۴۰۵ | 22:50

مثال چهارم: روش وردشی برای نوسانگر آنهارمونیک (پتانسیل x⁴)

در این مثال می‌خواهیم انرژی حالت پایه پتانسیل آنهارمونیک زیر را با روش تغییرات تقریب بزنیم:

V(x) = (1/2) m ω² x² + λ x⁴

یک تابع آزمون گاوسی انتخاب می‌کنیم:

ψ(x) = exp(-α x²)

که در آن α پارامتر وردشی است.

1) نرمال‌سازی

N = ∫_{-∞}^{∞} e^{-2αx²} dx

فرمول انتگرال گاوسی:

∫_{-∞}^{∞} e^{-ax²} dx = √(π/a)

پس:

N = √(π / 2α)

2) انرژی جنبشی

اپراتور:

T = - (ħ² / 2m) d²/dx²

نتیجه محاسبه استاندارد برای تابع گاوسی:

<T> = ħ² α / (2m)

3) انرژی پتانسیل هارمونیک

<x²> = 1 / (4α)

پس:

<V_h> = (1/2) m ω² (1 / 4α)

<V_h> = m ω² / (8α)

4) جمله آنهارمونیک

برای گاوسی داریم:

<x⁴> = 3 / (16 α²)

پس:

<V_{anh}> = λ (3 / 16 α²)

5) انرژی کل

E(α) =

ħ² α / (2m)

+ m ω² / (8α)

+ 3λ / (16 α²)

6) کمینه کردن انرژی

شرط مینیمم:

dE/dα = 0

مشتق:

dE/dα =

ħ² / (2m)

- m ω² / (8α²)

- 3λ / (8α³)

= 0

این معادله یک معادله جبری برای α بهینه است که معمولاً عددی حل می‌شود.

نتیجه

روش وردشی نشان می‌دهد که حضور جمله λx⁴ باعث افزایش انرژی حالت پایه نسبت به نوسانگر ساده می‌شود.

اگر λ = 0 باشد، نتیجه دقیق نوسانگر هماهنگ به دست می‌آید:

E₀ = (1/2) ħω

مثال سوم: روش وردشی برای اتم هیدروژن

توسط علی رضا نقش نیلچی | شنبه سی ام خرداد ۱۴۰۵ | 22:46

مثال سوم: روش وردشی برای اتم هیدروژن

در این مثال می‌خواهیم با استفاده از روش تغییرات (Variational Method) انرژی حالت پایه اتم هیدروژن را تقریب بزنیم.

هامیلتونی سیستم برابر است با:

H = - (ħ² / 2m) ∇² - e² / (4π ε₀ r)

یک تابع آزمون ساده انتخاب می‌کنیم:

ψ(r) = e^-βr

که در آن β پارامتر وردشی است.

1) نرمال‌سازی

N = ∫ |ψ|² d³r

در مختصات کروی:

d³r = 4π r² dr

پس داریم:

N = 4π ∫₀∞ r² e^-2βr dr

از رابطه انتگرال استفاده می‌کنیم:

∫₀∞ r² e^-ar dr = 2! / a³

در نتیجه:

N = 4π × 2 / (2β)³

N = π / β³

2) انرژی جنبشی

T = - (ħ² / 2m) ∇²

نتیجه محاسبه مقدار چشم‌داشتی انرژی جنبشی:

<T> = ħ² β² / (2m)

3) انرژی پتانسیل

پتانسیل کولنی:

V(r) = - e² / (4π ε₀ r)

انرژی پتانسیل متوسط:

<V> =

[ 4π ∫₀∞ r² e^-2βr (-e²/(4π ε₀ r)) dr ] / N

پس از ساده‌سازی:

<V> = - e² β / (4π ε₀)

4) انرژی کل

E(β) = <T> + <V>

E(β) = ħ² β² / (2m) - e² β / (4π ε₀)

5) کمینه کردن انرژی

dE/dβ = 0

مشتق انرژی:

dE/dβ = ħ² β / m - e² / (4π ε₀)

پس مقدار بهینه پارامتر:

β = m e² / (4π ε₀ ħ²)

6) انرژی حالت پایه

با جایگذاری مقدار β در انرژی:

E₀ = - m e⁴ / [2 (4π ε₀)² ħ²]

که مقدار عددی آن برابر است با:

E₀ = -13.6 eV

بنابراین روش وردشی با این تابع آزمون توانست انرژی حالت پایه اتم هیدروژن را به دست آورد.

مثال دوم: روش وردشی برای چاه پتانسیل بی‌نهایت یک‌بعدی

توسط علی رضا نقش نیلچی | شنبه سی ام خرداد ۱۴۰۵ | 22:43

مثال دوم: روش وردشی برای چاه پتانسیل بی‌نهایت یک‌بعدی

در این مثال می‌خواهیم انرژی حالت پایه یک ذره در چاه پتانسیل بی‌نهایت را با روش تغییرات تقریب بزنیم.

پتانسیل سیستم:

V(x) = 0 , 0 < x < L

V(x) = ∞ , x ≤ 0 or x ≥ L

یک تابع آزمون انتخاب می‌کنیم که شرایط مرزی را ارضا کند:

ψ(x) = x(L-x)

1) نرمال‌سازی

ابتدا باید مقدار زیر را حساب کنیم:

N = ∫₀ᴸ |ψ(x)|² dx

|ψ|² = x² (L-x)²

پس داریم:

N = ∫₀ᴸ x² (L-x)² dx

گسترش عبارت:

(L-x)² = L² -2Lx + x²

x²(L-x)² = L²x² -2Lx³ + x⁴

انتگرال:

∫₀ᴸ L²x² dx = L⁵/3

∫₀ᴸ 2Lx³ dx = L⁵/2

∫₀ᴸ x⁴ dx = L⁵/5

پس:

N = L⁵ ( 1/3 − 1/2 + 1/5 )

N = L⁵ / 30

2) انرژی جنبشی

در داخل چاه پتانسیل صفر است، بنابراین انرژی فقط جنبشی است.

H = - (ħ² / 2m) d²/dx²

مشتق اول تابع موج:

ψ'(x) = L - 2x

مشتق دوم:

ψ''(x) = -2

انرژی وردشی:

E = ∫ ψ* ( -ħ²/2m ) ψ'' dx / N

جایگذاری:

E = (ħ²/m) ∫₀ᴸ ψ(x) dx / N

انتگرال:

∫₀ᴸ x(L-x) dx

= ∫₀ᴸ (Lx - x²) dx

= L³/6

پس:

E = (ħ²/m) (L³/6) / (L⁵/30)

3) انرژی نهایی

E = 5 ħ² / (m L²)

انرژی دقیق حالت پایه این سیستم برابر است با:

E₁ = π² ħ² / (2mL²)

که مقدار به‌دست‌آمده با روش وردشی یک تقریب برای آن است.

مثال: الگوریتم محاسبه انرژی ذره در جعبه یک‌بعدی

توسط علی رضا نقش نیلچی | شنبه سی ام خرداد ۱۴۰۵ | 22:37

مثال: الگوریتم محاسبه انرژی ذره در جعبه یک‌بعدی

در این مثال می‌خواهیم انرژی‌های مجاز یک ذره با جرم μ را که در یک جعبه یک‌بعدی با طول L قرار دارد محاسبه کنیم.

مرحله 1: تعریف ناحیه حرکت ذره

0 < x < L

پتانسیل جعبه به صورت زیر تعریف می‌شود:

V(x) = 0 for 0 < x < L

V(x) = ∞ otherwise

در نتیجه تابع موج روی دیواره‌های جعبه صفر است:

ψ(0) = 0

ψ(L) = 0

مرحله 2: نوشتن معادله شرودینگر

-(ħ² / 2μ) d²ψ/dx² = Eψ

مرحله 3: حل معادله دیفرانسیل

حل عمومی معادله به صورت زیر است:

ψ(x) = A sin(kx) + B cos(kx)

با اعمال شرط مرزی:

ψ(0) = 0

نتیجه می‌گیریم:

B = 0

پس تابع موج به شکل زیر در می‌آید:

ψ(x) = A sin(kx)

مرحله 4: اعمال شرط مرزی دوم

ψ(L) = 0

بنابراین:

sin(kL) = 0

که نتیجه می‌دهد:

kL = nπ

k = nπ / L

که در آن:

n = 1 , 2 , 3 , ...

مرحله 5: محاسبه انرژی‌های مجاز

E = ħ²k² / 2μ

با قرار دادن مقدار k:

E_n = (n² π² ħ²) / (2μ L²)

مرحله 6: مثال عددی

اگر:

L = 2

انرژی حالت پایه:

E1 = (π² ħ²) / (8 μ)

انرژی اولین حالت برانگیخته:

E2 = (4 π² ħ²) / (8 μ)

E2 = (π² ħ²) / (2 μ)

روش وردشی برای نوسانگر هماهنگ سه‌بعدی

توسط علی رضا نقش نیلچی | شنبه سی ام خرداد ۱۴۰۵ | 22:31

حل بخش (الف): روش وردشی برای نوسانگر هماهنگ سه‌بعدی

در این مسئله می‌خواهیم با استفاده از روش تغییرات (Variational Method) انرژی حالت پایه نوسانگر هماهنگ سه‌بعدی را به دست آوریم.

تابع آزمون داده شده است:

ψ(r) = m rⁿ e^-βr

هامیلتونی نوسانگر هماهنگ سه‌بعدی:

H = - (ħ² / 2μ) ∇² + (1/2) μ ω² r²

انرژی وردشی به صورت زیر تعریف می‌شود:

E(β) = <ψ|H|ψ> / <ψ|ψ>

1) محاسبه نرمال‌سازی

ابتدا مقدار زیر را محاسبه می‌کنیم:

<ψ|ψ> = ∫ |ψ|² d³r

در مختصات کروی داریم:

d³r = 4π r² dr

بنابراین:

|ψ|² = m² r²ⁿ e^-2βr

N = 4π m² ∫₀∞ r²ⁿ⁺² e^-2βr dr

تغییر متغیر:

y = 2βr

dr = dy / (2β)

پس انتگرال به صورت زیر می‌شود:

N = 4π m² /(2β)²ⁿ⁺³ ∫₀∞ y²ⁿ⁺² e^-y dy

با استفاده از رابطه داده شده:

∫₀∞ y^k e^-y dy = k!

در نتیجه:

N = 4π m² (2n+2)! / (2β)²ⁿ⁺³

2) انرژی پتانسیل

V(r) = (1/2) μ ω² r²

مقدار چشم‌داشتی انرژی پتانسیل:

<V> = ( ∫ ψ* V ψ d³r ) / N

<V> = (1/2) μ ω² 4π m² ∫₀∞ r²ⁿ⁺⁴ e^-2βr dr / N

انتگرال برابر است با:

∫₀∞ r²ⁿ⁺⁴ e^-2βr dr =

(2n+4)! / (2β)²ⁿ⁺⁵

پس بعد از ساده‌سازی:

<V> = μ ω² (2n+4)(2n+3) / (8 β²)

3) انرژی جنبشی

اپراتور انرژی جنبشی:

T = - (ħ² / 2μ) ∇²

برای تابع شعاعی:

∇²ψ = (1/r²) d/dr ( r² dψ/dr )

با محاسبه مشتق‌ها و انجام انتگرال‌ها نتیجه به دست می‌آید:

<T> = (ħ² β²) / (2μ)

4) انرژی کل

E(β) = <T> + <V>

E(β) = (ħ² β²)/(2μ) + μ ω² (2n+4)(2n+3)/(8 β²)

5) کمینه کردن انرژی

dE/dβ = 0

(ħ²/μ) β − μ ω² (2n+4)(2n+3)/(4 β³) = 0

پس داریم:

β⁴ = μ² ω² (2n+4)(2n+3) / (4 ħ²)

6) انرژی مینیمم

E_min = (ħω / 2) √((2n+4)(2n+3))

این مقدار یک کران بالا برای انرژی حالت پایه نوسانگر هماهنگ سه‌بعدی است که از روش وردشی به دست آمده است.

الگوریتم محاسبه مسئله ذره در جعبه سه‌بعدی

توسط علی رضا نقش نیلچی | شنبه سی ام خرداد ۱۴۰۵ | 22:19

الگوریتم محاسبه مسئله ذره در جعبه سه‌بعدی

مرحله 1: دریافت داده‌ها

رقم یکان شماره دانشجویی را n و رقم دهگان را m در نظر می‌گیریم.

اگر هرکدام صفر باشند مقدار آن‌ها برابر 1 قرار داده می‌شود.

مرحله 2: تعیین پارامترهای مسئله

طول جعبه:

L = n

میدان الکتریکی:

E = m x̂

برای این مسئله:

n = 3

m = 4

L = 3

E = 4 x̂

مرحله 3: محاسبه پتانسیل الکتریکی

رابطه میدان و پتانسیل:

E = -∇φ

در راستای x داریم:

Ex = - dφ/dx = 4

بنابراین:

dφ/dx = -4

با انتگرال‌گیری:

φ(x) = -4x

مرحله 4: محاسبه انرژی پتانسیل ذره

V = qφ

پس:

V = -4qx

مرحله 5: نوشتن هامیلتونی سیستم

H = (1 / 2μ)(-iħ∇ − qA)² + qφ + Vbox

مرحله 6: معادله شرودینگر داخل جعبه

[(1 / 2μ)(-iħ∇ − qA)² − 4qx] ψ(x,y,z) = E ψ(x,y,z)

مرحله 7: انرژی‌های ذره در جعبه سه‌بعدی

E(nx,ny,nz) = (ħ²π² / 2μL²)(nx² + ny² + nz²)

با قرار دادن L = 3 :

E(nx,ny,nz) = (ħ²π² / 18μ)(nx² + ny² + nz²)

مرحله 8: حالت پایه

(nx,ny,nz) = (1,1,1)

E111 = ħ²π² / 6μ

مرحله 9: اولین حالت‌های برانگیخته

(2,1,1) , (1,2,1) , (1,1,2)

E = ħ²π² / 3μ

این سه حالت دارای انرژی یکسان هستند و بنابراین تبهگنی سه‌گانه دارند.

مرحله 10: تصحیح مرتبه اول انرژی

ΔE(1) = 〈n| -4qx |n〉

چون:

<x> = L / 2

و L = 3 :

<x> = 3 / 2

ΔE(1) = -4q × (3/2)

ΔE(1) = -6q

مرحله 11: دامنه گذار در تقریب دوقطبی

V = -d · E

d = q r

V = -4qx

حالت اولیه:

|111>

اولین حالت‌های برانگیخته:

|211> , |121> , |112>

تنها گذار مجاز:

|111> → |211>

المان ماتریسی:

<211|x|111> = -16 / (3π²)

دامنه گذار:

V211,111 = -4q <211|x|111>

V211,111 = 64q / (3π²)

حل مسئله ذره در جعبه سه‌بعدی در حضور میدان الکتریکی و مغناطیسی

توسط علی رضا نقش نیلچی | شنبه سی ام خرداد ۱۴۰۵ | 21:56

حل مسئله ذره در جعبه سه‌بعدی در حضور میدان الکتریکی و مغناطیسی

در این مسئله داده شده است:

m = 4 , n = 3

بنابراین طول ضلع جعبه برابر است با:

L = n = 3

و میدان الکتریکی ثابت در راستای محور x برابر است با:

E = 4 x-hat

فرض می‌کنیم ذره دارای جرم mu و بار q باشد.

الف) معادله شرودینگر ذره

ذره در یک جعبه مکعبی سه‌بعدی با طول ضلع 3 قرار دارد. بنابراین ناحیه مجاز حرکت ذره چنین است:

0 < x < 3

0 < y < 3

0 < z < 3

پتانسیل جعبه به صورت زیر است:

V_box(x,y,z) = 0 for 0 < x < 3 , 0 < y < 3 , 0 < z < 3

V_box(x,y,z) = infinity otherwise

یعنی داخل جعبه پتانسیل صفر است و خارج جعبه بی‌نهایت است. بنابراین تابع موج روی دیواره‌های جعبه صفر می‌شود:

psi = 0

at x = 0 , x = 3 , y = 0 , y = 3 , z = 0 , z = 3

در حضور میدان‌های الکترومغناطیسی، هامیلتونی ذره باردار به صورت زیر نوشته می‌شود:

H = (1 / 2mu) (-i hbar grad - q A)^2 + q phi + V_box

چون میدان الکتریکی در راستای x است و مقدار آن 4 می‌باشد:

E = 4 x-hat

از رابطه میدان الکتریکی و پتانسیل الکتریکی داریم:

E = - grad phi

پس:

E_x = - d phi / dx = 4

در نتیجه:

phi(x) = -4x

بنابراین انرژی پتانسیل الکتریکی برابر است با:

V_E = q phi = -4qx

برای میدان مغناطیسی ثابت می‌توان پتانسیل برداری را به شکل زیر انتخاب کرد:

A = (1/2) B × r

بنابراین معادله شرودینگر مستقل از زمان درون جعبه به صورت زیر است:

[(1 / 2mu) (-i hbar grad - q A)^2 - 4qx] psi(x,y,z) = E psi(x,y,z)

این معادله برای ناحیه زیر معتبر است:

0 < x < 3 , 0 < y < 3 , 0 < z < 3

و تابع موج روی دیواره‌های جعبه صفر است.

ب) بررسی امکان انتخاب میدان مغناطیسی مناسب

ابتدا حالت بدون میدان الکترومغناطیسی را بررسی می‌کنیم. اگر میدان الکتریکی و مغناطیسی وجود نداشته باشند، هامیلتونی ذره در جعبه سه‌بعدی برابر است با:

H_0 = - (hbar^2 / 2mu) grad^2

تابع موج‌های ذره در جعبه سه‌بعدی مکعبی به صورت زیر هستند:

psi_nx,ny,nz(x,y,z) =
(2/3)^(3/2)
sin(nx pi x / 3)
sin(ny pi y / 3)
sin(nz pi z / 3)

که در آن:

nx , ny , nz = 1 , 2 , 3 , ...

ترازهای انرژی در حالت بدون میدان برابرند با:

E_nx,ny,nz^(0) =
(hbar^2 pi^2 / 18mu) (nx^2 + ny^2 + nz^2)

حالت پایه

حالت پایه متناظر است با:

(nx,ny,nz) = (1,1,1)

پس انرژی حالت پایه:

E_111^(0) = hbar^2 pi^2 / 6mu

اولین حالت برانگیخته

اولین حالت‌های برانگیخته عبارت‌اند از:

(2,1,1) , (1,2,1) , (1,1,2)

برای هر سه حالت داریم:

nx^2 + ny^2 + nz^2 = 6

پس انرژی اولین حالت برانگیخته:

E_first^(0) = hbar^2 pi^2 / 3mu

این تراز انرژی سه‌گانه تبهگن است، زیرا سه حالت مختلف انرژی یکسان دارند.

اثر میدان الکتریکی

میدان الکتریکی یک جمله اختلالی به هامیلتونی اضافه می‌کند:

H_E = -4qx

این جمله فقط به مختصات x وابسته است، بنابراین تقارن بین سه جهت x، y و z را از بین می‌برد.

آیا میدان مغناطیسی می‌تواند اثر میدان الکتریکی را حذف کند؟ پاسخ کلی این است:

خیر. میدان مغناطیسی به طور کلی نمی‌تواند پتانسیل خطی ناشی از میدان الکتریکی را حذف کند.

زیرا میدان الکتریکی از طریق پتانسیل اسکالر وارد می‌شود:

q phi = -4qx

اما میدان مغناطیسی از طریق پتانسیل برداری در عبارت زیر وارد می‌شود:

(-i hbar grad - q A)^2

بنابراین ساختار اثر میدان الکتریکی و مغناطیسی متفاوت است. بهترین انتخاب برای اینکه جواب‌ها تا حد امکان شبیه حالت بدون برهم‌کنش باشند، این است که میدان مغناطیسی را صفر بگیریم:

B = 0

تصحیح انرژی در مرتبه اول

اگر میدان الکتریکی را به عنوان اختلال در نظر بگیریم، تصحیح مرتبه اول انرژی برابر است با:

Delta E^(1) = < nx,ny,nz | -4qx | nx,ny,nz >

در جعبه‌ای از 0 تا L مقدار میانگین x برابر است با:

<x> = L / 2

چون در این مسئله L = 3 است، داریم:

<x> = 3 / 2

پس:

Delta E^(1) = -4q(3/2) = -6q

بنابراین انرژی‌ها تا مرتبه اول تقریباً برابرند با:

E_nx,ny,nz =
(hbar^2 pi^2 / 18mu)(nx^2 + ny^2 + nz^2) - 6q

در مرتبه اول، چون تصحیح انرژی برای همه حالت‌ها یکسان است، تبهگنی اولین حالت برانگیخته همچنان باقی می‌ماند. اما در حل دقیق‌تر یا در مرتبه‌های بالاتر، چون میدان الکتریکی راستای x را از دو راستای دیگر متمایز می‌کند، تبهگنی می‌تواند شکسته شود.

ج) دامنه گذار به اولین حالت برانگیخته در تقریب دوقطبی الکتریکی

در تقریب دوقطبی الکتریکی، برهم‌کنش ذره با میدان الکتریکی به صورت زیر است:

V = - d . E

که در آن:

d = q r

چون میدان در راستای x است:

E = 4 x-hat

پس اختلال دوقطبی برابر می‌شود با:

V = -4qx

حالت اولیه، حالت پایه است:

|i> = |111>

اولین حالت‌های برانگیخته عبارت‌اند از:

|211> , |121> , |112>

چون اختلال فقط به x وابسته است، گذار مستقیم فقط به حالتی مجاز است که عدد کوانتومی مربوط به x تغییر کند. بنابراین گذار اصلی:

|111> → |211>

المان ماتریسی مورد نیاز:

<211|x|111>

توابع موج عبارت‌اند از:

psi_111 =
(2/3)^(3/2)
sin(pi x / 3)
sin(pi y / 3)
sin(pi z / 3)

psi_211 =
(2/3)^(3/2)
sin(2pi x / 3)
sin(pi y / 3)
sin(pi z / 3)

پس:

<211|x|111> =
∫∫∫ psi_211* x psi_111 dx dy dz

با محاسبه انتگرال‌ها نتیجه می‌شود:

<211|x|111> = -16 / (3 pi^2)

بنابراین المان ماتریسی اختلال دوقطبی:

V_211,111 = <211| -4qx |111>

پس:

V_211,111 = -4q <211|x|111>

با جایگذاری داریم:

V_211,111 = -4q (-16 / 3pi^2)

پس دامنه گذار یا المان کوپل‌شدن برابر است با:

V_211,111 = 64q / (3pi^2)

برای دو حالت دیگر داریم:

<121|x|111> = 0

<112|x|111> = 0

زیرا عملگر x فقط روی بخش x تابع موج اثر می‌گذارد و در جهت‌های y و z توابع موج متعامد می‌شوند.

جمع‌بندی نهایی

معادله شرودینگر ذره درون جعبه در حضور میدان‌ها:

[(1 / 2mu) (-i hbar grad - q A)^2 - 4qx] psi = E psi

انرژی‌های ذره در جعبه بدون میدان:

E_nx,ny,nz^(0) =
(hbar^2 pi^2 / 18mu)(nx^2 + ny^2 + nz^2)

تصحیح مرتبه اول ناشی از میدان الکتریکی:

Delta E^(1) = -6q

دامنه گذار از حالت پایه به اولین حالت برانگیخته مجاز:

|111> → |211>

V_211,111 = 64q / (3pi^2)

برهم‌کنش ارتعاش-چرخش و تغییر فاصله خطوط

توسط علی رضا نقش نیلچی | شنبه سی ام خرداد ۱۴۰۵ | 16:38

برهم‌کنش ارتعاش-چرخش و تغییر فاصله خطوط

در مدل‌های ساده فرض می‌شود که مقدار B برای تمام ترازهای ارتعاشی ثابت است. اما در واقعیت، با افزایش ارتعاش، میانگین فاصله اتم‌ها (طول پیوند) افزایش می‌یابد.

۱. وابستگی B به تراز ارتعاشی

ثابت چرخشی در تراز v از رابطه زیر به دست می‌آید:

B_v = B_e − α_e(v + 1/2)

که در آن:

B_e: ثابت چرخشی تعادلی.

α_e: ثابت برهم‌کنش ارتعاش-چرخش (عددی مثبت و کوچک).

نتیجه: ثابت چرخشی در تراز v=1 (یعنی B₁) همیشه کمتر از ثابت چرخشی در تراز v=0 (یعنی B₀) است.

۲. اثر بر شاخه‌های P و R

به دلیل اینکه B₁ < B₀، فاصله بین خطوط طیفی در دو شاخه متفاوت می‌شود:

در شاخه R: با افزایش J، خطوط به هم نزدیک‌تر می‌شوند (تراکم خطوط).

در شاخه P: با افزایش J، خطوط از هم دورتر می‌شوند.

۳. فرمول دقیق‌تر مکان خطوط

با در نظر گرفتن این تفاوت، فرکانس خطوط به صورت زیر محاسبه می‌شود:

ν_R(J) = ν₀ + (B₁ + B₀)(J+1) + (B₁ − B₀)(J+1)²

ν_P(J) = ν₀ − (B₁ + B₀)J + (B₁ − B₀)J²

چون (B₁ − B₀) عددی منفی است، در شاخه R این ترم باعث کاهش فرکانس و در شاخه P باعث افزایش فاصله می‌شود.

مثال مفهومی

اگر در یک طیف تجربی مشاهده کردید که در سمت فرکانس‌های بالا (سمت راست)، فاصله خطوط مدام در حال کم شدن است، این نشانه مستقیم افزایش طول پیوند در حالت برانگیخته است.

نتیجه‌گیری

طیف ارتعاشی-چرخشی نه تنها فرکانس ارتعاش و ثابت چرخشی را به ما می‌دهد، بلکه با بررسی تغییر فواصل، می‌توانیم بفهمیم که پیوند در اثر ارتعاش چقدر "شل" یا "طویل" می‌شود.

*ترکیب حرکت ارتعاشی + چرخشی و تشکیل طیف ارتعاشی–چرخشی

توسط علی رضا نقش نیلچی | شنبه سی ام خرداد ۱۴۰۵ | 16:37

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

**ترکیب حرکت ارتعاشی + چرخشی و تشکیل طیف ارتعاشی–چرخشی (Vibrational–Rotational Spectrum)**

در این مرحله می‌بینیم چرا باندهای IR شامل دو شاخه معروف **P** و **R** هستند و چرا در وسط‌شان یک شکاف (حالت Q) وجود دارد.

کد زیر آمادهٔ استفاده در بلاگفا است:

```html

طیف ارتعاشی–چرخشی: تشکیل شاخه‌های P و R

وقتی یک مولکول همزمان بچرخد و ارتعاش کند، انرژی کل آن برابر است با:

E(v,J) = G(v) + F(J)

که در آن:

G(v): انرژی ارتعاشی (هماهنگ یا ناهماهنگ)

F(J): انرژی چرخشی (B J(J+1))

۱. قانون گزینش در گذارهای IR

در طیف‌سنجی ارتعاشی–چرخشی، قوانین گزینش زیر برقرارند:

Δv = +1 (گذار پایه: v=0 → v=1)

ΔJ = ±1 (نه صفر!)

چرا ΔJ = 0 ممنوع است؟

چون فوتون گشتاور زاویه‌ای دارد و مولکول حتماً باید یک واحد کوانتومی گشتاور جذب یا دفع کند.

نتیجهٔ این قانون:

ΔJ = −1 → شاخه P (خطوط در سمت چپ، انرژی کمتر)

ΔJ = +1 → شاخه R (خطوط در سمت راست، انرژی بیشتر)

بنابراین طیف شامل یک شکاف مرکزی (Q-branch) است که برای مولکول‌های دو اتمی مجاز نیست.

۲. مکان خطوط شاخه P و R

اگر گذار v=0→1 را بررسی کنیم، مکان خطوط از رابطه زیر می‌آید:

ν_P(J) = ν₀ − 2B(J) (شاخه P: ΔJ = −1)

ν_R(J) = ν₀ + 2B(J+1) (شاخه R: ΔJ = +1)

که در آن ν₀ مقدار گذار ارتعاشی خالص است (مثلاً 2880 cm^-1).

۳. مثال عملی: طیف ارتعاشی–چرخشی CO

برای CO داریم:

ν₀ = 2143 cm^-1

B = 1.922 cm^-1

محاسبه چند خط نخست:

شاخه P (ΔJ = −1):

P(1): 2143 − 2(1.922) = 2139.16

P(2): 2143 − 2(2×1.922) = 2135.32

P(3): 2143 − 2(3×1.922) = 2131.48

شاخه R (ΔJ = +1):

R(0): 2143 + 2(1)×1.922 = 2146.84

R(1): 2143 + 2(2)×1.922 = 2150.68

R(2): 2143 + 2(3)×1.922 = 2154.52

نتیجه: طیف شامل دو «بال» (P و R) و یک شکاف مرکزی است.

۴. شدت خطوط

شدت خطوط در هر شاخه به جمعیت اولیه J بستگی دارد:

N_J ∝ (2J+1) exp(−BJ(J+1)/kT)

در دمای اتاق معمولاً بیشترین شدت متعلق به گذارهایی است که از J≈6 یا 7 شروع می‌شوند.

```

---

مرحله بعدی طبق ترتیب آموزش:

### ✅ **اثر ناهماهنگی و اعوجاج چرخشی بر شکل شاخه‌های P و R**

(چرا فاصله خطوط یکنواخت نیست؟ چرا شاخه R سریع‌تر باز می‌شود؟)

ادامه بدهم؟

h3>ناهماهنگی در طیف‌سنجی ارتعاشی: مدل پتانسیل مورس

توسط علی رضا نقش نیلچی | شنبه سی ام خرداد ۱۴۰۵ | 16:35

ناهماهنگی در طیف‌سنجی ارتعاشی: مدل پتانسیل مورس

در مدل نوسانگر هماهنگ، فاصله‌ی بین ترازها دقیقاً ثابت است. اما در دنیای واقعی با افزایش v (عدد کوانتومی ارتعاشی)، مولکول به سمت گسستگی پیوند نزدیک می‌شود و فاصله ترازها کم‌کم کمتر می‌شود. این رفتار را ناهماهنگی (Anharmonicity) می‌نامند.

۱. فرمول انرژی ناهماهنگ

فرمول انرژی ارتعاشی واقعی‌تر به شکل زیر است:

G(v) = ω_e (v + 1/2) − ω_e x_e (v + 1/2)²

که در آن:

ω_e: فرکانس ارتعاشی هماهنگ

x_e: ثابت ناهماهنگی

ω_e x_e: اندازه اثر ناهماهنگی

۲. پیامدهای ناهماهنگی

فاصله‌ی بین ترازها با افزایش v کمتر می‌شود.

فرکانس گذار v=0→1 از v=1→2 کمی بزرگ‌تر است.

امکان مشاهدهٔ گذارهای فراارتعاشی (overtones) مانند v=0→2 وجود دارد.

در نهایت، انرژی کافی باعث شکستن پیوند می‌شود (گسستگی پیوند).

۳. گذارهای ارتعاشی با ناهماهنگی

در نوسانگر هماهنگ فقط گذار Δv = ±1 مجاز بود.

اما در مدل مورس (ناهماهنگ):

Δv = ±1 (قوی‌ترین گذار)

Δv = ±2 (overtone اول)

Δv = ±3 (overtone دوم)

شدت گذارهای overtone بسیار کمتر است، اما در طیف IR قابل مشاهده‌اند.

۴. مثال عددی: محاسبه موقعیت گذار v=0→1

اگر برای مولکولی داشته باشیم:

ω_e = 3000 cm^-1

x_e = 0.02

مکان خط v=0→1 چقدر است؟

از فرمول اختلاف انرژی استفاده می‌کنیم:

ΔG = G(1) − G(0)

محاسبه:

G(1) = ω_e(3/2) − ω_ex_e(3/2)²

G(0) = ω_e(1/2) − ω_ex_e(1/2)²

اختلاف:

ΔG = ω_e − 2ω_ex_e

ΔG = 3000 − 2(3000)(0.02)

ΔG = 3000 − 120 = 2880 cm^-1

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

۵. پتانسیل مورس

پتانسیل مورس شکل دقیق‌تری برای توصیف انرژی پیوند است:

V(r) = D_e (1 − e^{−a(r − r_e)})²

ویژگی‌های مهم:

در نزدیکی r_e شبیه یک نوسانگر هماهنگ است.

برای کشیدگی زیاد، انرژی اشباع شده و در نهایت پیوند می‌شکند.

عمق چاه پتانسیل (D_e) انرژی گسستگی مولکول است.

طیف‌سنجی ارتعاشی: مدل نوسانگر هماهنگ (Harmonic Oscillator)

توسط علی رضا نقش نیلچی | شنبه سی ام خرداد ۱۴۰۵ | 16:32

طیف‌سنجی ارتعاشی: مدل نوسانگر هماهنگ (Harmonic Oscillator)

وقتی اتم‌ها در یک مولکول دو اتمی در طول محور پیوند نوسان می‌کنند، انرژی ارتعاشی پیدا می‌کنند. ساده‌ترین مدل برای توصیف این حرکت، مدل نوسانگر هماهنگ است.

۱. فرمول انرژی ترازهای ارتعاشی

انرژی ترازهای ارتعاشی بر حسب عدد موجی (cm^-1) از رابطه زیر به‌دست می‌آید:

G(v) = ω_e (v + 1/2)

در این رابطه:

v: عدد کوانتومی ارتعاشی است که مقادیر 0, 1, 2, ... را می‌پذیرد.

ω_e: فرکانس ارتعاشی تعادلی (بر حسب cm^-1).

۲. انرژی نقطه صفر (Zero Point Energy)

حتی در پایین‌ترین تراز انرژی (v = 0)، مولکول همچنان دارای انرژی ارتعاشی است. این انرژی را "انرژی نقطه صفر" می‌نامند:

G(0) = 1/2 ω_e

این یعنی مولکول هرگز از ارتعاش باز نمی‌ایستد، حتی در دمای صفر مطلق!

۳. رابطه فرکانس ارتعاشی با ثابت فنر (k)

فرکانس ارتعاش به سختی پیوند (ثابت فنر k) و جرم اتم‌ها (جرم کاهش‌یافته μ) بستگی دارد:

ω_e = (1 / 2πc) × √(k / μ)

در این رابطه:

k: ثابت نیروی پیوند (Force Constant) بر حسب N/m.

μ: جرم کاهش یافته بر حسب kg.

c: سرعت نور بر حسب cm/s.

مثال عددی: محاسبه ثابت نیروی پیوند (k)

اگر فرکانس ارتعاشی مولکول HCl برابر با 2886 cm^-1 باشد و جرم کاهش یافته آن μ = 1.62 × 10^-27 kg باشد، ثابت فنر چقدر است؟

ابتدا فرمول را برای k مرتب می‌کنیم:

k = (2πc ω_e)² × μ

k = (2 × 3.14 × 3×10¹⁰ × 2886)² × 1.62×10^-27

k ≈ (5.44 × 10¹⁴)² × 1.62×10^-27

k ≈ 2.96 × 10²⁹ × 1.62×10^-27

k ≈ 480 N/m

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

شدت خطوط طیفی و توزیع جمعیت در ترازهای چرخشی

توسط علی رضا نقش نیلچی | شنبه سی ام خرداد ۱۴۰۵ | 16:31

شدت خطوط طیفی و توزیع جمعیت در ترازهای چرخشی

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

مرحله ۱: عوامل موثر بر جمعیت تراز J

جمعیت هر تراز (N_J) نسبت به تراز پایه (N₀) از رابطه زیر تبعیت می‌کند:

N_J / N₀ = (2J + 1) exp( −E_J / kT )

در این رابطه:

(2J + 1): تبهگنی تراز (تعداد حالت‌های هم‌انرژی). این عامل باعث افزایش شدت با بالا رفتن J می‌شود.

exp(−E_J / kT): فاکتور بولتزمن. این عامل با افزایش انرژی (افزایش J) باعث کاهش جمعیت می‌شود.

مرحله ۲: یافتن تراز با بیشترین جمعیت (J_max)

ترازی که بیشترین جمعیت را دارد و قوی‌ترین خط طیفی را ایجاد می‌کند، از رابطه زیر محاسبه می‌شود:

J_max ≈ √( kT / 2hBc ) − 1/2

مثال عددی

برای مولکولی با ثابت چرخشی B = 2 cm^-1 در دمای اتاق (T = 300 K)، کدام تراز J بیشترین جمعیت را دارد؟

k = 1.38 × 10^-23 J/K (ثابت بولتزمن)

h = 6.626 × 10^-34 J.s

c = 3.00 × 10¹⁰ cm/s (سرعت نور بر حسب سانتی‌متر)

kT / 2hBc = (1.38×10^-23 × 300) / (2 × 6.626×10^-34 × 3.00×10¹⁰ × 2)

kT / 2hBc ≈ 4.14×10^-21 / 7.95×10^-23 ≈ 52

J_max ≈ √52 − 0.5 ≈ 7.2 − 0.5 ≈ 6.7

بنابراین تراز J = 7 بیشترین جمعیت را دارد و گذار از J=7 به J=8 قوی‌ترین خط در طیف خواهد بود.

نتیجه‌گیری

۱. با افزایش دما، J_max به مقادیر بالاتر منتقل می‌شود (خطوط J بالا قوی‌تر می‌شوند).

۲. با افزایش ثابت چرخشی B، J_max کاهش می‌یابد.

محاسبه ثابت اعوجاج گریز از مرکز (D) از داده‌های طیفی

توسط علی رضا نقش نیلچی | شنبه سی ام خرداد ۱۴۰۵ | 16:29

محاسبه ثابت اعوجاج گریز از مرکز (D) از داده‌های طیفی

در روتور غیرصلب انرژی ترازها از رابطه زیر به‌دست می‌آید:

E(J) = B J(J+1) − D J²(J+1)²

مکان خطوط طیفی (گذار J → J+1):

ν(J) = 2B(J+1) − 4D(J+1)³

مثال عددی

فرض کنید از طیف تجربی داریم:

ν(0→1) = 3.96 cm^-1

ν(1→2) = 7.84 cm^-1

می‌خواهیم B و D را محاسبه کنیم.

مرحله ۱: نوشتن معادلات

برای J=0:

ν₀ = 2B(1) − 4D(1)³

3.96 = 2B − 4D

برای J=1:

ν₁ = 2B(2) − 4D(2)³

7.84 = 4B − 32D

مرحله ۲: حل دستگاه معادلات

معادله اول:

2B − 4D = 3.96

معادله دوم:

4B − 32D = 7.84

معادله اول را در 2 ضرب می‌کنیم:

4B − 8D = 7.92

تفریق از معادله دوم:

(4B − 32D) − (4B − 8D) = 7.84 − 7.92

−24D = −0.08

D = 0.0033 cm^-1

مرحله ۳: محاسبه B

2B − 4(0.0033) = 3.96

2B − 0.0132 = 3.96

2B = 3.9732

B = 1.9866 cm^-1

نتیجه نهایی

B ≈ 1.99 cm^-1

D ≈ 0.0033 cm^-1

می‌بینیم که مقدار D بسیار کوچک‌تر از B است،

اما در J های بزرگ اثر آن کاملاً قابل مشاهده می‌شود.

اثر گریز از مرکز در طیف چرخشی (Non-Rigid Rotor)

توسط علی رضا نقش نیلچی | شنبه سی ام خرداد ۱۴۰۵ | 16:27

اثر گریز از مرکز در طیف چرخشی (Non-Rigid Rotor)

در مدل روتور صلب فرض کردیم فاصله پیوند ثابت است.

اما در واقعیت، با افزایش سرعت چرخش (افزایش J)، نیروی گریز از مرکز باعث کشیده شدن پیوند می‌شود.

در نتیجه گشتاور لختی کمی افزایش می‌یابد و انرژی واقعی کمی کمتر از مقدار مدل صلب می‌شود.

مرحله ۱: فرمول انرژی در روتور غیرصلب

E(J) = B J(J+1) − D J²(J+1)²

که در آن:

B = rotational constant

D = centrifugal distortion constant

مقدار D بسیار کوچک‌تر از B است.

مرحله ۲: محاسبه نمونه عددی

فرض کنید:

B = 2 cm^-1

D = 0.01 cm^-1

انرژی ترازها:

J=0:

E = 0

J=1:

E = 2×1×2 − 0.01×1²×2²

E = 4 − 0.04 = 3.96 cm^-1

J=2:

E = 2×2×3 − 0.01×4×9

E = 12 − 0.36 = 11.64 cm^-1

J=3:

E = 2×3×4 − 0.01×9×16

E = 24 − 1.44 = 22.56 cm^-1

مرحله ۳: مقایسه با روتور صلب

J	Rigid Rotor	Non-Rigid Rotor
1	4	3.96
2	12	11.64
3	24	22.56

نتیجه مهم

با افزایش J:

Energy decreases slightly compared to rigid model

Spacing between lines becomes smaller

بنابراین در طیف واقعی، فاصله خطوط در J های بالا کمی کاهش می‌یابد.

این انحراف نشانه‌ای از کش آمدن پیوند در اثر چرخش سریع است.

مثال: بررسی ترازهای انرژی چرخشی و الگوی آن‌ها

توسط علی رضا نقش نیلچی | شنبه سی ام خرداد ۱۴۰۵ | 16:24

مثال: بررسی ترازهای انرژی چرخشی و الگوی آن‌ها

در این مثال می‌خواهیم ببینیم ترازهای انرژی (E_J) چگونه چیده می‌شوند و چرا فاصله بین خطوط طیفی (که از تفاضل این ترازها به‌دست می‌آید) ثابت است.

مرحله ۱: فرمول ترازهای انرژی چرخشی

انرژی هر تراز با عدد کوانتومی J برابر است با:

E_J = B J(J+1)

(در اینجا انرژی را بر حسب cm⁻¹ فرض می‌کنیم).

مرحله ۲: محاسبه مقادیر انرژی برای چند تراز اول

فرض می‌کنیم ثابت چرخشی B = 2 cm⁻¹ باشد:

J=0: E₀ = 2 × 0(1) = 0 cm⁻¹

J=1: E₁ = 2 × 1(2) = 4 cm⁻¹

J=2: E₂ = 2 × 2(3) = 12 cm⁻¹

J=3: E₃ = 2 × 3(4) = 24 cm⁻¹

J=4: E₄ = 2 × 4(5) = 40 cm⁻¹

مرحله ۳: محاسبه فاصله بین ترازهای متوالی (خطوط طیفی)

فاصله بین ترازها (ΔE) همان مکان خطوط جذبی است:

ΔE_0→1 = E₁ - E₀ = 4 - 0 = 4 cm⁻¹

ΔE_1→2 = E₂ - E₁ = 12 - 4 = 8 cm⁻¹

ΔE_2→3 = E₃ - E₂ = 24 - 12 = 12 cm⁻¹

ΔE_3→4 = E₄ - E₃ = 40 - 24 = 16 cm⁻¹

مرحله ۴: تحلیل و مشاهده الگو

همانطور که مشاهده می‌کنید، خودِ ترازهای انرژی با افزایش J سریع‌تر بالا می‌روند، اما فاصله بین هر خط طیفی متوالی همواره مقدار ثابتی (یعنی 2B) را به نمایش می‌گذارد.

گذار	انرژی جذب شده (cm⁻¹)
J=0 → 1	4
J=1 → 2	8
J=2 → 3	12
J=3 → 4	16

نکته: این رفتار نردبانی ترازها، ویژگی اصلی روتور صلب است.

مثال: پیش‌بینی موقعیت خطوط چرخشی ایزوتوپ جدید و مقایسه با مولکول اصلی

توسط علی رضا نقش نیلچی | شنبه سی ام خرداد ۱۴۰۵ | 16:22

مثال: پیش‌بینی موقعیت خطوط چرخشی ایزوتوپ جدید و مقایسه با مولکول اصلی

در این مثال می‌خواهیم موقعیت خطوط طیفی دو ایزوتوپ

12CO

و

13CO

را با هم مقایسه کنیم.

فرض می‌کنیم ثابت‌های چرخشی دو مولکول برابر باشند با:

B(12CO) = 1.922 cm^-1

B(13CO) = 1.840 cm^-1

می‌خواهیم موقعیت چند خط اول طیف چرخشی را برای هر دو ایزوتوپ به‌دست آوریم.

مرحله ۱: رابطه مکان خطوط چرخشی

برای گذارهای چرخشی مجاز:

ΔJ = +1

و مکان هر خط از رابطه زیر به دست می‌آید:

ν = 2B(J+1)

اکنون برای چند مقدار اول

J

محاسبه می‌کنیم.

مرحله ۲: خطوط ایزوتوپ 12CO

برای

J = 0

:

ν = 2 × 1.922 × (0+1)

ν = 3.844 cm^-1

برای

J = 1

:

ν = 2 × 1.922 × (1+1)

ν = 7.688 cm^-1

برای

J = 2

:

ν = 2 × 1.922 × (2+1)

ν = 11.532 cm^-1

مرحله ۳: خطوط ایزوتوپ 13CO

برای

J = 0

:

ν = 2 × 1.840 × (0+1)

ν = 3.680 cm^-1

برای

J = 1

:

ν = 2 × 1.840 × (1+1)

ν = 7.360 cm^-1

برای

J = 2

:

ν = 2 × 1.840 × (2+1)

ν = 11.040 cm^-1

مرحله ۴: مقایسه خطوط دو ایزوتوپ

موقعیت خطوط برای

13CO

همیشه کمی کمتر از

12CO

است.

زیرا:

Heavier isotope → larger μ → larger I → smaller B

پس:

Smaller B → smaller line positions

جدول نهایی

J	Line of 12CO	Line of 13CO
0 → 1	3.844 cm^-1	3.680 cm^-1
1 → 2	7.688 cm^-1	7.360 cm^-1
2 → 3	11.532 cm^-1	11.040 cm^-1

نتیجه نهایی

ایزوتوپ سنگین‌تر یعنی

13CO

خطوط طیفی را به عدد موج‌های کوچکتر منتقل می‌کند.

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

مثال: اثر ایزوتوپی در طیف چرخشی مولکول

توسط علی رضا نقش نیلچی | شنبه سی ام خرداد ۱۴۰۵ | 16:21

مثال: اثر ایزوتوپی در طیف چرخشی مولکول CO

دو ایزوتوپ مختلف مولکول CO را در نظر بگیرید:

13CO (با کربن 13)

12CO (با کربن معمولی 12)

اگر ثابت چرخشی ایزوتوپ 12CO برابر باشد با:

B(12CO) = 1.922 cm^-1

ثابت چرخشی ایزوتوپ سنگین‌تر یعنی 13CO را محاسبه کنید.

مرحله ۱: رابطه جرم کاهیده و ثابت چرخشی

ثابت چرخشی وابسته به جرم کاهیده است:

B ∝ 1 / μ

پس نسبت دو ثابت چرخشی:

B(13CO) / B(12CO) = μ(12CO) / μ(13CO)

مرحله ۲: محاسبه جرم‌های کاهیده

برای 12CO:

μ(12CO) = (12 × 16) / (12 + 16) = 6.857 amu

برای 13CO:

μ(13CO) = (13 × 16) / (13 + 16) = 7.172 amu

مرحله ۳: محاسبه نسبت ثابت‌ها

B(13CO) = B(12CO) × μ(12CO) / μ(13CO)

B(13CO) = 1.922 × (6.857 / 7.172)

B(13CO) = 1.922 × 0.9561

B(13CO) ≈ 1.84 cm^-1

نتیجه نهایی

B(13CO) ≈ 1.84 cm^-1

ایزوتوپ سنگین‌تر → جرم کاهیده بزرگ‌تر → ثابت چرخشی کوچکتر → خطوط طیفی به عدد موج کمتر جابه‌جا می‌شوند.

h3>مثال: محاسبه فاصله تعادلی پیوند (r) از گشتاور لختی

توسط علی رضا نقش نیلچی | شنبه سی ام خرداد ۱۴۰۵ | 16:18

در این مرحله، از گشتاور لختی که در مرحله قبل به دست آوردیم استفاده می‌کنیم تا **فاصله پیوند (Bond Length)** را محاسبه کنیم. فرض می‌کنیم مولکول مورد نظر **Carbon Monoxide (CO)** باشد.

:::writing

مثال: محاسبه فاصله تعادلی پیوند (r) از گشتاور لختی

با فرض اینکه مولکول مورد نظر CO باشد و گشتاور لختی آن در مرحله قبل محاسبه شده باشد:

I = 1.45 × 10^-46 kg·m²

مرحله ۱: محاسبه جرم کاهیده (Reduced Mass)

برای مولکول CO با جرم‌های اتمی C=12 و O=16:

μ = (m₁ × m₂) / (m₁ + m₂)

μ = (12 × 16) / (12 + 16) = 192 / 28 = 6.857 amu

تبدیل واحد amu به کیلوگرم:

1 amu = 1.66 × 10^-27 kg

μ = 6.857 × 1.66 × 10^-27 = 1.138 × 10^-26 kg

مرحله ۲: رابطه گشتاور لختی و فاصله پیوند

I = μ r²

بنابراین برای یافتن فاصله پیوند (r):

r = √(I / μ)

مرحله ۳: جایگذاری و محاسبه نهایی

r = √(1.45 × 10^-46 / 1.138 × 10^-26)

r = √(1.274 × 10^-20)

r ≈ 1.128 × 10^-10 m

تبدیل به واحد آنگستروم:

r ≈ 1.128 Å

نتیجه نهایی

Bond Length (r) = 1.128 Å

:::

در مرحله بعد، اگر تمایل داشته باشید، می‌توانیم به سراغ **اثر ایزوتوپی (Isotope Effect)** برویم که نشان می‌دهد چگونه تغییر یک اتم (مثلاً از C12 به C13) باعث جابجایی خطوط طیفی می‌شود. مایل هستید؟

مثال: محاسبه گشتاور لختی مولکول از ثابت چرخشی

توسط علی رضا نقش نیلچی | شنبه سی ام خرداد ۱۴۰۵ | 16:16

مثال: محاسبه گشتاور لختی مولکول از ثابت چرخشی B

در مثال قبل به‌دست آوردیم که ثابت چرخشی مولکول برابر است با:

B = 10 cm^-1

حال می‌خواهیم گشتاور لختی مولکول را به‌دست آوریم.

مرحله ۱: رابطه بین B و I

B = h / (8π² c I)

بنابراین:

I = h / (8π² c B)

مرحله ۲: تبدیل B به واحد SI

باید B را از cm^-1 به m^-1 تبدیل کنیم:

10 cm^-1 = 1000 m^-1

مرحله ۳: قرار دادن مقادیر در فرمول

h = 6.626 × 10^-34 J·s

c = 3.00 × 10⁸ m/s

I = 6.626×10^-34 / (8π² × 3.00×10⁸ × 1000)

I ≈ 2.8 × 10^-47 kg·m²

نتیجه نهایی

I = 2.8 × 10^-47 kg·m²

در مرحله بعد (اگر بخواهید) فاصله پیوند مولکول r را از رابطه زیر محاسبه می‌کنیم:

I = μ r²

مثال: محاسبه ثابت چرخشی از خطوط متوالی طیف چرخشی

توسط علی رضا نقش نیلچی | شنبه سی ام خرداد ۱۴۰۵ | 16:14

مثال: محاسبه ثابت چرخشی از خطوط متوالی طیف چرخشی

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

20 cm^-1 , 40 cm^-1 , 60 cm^-1 , 80 cm^-1

ثابت چرخشی مولکول را محاسبه کنید.

مرحله ۱: رابطه ترازهای چرخشی

E(J) = B J(J+1)

اختلاف انرژی بین دو تراز متوالی:

ΔE = E(J+1) − E(J)

ΔE = B[(J+1)(J+2) − J(J+1)]

ΔE = 2B(J+1)

مرحله ۲: موقعیت خطوط طیفی

ν = 2B(J+1)

بنابراین فاصله بین خطوط متوالی برابر است با:

Δν = 2B

مرحله ۳: محاسبه فاصله خطوط

40 − 20 = 20

60 − 40 = 20

80 − 60 = 20

پس فاصله خطوط طیفی:

Δν = 20 cm^-1

مرحله ۴: محاسبه ثابت چرخشی

2B = 20

B = 10 cm^-1

نتیجه

Rotational constant:

B = 10 cm^-1

:::

اگر موافق باشید، مثال بعدی که به ترتیب معمول در درس طیف‌سنجی می‌آید این است:

**محاسبه گشتاور لختی از ثابت چرخش(B** (مرحله بعد از همین مثال).

مثال: محاسبه ثابت چرخشی و فاصله پیوند مولکول

توسط علی رضا نقش نیلچی | شنبه سی ام خرداد ۱۴۰۵ | 16:5

HCl

خط جذبی گذار چرخشی مولکول HCl از J=2 → J=3 در عدد موج

300 cm^-1 مشاهده شده است.

ثابت چرخشی B و فاصله تعادلی پیوند را محاسبه کنید.

مرحله ۱: رابطه گذار چرخشی

For a rigid rotor:

E(J) = B J(J+1)

Transition energy:

ΔE = E(J+1) − E(J) = 2B(J+1)

برای J=2:

ΔE = 2B(3) = 6B

از آنجا که عدد موج داده شده برابر 300 cm^-1 است:

6B = 300

B = 50 cm^-1

مرحله ۲: محاسبه گشتاور لختی

B = h / (8π²cI)

I = h / (8π²cB)

تبدیل B به واحد SI:

B = 50 cm^-1 = 5000 m^-1

h = 6.626 × 10^-34 J.s

c = 3.00 × 10⁸ m/s

I = 6.626×10^-34 / (8π² × 3.00×10⁸ × 5000)

I ≈ 5.6 × 10^-47 kg.m²

مرحله ۳: محاسبه فاصله تعادلی

I = μ r²

جرم کاهیده:

μ = (m_H m_Cl) / (m_H + m_Cl)

m_H = 1u

m_Cl = 35.5u

μ = (1 × 35.5) / (36.5) u

μ ≈ 0.973 u

1u = 1.6605 × 10^-27 kg

μ ≈ 1.61 × 10^-27 kg

r = √(I / μ)

r = √(5.6×10^-47 / 1.61×10^-27)

r ≈ 1.86 × 10^-10 m

r ≈ 1.86 Å

نتیجه نهایی

B = 50 cm^-1

I ≈ 5.6 × 10^-47 kg.m²

r ≈ 1.86 Å

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

توسط علی رضا نقش نیلچی | شنبه سی ام خرداد ۱۴۰۵ | 15:57

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

در این مثال، می‌خواهیم انرژی حالت پایه یک نوسانگر هماهنگ یک‌بعدی را با استفاده از یک تابع آزمون نمایی (گاوسی) تقریب بزنیم.

ورودی‌های مسئله:

Hamiltonian: H = -(ħ² / 2m) (d²/dx²) + (1/2) m ω² x²

Trial Function: ψ(x) = exp(-α x²)

گام 1: محاسبه نرمال‌سازی (Normalization)

ابتدا انتگرال مربع تابع موج را محاسبه می‌کنیم تا مخرج کسر انرژی مشخص شود:

<ψ|ψ> = ∫_{-∞}^{+∞} exp(-2α x²) dx = sqrt(π / 2α)

گام 2: محاسبه مقدار چشم‌داشتی انرژی پتانسیل

پتانسیل ما V = (1/2)mω²x² است. انتگرال آن روی تابع آزمون برابر است با:

<V> = [ ∫_{-∞}^{+∞} exp(-α x²) (1/2 m ω² x²) exp(-α x²) dx ] / <ψ|ψ>

<V> = m ω² / (8α)

گام 3: محاسبه مقدار چشم‌داشتی انرژی جنبشی

با مشتق‌گیری از تابع آزمون و جایگذاری در فرمول انرژی جنبشی:

<T> = [ ∫_{-∞}^{+∞} ψ* (-(ħ²/2m) d²/dx²) ψ dx ] / <ψ|ψ>

<T> = ħ² α / (2m)

گام 4: تشکیل تابع انرژی کل E(α)

انرژی کل تابعی از پارامتر α خواهد بود:

E(α) = <T> + <V> = (ħ² α / 2m) + (m ω² / 8α)

گام 5: کمینه‌سازی (یافتن α بهینه)

از تابع انرژی نسبت به α مشتق می‌گیریم و برابر صفر قرار می‌دهیم:

dE/dα = (ħ² / 2m) - (m ω² / 8α²) = 0

8α² ħ² = 2m² ω²

α² = m² ω² / (4ħ²)

α_opt = m ω / (2ħ)

گام 6: محاسبه انرژی نهایی (E_min)

حالا مقدار α_opt را در تابع E(α) جایگذاری می‌کنیم:

E_min = (ħ² / 2m) * (m ω / 2ħ) + (m ω² / 8) * (2ħ / m ω)

E_min = (1/4) ħω + (1/4) ħω

E_min = (1/2) ħω

تحلیل نتیجه ردگیری

در این مثال خاص، چون تابع آزمون ما (گاوسی) دقیقاً هم‌شکل با تابع موج واقعی نوسانگر هماهنگ بود، انرژی به‌دست آمده (1/2 ħω) دقیقاً با جواب واقعی برابر شد.

Step	Input	Operation	Output
1	ψ = exp(-αx²)	∫\|ψ\|² dx	N = sqrt(π/2α)
2	V = 1/2 mω²x²	<V> Calculation	mω² / 8α
3	T = -ħ²/2m ∇²	<T> Calculation	ħ²α / 2m
4	T and V	Summation	E(α)
5	E(α)	dE/dα = 0	α_opt = mω/2ħ
6	α_opt	Substitution	E = 1/2 ħω

آموزش ریاضی

طیف چرخشی رامان و فاصله خطوط طیفی (4B)

طیف چرخشی رامان و فاصله خطوط 4B

۱. قانون گزینش در رامان چرخشی

۲. انرژی ترازهای چرخشی

۳. محاسبه مکان خطوط رامان

۴. موقعیت چند خط اول

۵. فاصله خطوط طیفی

۶. نتیجه مهم

طیف‌سنجی رامان (Raman Spectroscopy)

طیف‌سنجی رامان (Raman Spectroscopy)

۱. انواع پراکندگی نور

۲. شرط فعال بودن در رامان

۳. مثال مهم

۴. تفاوت خطوط Stokes و Anti‑Stokes

۵. شیفت رامان

بررسی حل تحلیلی معادله شرودینگر برای پتانسیل خطی

بررسی حل تحلیلی معادله شرودینگر برای پتانسیل خطی

کاربرد نظریه اختلال در اتم هلیوم (حالت سینگلت)

برهم‌کنش بین دو الکترون

مسئله برهم‌کنش بین دو الکترون در اتم هلیوم

حل قسمت (ب)

1) تابع موج حالت پایه

2) حالت اسپینی سینگلت

3) تصحیح انرژی مرتبه اول

4) محاسبه انتگرال

نتیجه

اتم هلیوم با برهم‌کنش خطی

مثال چهارم: روش وردشی برای نوسانگر آنهارمونیک (پتانسیل x⁴)

مثال چهارم: روش وردشی برای نوسانگر آنهارمونیک (پتانسیل x⁴)

1) نرمال‌سازی

2) انرژی جنبشی

3) انرژی پتانسیل هارمونیک

4) جمله آنهارمونیک

5) انرژی کل

6) کمینه کردن انرژی

نتیجه

مثال سوم: روش وردشی برای اتم هیدروژن

مثال سوم: روش وردشی برای اتم هیدروژن

1) نرمال‌سازی

2) انرژی جنبشی

3) انرژی پتانسیل

4) انرژی کل

5) کمینه کردن انرژی

6) انرژی حالت پایه

مثال دوم: روش وردشی برای چاه پتانسیل بی‌نهایت یک‌بعدی

مثال دوم: روش وردشی برای چاه پتانسیل بی‌نهایت یک‌بعدی

1) نرمال‌سازی

2) انرژی جنبشی

3) انرژی نهایی

مثال: الگوریتم محاسبه انرژی ذره در جعبه یک‌بعدی

مثال: الگوریتم محاسبه انرژی ذره در جعبه یک‌بعدی

روش وردشی برای نوسانگر هماهنگ سه‌بعدی

حل بخش (الف): روش وردشی برای نوسانگر هماهنگ سه‌بعدی

1) محاسبه نرمال‌سازی

2) انرژی پتانسیل

3) انرژی جنبشی

4) انرژی کل

5) کمینه کردن انرژی

6) انرژی مینیمم

الگوریتم محاسبه مسئله ذره در جعبه سه‌بعدی

الگوریتم محاسبه مسئله ذره در جعبه سه‌بعدی

حل مسئله ذره در جعبه سه‌بعدی در حضور میدان الکتریکی و مغناطیسی

حل مسئله ذره در جعبه سه‌بعدی در حضور میدان الکتریکی و مغناطیسی

الف) معادله شرودینگر ذره

ب) بررسی امکان انتخاب میدان مغناطیسی مناسب

حالت پایه

اولین حالت برانگیخته

اثر میدان الکتریکی

تصحیح انرژی در مرتبه اول

ج) دامنه گذار به اولین حالت برانگیخته در تقریب دوقطبی الکتریکی

جمع‌بندی نهایی

برهم‌کنش ارتعاش-چرخش و تغییر فاصله خطوط

برهم‌کنش ارتعاش-چرخش و تغییر فاصله خطوط

۱. وابستگی B به تراز ارتعاشی

۲. اثر بر شاخه‌های P و R

۳. فرمول دقیق‌تر مکان خطوط

مثال مفهومی

نتیجه‌گیری

مرحله ۲: یافتن تراز با بیشترین جمعیت (J_max)