קירוב WKB

בפיזיקה מתמטית, קירוב WKB או **שיטת WKB** היא שיטת קירוב ליניארי לצורך מציאת פתרונות של משוואות דיפרנציאליות חלקיות עם מקדמים המשתנים במרחב. בדרך כלל, נעשה שימוש בשיטה זו עבור חישוב סמי-קלאסי במכניקה קוונטית, שבו פונקציית הגל מתנהגת כפונקציה אקספוננציאלית, המתפשטת סמי-קלאסית במרחב ואשר או המשרעת או הפאזה של פונקציה זו משתנים באופן איטי.

מקור השם WKB הוא בראשי תיבות של המדענים ונצל-קרמרס-ברילואן (**W**entzel–**K**ramers–**B**rillouin). קירוב זה ידוע גם כקירוב LG או קירוב ליוביל-גרין (**L**iouville–**G**reen). לעיתים הוא נקרא  JWKB או WKBJ כאשר ה-J מתייחסת לג'פרי (**J**effreys).

## היסטוריה
קירוב WKB קרוי על שם הפיזיקאים: גרגור ונצל, הנס קרמרס ולאון ברילואן (**W**entzel–**K**ramers–**B**rillouin), אשר פיתחו אותה בשנת 1926.

בשנת 1923, פיתח המתמטיקאי הבריטי הרולד ג'פרי שיטה כללית של קירוב לפתרונות משוואות ליניאריות דיפרנציאליות מסדר שני. שנתיים מאוחר יותר פותחה משוואת שרדינגר וונצל, קרמרס, וברילואן שלא היו מודעים לעבודה המוקדמת של ג'פרי, פיתחו קירוב ליניארי כללי שניתן ליישמו כפתרון למשוואת שרדינגר.

מקורות מוקדמים יותר לשיטה ניתן למצוא בעבודתם של : קרליני (1817), ליוביל (1837), גרין 1837, ריילי (1912) וגנז בשנת 1915. יש הסבורים כי ניתן לזהות את ז'וזף ליוביל וג'ורג' גרין כאבות השיטה ולכן קירוב זה ידוע גם כקירוב LG או קירוב ליוביל-גרין (**L**iouville–**G**reen).

תרומתם של ונצל, קרמרס, וברילואן לשיטה הייתה הכללת הטיפול של נקודות מפנה , על ידי "תפירת" הפתרונות התנודתיים והזמניים  משני צדי נקודת המפנה. לדוגמה, פתרון משוואת שרדינגר עם גבעת אנרגיה פוטנציאלית (אזור סביב מקסימום מקומי של אנרגיה פוטנציאלית).

## שיטת WKB
באופן כללי, תאוריית WKB היא שיטה למציאת קירוב לפתרון משוואות דיפרנציאליות אשר נגזרתה הגבוהה ביותר מוכפלת בפרמטר קטן, \epsilon:

\epsilon \frac{\mathrm{d}^ny}{\mathrm{d}x^n} + a(x)\frac{\mathrm{d}^{n-1}y}{\mathrm{d}x^{n-1}} + \cdots + k(x)\frac{\mathrm{d}y}{\mathrm{d}x} + m(x)y= 0 ~\quad(1)

מניחים פתרון בצורת הרחבה לטור אסימפטוטי:

 y(x) \sim \exp\left[\frac{1}{\delta}\sum_{n=0}^{\infty}\delta^nS_n(x)\right] \quad(2)

בגבול שבו \delta\rightarrow0.

## דוגמה

מניחים משוואה דיפרנציאלית ליניארית מסדר שני כגון:

\epsilon^2 \frac{d^2 y}{dx^2} = Q(x) y\quad(3)

כאשר Q(x) \neq 0.

מציבים:

y(x) = \exp\left[\frac{1}{\delta}\sum_{n=0}^\infty \delta^nS_n(x)\right] \quad(4)

ומקבלים:

