בפיזיקה מתמטית, קירוב 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 :



ϵ


d

n

y

d

x

n


+
a
(
x
)


d

n
−
1

y

d

x

n
−
1


+
⋯
+
k
(
x
)

d

y

d

x

+
m
(
x
)
y
=
0
 

(
1
)

{\displaystyle \epsilon {\frac {\mathrm {d} ^{n}y}{\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)}

{\displaystyle \epsilon {\frac {\mathrm {d} ^{n}y}{\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
)
∼
exp
⁡

[

1
δ

∑

n
=
0

∞

δ

n

S

n

(
x
)

]

(
2
)

{\displaystyle y(x)\sim \exp \left[{\frac {1}{\delta }}\sum \_{n=0}^{\infty }\delta ^{n}S\_{n}(x)\right]\quad (2)}

{\displaystyle y(x)\sim \exp \left[{\frac {1}{\delta }}\sum _{n=0}^{\infty }\delta ^{n}S_{n}(x)\right]\quad (2)}
בגבול שבו {\displaystyle \delta \rightarrow 0}.

דוגמה
-----

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


ϵ

2



d

2

y

d

x

2


=
Q
(
x
)
y

(
3
)

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

{\displaystyle \epsilon ^{2}{\frac {d^{2}y}{dx^{2}}}=Q(x)y\quad (3)}
כאשר {\displaystyle Q(x)\neq 0}.

מציבים:



y
(
x
)
=
exp
⁡

[

1
δ

∑

n
=
0

∞

δ

n

S

n

(
x
)

]

(
4
)

{\displaystyle y(x)=\exp \left[{\frac {1}{\delta }}\sum \_{n=0}^{\infty }\delta ^{n}S\_{n}(x)\right]\quad (4)}

{\displaystyle y(x)=\exp \left[{\frac {1}{\delta }}\sum _{n=0}^{\infty }\delta ^{n}S_{n}(x)\right]\quad (4)}
ומקבלים:

{\displaystyle \epsilon ^{2}\left[{\frac {1}{\delta ^{2}}}\left(\sum \_{n=0}^{\infty }\delta ^{n}S\_{n}'\right)^{2}+{\frac {1}{\delta }}\sum \_{n=0}^{\infty }\delta ^{n}S\_{n}''\right]=Q(x)\quad (5)}

עבור הסדר המוביל, אפשר לקרב את משוואה {\displaystyle \quad (5)} :

{\displaystyle {\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)}

בגבול שבו {\displaystyle \delta \rightarrow 0}, האיבר הדומיננטי ניתן על ידי:



ϵ

2

δ

2


S

0

′

2

∼
Q
(
x
)

(
7
)

{\displaystyle {\frac {\epsilon ^{2}}{\delta ^{2}}}S\_{0}'^{2}\sim Q(x)\quad (7)}

{\displaystyle {\frac {\epsilon ^{2}}{\delta ^{2}}}S_{0}'^{2}\sim Q(x)\quad (7)}
כך ש-\delta פרופורציונלי ל-\epsilon .

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


ϵ

0

:

S

0

′

2

=
Q
(
x
)

(
8
)

{\displaystyle \epsilon ^{0}:\quad S\_{0}'^{2}=Q(x)\quad (8)}

{\displaystyle \epsilon ^{0}:\quad S_{0}'^{2}=Q(x)\quad (8)}
המוכרת כמשוואה איקונלית (אנ') עם הפתרון:

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


ϵ

1

:

2

S

0

′

S

1

′

+

S

0

″

=
0

(
10
)

{\displaystyle \epsilon ^{1}:\quad 2S\_{0}'S\_{1}'+S\_{0}''=0\quad (10)}

{\displaystyle \epsilon ^{1}:\quad 2S_{0}'S_{1}'+S_{0}''=0\quad (10)}
המוכרת כמשוואת מעבר (אנ') עם הפתרון:


S

1

(
x
)
=
−

1
4

ln
⁡
Q
(
x
)
+

k

1

(
11
)

{\displaystyle S\_{1}(x)=-{\frac {1}{4}}\ln Q(x)+k\_{1}\quad (11)}

{\displaystyle S_{1}(x)=-{\frac {1}{4}}\ln Q(x)+k_{1}\quad (11)}
כאשר {\displaystyle k\_{1}} הוא קבוע שרירותי.

מאחר ש-{\displaystyle S\_{0}} יכול להיות הן חיובי והן שלילי, קירוב WKB של {\displaystyle \quad (3)} בסדר ראשון יהיה צירוף ליניארי של זוג פתרונות כדלהלן:

{\displaystyle y(x)\approx c\_{1}Q^{-{\frac {1}{4}}}(x)\exp \left[{\frac {1}{\epsilon }}\int \_{x\_{0}}^{x}{\sqrt {Q(t)}}dt\right]+c\_{2}Q^{-{\frac {1}{4}}}(x)\exp \left[-{\frac {1}{\epsilon }}\int \_{x\_{0}}^{x}{\sqrt {Q(t)}}dt\right]\quad (12)}

ניתן לקבל פתרונות לסדרים גבוהים יותר {\displaystyle (n\geq 2)} על ידי חישוב דומה של:

{\displaystyle 2S\_{0}'S\_{n}'+S''\_{n-1}+\sum \_{j=1}^{n-1}S'\_{j}S'\_{n-j}=0\quad (13)}

יישום עבור משוואת שרדינגר
-------------------------

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

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

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

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

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



Ψ
(
x
)
=

e

Φ
(
x
)

 

(
16
)

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

{\displaystyle \Psi (x)=e^{\Phi (x)}~\quad (16)}
כך ש:

{\displaystyle \Phi ''(x)+\left[\Phi '(x)\right]^{2}={\frac {2m}{\hbar ^{2}}}\left(V(x)-E\right)\quad (17)}

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

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

{\displaystyle \Phi '(x)=A(x)+iB(x)~\quad (18)}

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

משוואת שרדינגר מתחלקת ל-2 משוואות: אחת עבור החלק הממשי {\displaystyle \quad (19)} ואחת עבור החלק הדמיוני {\displaystyle \quad (20)} כדלהלן:

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



B
′

(
x
)
+
2
A
(
x
)
B
(
x
)
=
0
 

(
20
)

{\displaystyle B'(x)+2A(x)B(x)=0~\quad (20)}

{\displaystyle B'(x)+2A(x)B(x)=0~\quad (20)}
לאחר מכן, משתמשים בקירוב הסמי-קלאסי, קירוב WKB. משמעות הדבר היא כי כל פונקציה מפותחת כטור חזקות של קבוע פלאנק המצומצם, \hbar. מתוך המשוואות לעיל, ניתן לראות כי טור החזקות חייב להתחיל עם סדר של {\displaystyle {\frac {1}{\hbar }}} כדי להתאים לחלק הממשי של המשוואה.

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

{\displaystyle A(x)={\frac {1}{\hbar }}\sum \_{n=0}^{\infty }\hbar ^{n}A\_{n}(x)~\quad (21)}

B
(
x
)
=

1
ℏ

∑

n
=
0

∞

ℏ

n

B

n

(
x
)
 

(
22
)

{\displaystyle B(x)={\frac {1}{\hbar }}\sum \_{n=0}^{\infty }\hbar ^{n}B\_{n}(x)~\quad (22)}

{\displaystyle B(x)={\frac {1}{\hbar }}\sum _{n=0}^{\infty }\hbar ^{n}B_{n}(x)~\quad (22)}
עבור "סדר 0" ניתן להביע את *{\displaystyle A}* ו-*{\displaystyle B}* כך:

{\displaystyle A\_{0}(x)^{2}-B\_{0}(x)^{2}=2m\left(V(x)-E\right)\quad (23)}



A

0

(
x
)

B

0

(
x
)
=
0

(
24
)

{\displaystyle A\_{0}(x)B\_{0}(x)=0\;\quad (24)}

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

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


B

0

(
x
)
=
±

2
m

(

E
−
V
(
x
)

)

(
24
)

{\displaystyle B\_{0}(x)=\pm {\sqrt {2m\left(E-V(x)\right)}}\quad (24)}

{\displaystyle B_{0}(x)=\pm {\sqrt {2m\left(E-V(x)\right)}}\quad (24)}
המשוואה שלעיל, {\displaystyle \quad (24)}, תקפה רק כאשר האנרגיה הכללית גדולה יותר מהאנרגיה הפוטנציאלית {\displaystyle (E>V(x))}, כפי שקורה תמיד בתנועה קלאסית.

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



Ψ
(
x
)
≈

C

0


e

i
∫

d

x


2
m

ℏ

2


(

E
−
V
(
x
)

)

+
θ


2
m

ℏ

2


(

E
−
V
(
x
)

)

4

 

(
25
)

{\displaystyle \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)}

{\displaystyle \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)}
מצד שני, אם הפאזה משתנה לאט יחסית בהשוואה למשרעת {\displaystyle (B\_{0}(x)=0)}, אזי מקבלים:


A

0

(
x
)
=
±

2
m

(

V
(
x
)
−
E

)

(
26
)

{\displaystyle A\_{0}(x)=\pm {\sqrt {2m\left(V(x)-E\right)}}\quad (26)}

{\displaystyle A_{0}(x)=\pm {\sqrt {2m\left(V(x)-E\right)}}\quad (26)}
המשוואה שלעיל, {\displaystyle \quad (26)}, תקפה רק כאשר האנרגיה הפוטנציאלית גדולה יותר מהאנרגיה הכללית {\displaystyle (E<V(x))}, כפי שקורה תמיד במנהור קוונטי.

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

{\displaystyle \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 הפתרונות המקורבים למשוואת שרדינגר, הן {\displaystyle \quad (25)} והן {\displaystyle \quad (27)}, נעשים סינגולריים בנקודת המפנה או נקודת ה"תפר", כאשר {\displaystyle E=V(x)} ולכן אינם תקפים בנקודה זו.

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

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

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

{\displaystyle {\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)}

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

{\displaystyle {\frac {\mathrm {d} ^{2}}{\mathrm {d} x^{2}}}\Psi (x)=U\_{1}\cdot (x-x\_{1})\cdot \Psi (x)\quad (29)}

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



Ψ
(
x
)
=

C

A

Ai

(

U

1

3

⋅
(
x
−

x

1

)

)

+

C

B

Bi

(

U

1

3

⋅
(
x
−

x

1

)

)

(
30
)

{\displaystyle \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)}

{\displaystyle \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 הפתרונות המקורבים מעל ומתחת גבעת הפוטנציאל נעשית על ידי מציאת היחס בין המקדמים השונים, {\displaystyle C\_{0},\theta } ו-{\displaystyle C\_{+},C\_{-}}. מאחר שפונקציות איירי בגבול האסימפטוטי, הופכות לפונקציות אקספוננטציאליות כגון: סינוס וקוסינוס, ניתן להביע את היחס בין המקדמים השונים דרך "פונקציות הקשר":

{\displaystyle C\_{+}=+{\frac {1}{2}}C\_{0}\cos {\left(\theta -{\frac {\pi }{4}}\right)}\quad (31)}



C

−

=
−

1
2

C

0

sin
⁡

(

θ
−

π
4

)

(
32
)

{\displaystyle C\_{-}=-{\frac {1}{2}}C\_{0}\sin {\left(\theta -{\frac {\pi }{4}}\right)}\quad (32)}

{\displaystyle 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. מסת"ב 0-07-004452-X
* Child, M. S. (1991). Semiclassical mechanics with molecular applications. Oxford: Clarendon Press. מסת"ב 0-19-855654-3.
* Griffiths, David J. (2004). Introduction to Quantum Mechanics (2nd ed.). Prentice Hall. מסת"ב 0-13-111892-7.
* Hall, B.C. (2013). Quantum Theory for Mathematicians. Springer.
* Liboff, Richard L. (2003). Introductory Quantum Mechanics (4th ed.). Addison-Wesley. מסת"ב 0-8053-8714-5.
* Olver, Frank William John (1974). Asymptotics and Special Functions. Academic Press. מסת"ב 0-12-525850-X.
* Razavy, Mohsen (2003). Quantum Theory of Tunneling. World Scientific. מסת"ב 981-238-019-1.
* Sakurai, J. J. (1993). Modern Quantum Mechanics. Addison-Wesley. מסת"ב 0-201-53929-2.

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

1. ↑  R.B. Dingle, **Asymptotic Expansions: Their Derivation and Interpretation**, Academic Press, 1973
2. ↑  Adrian E. Gill, **Atmosphere-ocean dynamics**, Academic Press, 1982, עמ' 297, מסת"ב 978-0-12-283522-3