\epsilon^2\left[\frac{1}{\delta^2}\left(\sum_{n=0}^\infty \delta^nS_n'\right)^2 + \frac{1}{\delta}\sum_{n=0}^{\infty}\delta^nS_n*\right] = Q(x)\quad(5)</math>
</center>
עבור הסדר המוביל, אפשר לקרב את משוואה <math>\quad(5)</math> :

<center><math>\frac{\epsilon^2}{\delta^2}S_0'^2 + \frac{2\epsilon^2}{\delta}S_0'S_1' + \frac{\epsilon^2}{\delta}S_0* = Q(x)\quad(6)

בגבול שבו \delta\rightarrow0, האיבר הדומיננטי ניתן על ידי:

\frac{\epsilon^2}{\delta^2}S_0'^2 \sim Q(x) \quad(7)

כך ש-\delta פרופורציונלי ל-\epsilon.

הצבת 2 גדלים אלו כשווים והשוואת החזקות של \epsilon ב"סדר 0" מניבה:

\epsilon^0: \quad S_0'^2 = Q(x) \quad(8)

המוכרת כמשוואה איקונלית  עם הפתרון:

S_0(x) = \pm \int_{x_0}^x \sqrt{Q(t)}\,dt \quad(9)השוואת החזקות של \epsilon בסדר הראשון קובעת:

\epsilon^1: \quad 2S_0'S_1' + S_0* = 0 \quad(10)</math>

המוכרת כ[[משוואת מעבר]]  עם הפתרון:

<math display="block">S_1(x) = -\frac{1}{4}\ln Q(x) + k_1 \quad(11)</math>

כאשר <math>k_1 </math> הוא קבוע שרירותי.

מאחר ש-<math>S_0</math> יכול להיות הן חיובי והן שלילי, קירוב WKB של <math>\quad(3)</math> בסדר ראשון יהיה [[צירוף ליניארי]] של זוג פתרונות כדלהלן:

<center><math>y(x) \approx c_1Q^{-\frac{1}{4}}(x)\exp\left[\frac{1}{\epsilon}\int_{x_0}^x\sqrt{Q(t)}dt\right] + c_2Q^{-\frac{1}{4}}(x)\exp\left[-\frac{1}{\epsilon}\int_{x_0}^x\sqrt{Q(t)}dt\right]\quad(12)</math>
</center>

ניתן לקבל פתרונות לסדרים גבוהים יותר <math>(n \ge 2)</math> על ידי חישוב דומה של:

<center><math> 2S_0'S_n' + S*_{n-1} + \sum_{j=1}^{n-1}S'_jS'_{n-j} = 0\quad(13)


## יישום עבור משוואת שרדינגר
משוואת שרדינגר שאינה תלויה בזמן במימד אחד היא:

-\frac{\hbar^2}{2m} \frac{\mathrm{d}^2}{\mathrm{d}x^2} \Psi(x) + V(x) \Psi(x) = E \Psi(x) \quad(14)

ניתן לכתוב משוואה זו בצורה אחרת:

 \frac{\mathrm{d}^2}{\mathrm{d}x^2} \Psi(x) = \frac{2m}{\hbar^2} \left( V(x) - E \right) \Psi(x) \quad(15)


פונקציית הגל ניתנת להגדרה על ידי אקספוננט של פונקציה אחרת, *<math>\Phi </math>* (המשויכת לפעולה), אשר יכולה להיות גם פונקציה מרוכבת:

\Psi(x) = e^{\Phi(x)} ~\quad(16)

כך ש:

\Phi*(x) + \left[\Phi'(x)\right]^2 = \frac{2m}{\hbar^2} \left( V(x) - E \right)\quad(17)</math> 
</center>

כאשר *\Phi' * היא הנגזרת של *\Phi * לפי *x *.

הנגזרת *\Phi' *, יכולה להתפצל לחלק ממשי ולחלק מדומה הכוללים את הפונקציות הממשיות *A * ו-*B *:

<center><math>\Phi'(x) = A(x) + i B(x)~\quad(18)</math> 
</center>

ואז [[משרעת]] [[פונקציית גל|פונקציית הגל]] היא <math>\exp\left[\int^x A(x')dx'\right]\,\!~</math> ואילו ה[[פאזה (גלים)|פאזה]] היא: <math>\int^x B(x')dx'\,\!~</math>.

[[משוואת שרדינגר]] מתחלקת ל-2 משוואות: אחת עבור החלק ה[[מספר ממשי|ממשי]] <math>\quad(19)</math> ואחת עבור החלק ה[[מספר מדומה|דמיוני]] <math>\quad(20)</math> כדלהלן:

<center><math>A'(x) + A(x)^2 - B(x)^2 = \frac{2m}{\hbar^2} \left( V(x) - E \right) ~\quad(19)</math> 
</center>

<math display="block">B'(x) + 2 A(x) B(x) = 0 ~\quad(20)</math>

לאחר מכן, משתמשים בקירוב ה[[פיזיקה סמי-קלאסית|סמי-קלאסי]], קירוב WKB. משמעות הדבר היא כי כל פונקציה מפותחת כ[[טור חזקות]] של [[קבוע פלאנק המצומצם]], <math>\hbar</math>. מתוך המשוואות לעיל, ניתן לראות כי טור החזקות חייב להתחיל עם סדר של <math>\frac{1}{ \hbar}</math> כדי להתאים לחלק הממשי של המשוואה.

על מנת להשיג גבול קלאסי טוב, יש צורך להתחיל עם חזקה גבוהה של [[קבוע פלאנק המצומצם]], <math>\hbar</math> ככל האפשר:

<center><math>A(x) = \frac{1}{\hbar} \sum_{n=0}^\infty \hbar^n A_n(x)~\quad(21)</math>

</center><math display="block">B(x) = \frac{1}{\hbar} \sum_{n=0}^\infty \hbar^n B_n(x)~\quad(22)</math>

עבור "סדר 0" ניתן להביע את *A * ו-*B * כך:

<center><math>A_0(x)^2 - B_0(x)^2 = 2m \left( V(x) - E \right)\quad(23)</math>
</center>

<math display="block">A_0(x) B_0(x) = 0 \;\quad(24)</math>

הנגזרות הראשונות *A'(x) * ו-*B'(x) '' נמחקו מאחר שהן כוללות פקטורים מסדר גודל של \frac{1}{ \hbar} הגבוה יחסית לסדר הדומיננטי \frac{1}{\hbar^2}.

מצד אחד, אם המשרעת משתנה לאט יחסית בהשוואה לפאזה (A_0(x) = 0), אזי מקבלים:

B_0(x) = \pm \sqrt{ 2m \left( E - V(x) \right) }\quad(24)

המשוואה שלעיל, \quad(24), תקפה רק כאשר האנרגיה הכללית גדולה יותר מהאנרגיה הפוטנציאלית (E>V(x)), כפי שקורה תמיד בתנועה קלאסית.

אחרי ביצוע תהליך דומה לסדר הבא בפיתוח לטור חזקות, מקבלים כפתרון למשוואת שרדינגר את פונקציית הגל:

\Psi(x) \approx C_0 \frac{ e^{i \int \mathrm{d}x \sqrt{\frac{2m}{\hbar^2} \left( E - V(x) \right)} + \theta} }{\sqrt[4]{\frac{2m}{\hbar^2} \left( E - V(x) \right)}}~\quad(25)

מצד שני, אם הפאזה משתנה לאט יחסית בהשוואה למשרעת (B_0(x) = 0), אזי מקבלים:

A_0(x) = \pm \sqrt{ 2m \left( V(x) - E \right) }\quad(26)

המשוואה שלעיל, \quad(26), תקפה רק כאשר האנרגיה הפוטנציאלית גדולה יותר מהאנרגיה הכללית (E, כפי שקורה תמיד במנהור קוונטי.

אחרי ביצוע תהליך דומה לסדר הבא בפיתוח לטור חזקות, כמקודם, מקבלים כפתרון למשוואת שרדינגר את פונקציית הגל:

 \Psi(x) \approx \frac{ C_{+} e^{+\int \mathrm{d}x \sqrt{\frac{2m}{\hbar^2} \left( V(x) - E \right)}} 
+ C_{-} e^{- \int \mathrm{d}x \sqrt{\frac{2m}{\hbar^2} \left( V(x) - E \right)}}}{\sqrt[4]{\frac{2m}{\hbar^2} \left( V(x) - E \right)}}~\quad(27)

2 הפתרונות המקורבים למשוואת שרדינגר, הן \quad(25) והן \quad(27), נעשים סינגולריים בנקודת המפנה או נקודת ה"תפר", כאשר E=V(x) ולכן אינם תקפים בנקודה זו.

תקפות הקירובים היא מעל ומתחת גבעת הפוטנציאל. מעל גבעת הפוטנציאל החלקיק מתנהג כגל חופשי - פונקציית הגל אוסצילטורית. מתחת גבעת הפוטנציאל החלקיק עובר שינויים אקספוננטציאלים במשרעת.

על מנת להבטיח את הגלובליות של קירוב WKB עבור משוואת שרדינגר הכוללת גבעת פוטנציאל נדרש למצוא פתרון מקורב בנקודות המפנה שבהן E=V(x). הדבר נעשה בדרך הבאה:

עבור נקודת המפנה הקלאסית, x_1 וקרוב ל- E=V(x_1), הביטוי \frac{2m}{\hbar^2}\left(V(x)-E\right) ניתן לפיתוח כטור חזקות:

\frac{2m}{\hbar^2}\left(V(x)-E\right) = U_1 \cdot (x - x_1) + U_2 \cdot (x - x_1)^2 + \cdots\;\quad(28)


עבור הסדר הראשון מתקבל:

\frac{\mathrm{d}^2}{\mathrm{d}x^2} \Psi(x) = U_1 \cdot (x - x_1) \cdot \Psi(x)\quad(29)


משוואה דיפרנציאלית זו ידועה כמשוואת איירי וניתן להביע את פתרונה במונחים של פונקציית איירי:

\Psi(x) = C_A \textrm{Ai}\left( \sqrt[3]{U_1} \cdot (x - x_1) \right) + C_B \textrm{Bi}\left( \sqrt[3]{U_1} \cdot (x - x_1) \right)\quad(30)

"תפירת" פתרון זה עם 2 הפתרונות המקורבים מעל ומתחת גבעת הפוטנציאל נעשית על ידי מציאת היחס בין המקדמים השונים, C_0,\theta ו-C_{+},C_{-}. מאחר שפונקציות איירי בגבול האסימפטוטי, הופכות לפונקציות אקספוננטציאליות כגון: סינוס וקוסינוס, ניתן להביע את היחס בין המקדמים השונים דרך "פונקציות הקשר":

 C_{+} = + \frac{1}{2} C_0 \cos{\left(\theta - \frac{\pi}{4}\right)}\quad(31)


 C_{-} = - \frac{1}{2} C_0 \sin{\left(\theta - \frac{\pi}{4}\right)}\quad(32)

לסיכום, מתקבל פתרון מקורב גלובלי עבור משוואת שרדינגר הכוללת גבעת פוטנציאל.

## ראו גם
* תורת ההפרעות
* מקדם העברה
* מנהור קוונטי

## לקריאה נוספת


* Bender, Carl; Orszag, Steven (1978). Advanced Mathematical Methods for Scientists and Engineers. McGraw-Hill. 
* Child, M. S. (1991). Semiclassical mechanics with molecular applications. Oxford: Clarendon Press. .
* Griffiths, David J. (2004). Introduction to Quantum Mechanics (2nd ed.). Prentice Hall. .
* Hall, B.C. (2013). Quantum Theory for Mathematicians. Springer.
* Liboff, Richard L. (2003). Introductory Quantum Mechanics (4th ed.). Addison-Wesley. .
* Olver, Frank William John (1974). Asymptotics and Special Functions. Academic Press. .
* Razavy, Mohsen (2003). Quantum Theory of Tunneling. World Scientific. .
* Sakurai, J. J. (1993). Modern Quantum Mechanics. Addison-Wesley. .


## הערות שוליים

קטגוריה:פיזיקה תאורטית
קטגוריה:מכניקת הקוונטים