הצגת המקור של הדף "תורת ההפרעות (מכניקת הקוונטים)"

'''תורת ההפרעות''' (בתרגומים ישנים מכונה גם '''טרידה''') היא שיטה לפתרון מקורב של בעיות ב[[מכניקת הקוונטים]], על ידי שימוש ב[[תורת ההפרעות]] המתמטית. עיקר השיטה הוא הפרדת הבעיה הפיזיקלית לשתי תתי-בעיות:<ref>{{קישור כללי|כתובת=https://phys.libretexts.org/Bookshelves/Quantum_Mechanics/Quantum_Mechanics_(Fowler)/09%3A_Perturbation_Theory/9.01%3A_Time-Independent_Perturbation_Theory|כותרת=9.1: Time-Independent Perturbation Theory|אתר=Physics LibreTexts|תאריך=2016-12-16|שפה=en|תאריך_וידוא=2023-04-07}}</ref>

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

על פי השיטה ניתן להרכיב את הפתרונות של הבעיה הכללית (כולל ההפרעה) באמצעות [[סופרפוזיציה]] של פתרונות של בעיית הבסיס.

==מבוא==

במסגרת תורת ההפרעות ניתן לנתח מערכות פיזיקלית עם [[המילטוניאן]] מהצורה:

<math>\hat{H}=\hat{H}_0+\lambda\hat{V}</math>

כך ש-<math>\lambda\ll1</math> ול-<math>\hat{H}_0</math> אוסף [[מצב קוונטי|מצבים עצמיים]] בדידים <math>|n^{(0)}\rang</math> עם אנרגיות עצמיות <math>E_n^{(0)}</math>. על פי השיטה ניתן למצוא אנרגיות עצמיות ומצבים עצמיים של ההמילטוניאן הכללי <math>\hat{H}</math> באמצעות הפתרונות הידועים של <math>\hat{H}_0</math>. פתרונות אלו יהיו מהצורה:

<math>|n\rang=|n^{(0)}\rang+\lambda|n^{(1)}\rang+\lambda^2|n^{(2)}\rang+\dots</math>

עם אנרגיות עצמיות מתאימות מהצורה:

<math>E_n=E_n^{(0)}+\lambda E_n^{(1)}+\lambda^2 E_n^{(2)}+\dots</math>

האיברים <math>E_n^{(k)}</math> ו-<math>|n^{(k)}\rang</math> נקראים התיקון מסדר <math>k</math> לאנרגיה ולמצב העצמי בהתאמה.

בעיות מתחום תורת ההפרעות ניתן לחלק לשני סוגים:
* '''תורת ההפרעות שאינה תלויה בזמן''': זהו המקרה בו ההפרעה <math>\hat{V}</math> בעלת תלות מרחבית בלבד ואינה תלויה בזמן. במקרה זה ישנו [[חוק שימור האנרגיה|שימור אנרגיה]] ועל כן ניתן לפתור את הבעיה כפתרון בעיית [[ערך עצמי|ערכים עצמיים]].
* '''תורת ההפרעות שתלויה בזמן''': זהו המקרה בו ההפרעה תלויה במרחב ובזמן, כלומר <math>\hat{V}\equiv \hat{V}(t)</math>. במקרה זה האנרגיה אינה נשמרת ולמעשה הפתרון מיוצג על-ידי מעברים משתנים בזמן בין המצבים העצמיים של בעיית הבסיס.

את תורת ההפרעות שאינה תלויה בזמן פיתח [[ארווין שרדינגר]] ב-[[1926]] כהמשך ישיר לעבודתו על [[משוואת שרדינגר]].<ref>{{צ-מאמר|מחבר=E. Schrödinger|שם=Quantisierung als Eigenwertproblem|כתב עת=Annalen der Physik|כרך=385|שנת הוצאה=1926|עמ=437–490|doi=10.1002/andp.19263851302|קישור=https://onlinelibrary.wiley.com/doi/10.1002/andp.19263851302}}</ref> את [[תורת ההפרעות התלויה בזמן]] פיתח [[פול דיראק]].<ref>{{צ-ספר|שם=Quantum mechanics|קישור=https://www.worldcat.org/oclc/42752492|מו"ל=Dover Publications|שנת הוצאה=1999|מקום הוצאה=Mineola, N.Y.|ISBN=0-486-40924-4|oclc=42752492|מחבר=Albert Messiah}}</ref>

==תורת ההפרעות שאינה תלויה בזמן==
=== ללא ניוון ===
כאשר פתרונות בעיית הבסיס הם ללא [[ניוון (פיזיקה)|ניוון]], ניתן לחשב את התיקון מסדר ראשון ושני לאנרגיה באופן הבא:

<math> E_n^{(1)} = \langle n^{(0)} | V | n^{(0)} \rangle </math>

<math>E_n^{(2)} = \sum_{k \ne n} \frac{|\langle k^{(0)}|V|n^{(0)} \rangle|^2} {E_n^{(0)} - E_k^{(0)}} </math>

את התיקון למצב העצמי מסדר ראשון ניתן לחשב באמצעות:

<math> |n^{(1)}\rang = \sum_{k \ne n} \frac{\langle k^{(0)}|V|n^{(0)} \rangle}{E_n^{(0)} - E_k^{(0)}} |k^{(0)}\rang </math>

פתרונות אלו נכונים אך ורק במידה ואין ניוון. במקרה ויש ניוון המכנה בתיקון מסדר שני לאנרגיה ומסדר ראשון למצב מתאפס וישנה [[חילוק באפס|חלוקה באפס]]. במקרה זה נדרשת שיטה אחרת.

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

עבור אוסף מצבים עצמיים <math>|1\rang,|2\rang,\dots,|k\rang</math> בעלי אותה אנרגיה עצמית <math>E_n</math> בבעיית הבסיס, ניתן לצמצם את ההפרעה <math>\hat{V}</math> למרחב הנפרס על ידי מצבים אלו ולהציג אותה בייצוג [[מטריצה|מטריצי]]:

<math>V=\left(\begin{matrix}
V_{11} & V_{12} & \dots &V_{1k} \\
V_{21} & V_{22} & \dots & V_{2k} \\
\vdots & \vdots & \ddots & \vdots \\
V_{k1} & V_{k2} & \dots & V_{kk}
\end{matrix}\right)</math>

כך ש:

<math>V_{ij}=\lang i|\hat{V}|j\rang</math>

עבור וקטור עצמי <math>(\alpha_1,\dots,\alpha_k)</math> וערך עצמי מתאים <math>\Delta E</math> של המטריצה <math>V</math>, מתקבל מצב עצמי ואנרגיה עצמית מתוקנים כך שהאנרגיה המתוקנת היא <math>E + \lambda\Delta E</math> והמצב המתוקן הוא:

<math>|\psi\rang=\sum_{i=1}^k{\alpha_i|i\rang}</math>

ובכך ניתן לקבל <math>k</math> מצבים עצמיים מתוקנים עם אנרגיות מתוקנות.

=== דוגמאות ===
==== אלקטרונים כמעט חופשיים בגביש חד-ממדי ====
ניתן לתאר [[גביש]] חד־ממדי כמערכת עם המילטוניאן מהצורה <math>H=\frac{p^2}{2m}+V</math>, כאשר <math>V</math> הוא הפרעה מחזורית עם מחזור מאורך <math>a</math>. כלומר, לכל <math>x\in\mathbb{R}</math> מתקיים <math>V(x+a)=V(x)</math>. במקרה זה ההמילטוניאן הבלתי מופרע הוא <math>\frac{p^2}{2m}</math> והמצבים העצמיים שלו הם ה[[חלקיק חופשי|אלקטרונים החופשיים]] <math>|k\rang</math> עם האנרגיות <math>E_k=\frac{\hbar^2k^2}{2m}</math> כאשר <math>k\in\mathbb{R}</math> הוא [[מספר גל|מספר הגל]].

ניתן לתאר את הפוטנציאל המחזורי באמצעות [[טור פורייה]]:

<math>V(x)=\sum_{n\in\mathbb{Z}}{\tilde{V}_ne^{nKx}}</math> כאשר <math>K=\frac{2\pi}{a}</math> ו-<math>\tilde{V}_{-n}=\tilde{V}_n^*</math> (זאת מכיוון ש-<math>V</math> חייב להיות [[אופרטור הרמיטי|הרמיטי]])

מאחר שהוספת קבוע להמילטוניאן שקול להזזת [[מופע|פאזה]], מה שאינו משנה את סטטיסטיקת מיקום ה[[אלקטרון|אלקטרונים]], ניתן להניח בלי הגבלת הכלליות כי <math>\tilde{V}_0=0</math>.

בבעיית הבסיס אין ניוון באנרגיה, לכן בעזרת תורת ההפרעות מתקבל כי התיקון הראשון לאנרגיה מתאפס והתיקון השני לאנרגיה הוא:

<math>E_k^{(2)}=\sum_{0\ne n\in\mathbb{Z}}{\frac{|\tilde{V}_n|^2}{E_k-E_{k+nK}}}</math>

פתרון זה נקרא [[מודל אלקטרונים כמעט חופשיים]] והוא ניתן להרחבה גם לשלושה ממדים.

==== אפקט שטארק באטום מימן ====
[[אפקט שטארק]] באטום [[מימן]] מתקבל כאשר המימן מצוי ב[[שדה חשמלי]]. במקרה זה ההמילטוניאן הוא מהצורה <math>H=H_0+eEz</math> כאשר <math>H_0</math> הוא ההמילטוניאן של [[אטום המימן]] ו-<math>V=eEz</math> מתאר שדה חשמלי מגודל <math>E</math> בכיוון ציר <math>z</math> (<math>e</math> הוא [[מטען חשמלי|מטען]] האלקטרון).

בבעיה הבלתי מופרעת של אטום המימן ישנו ניוון מגודל <math>n^2</math> לכל רמת אנרגיה <math>n</math> ולכן יש להשתמש בתורת ההפרעות במקרה המנוון.

עבור <math>n=2</math>, המצבים היחידים המקיימים אינטראקציה זה עם זה הם <math>|200\rang</math> ו-<math>|210\rang</math> עם <math>\lang200|V|210\rang=3eEa_0</math> כאשר <math>a_0</math> הוא [[רדיוס בוהר]].

מכל זה, המצבים <math>|211\rang</math> ו-<math>|21{-1}\rang</math> נשארים ללא שינוי, ומתקבלים מצבים מלוכסנים (לפי ההפרעה), <math>\frac{\sqrt{2}}{2}\left(|200\rang\pm|210\rang\right)</math>, בעלי אנרגיות <math>E_2\pm 3eEa_0</math> בהתאמה.

==תורת ההפרעות התלויה בזמן==
=== מעבר בין מצבים עצמיים של ההמילטוניאן הבלתי מופרע ===
לרוב בבעיות של תורת ההפרעות התלויה בזמן יש למצוא את הסיכוי שחלקיק ימצא במצב קוונטי סופי <math>|f\rang</math> בהינתן מצב התחלתי <math>|b\rang</math>, כאשר <math>|b\rang</math> ו-<math>|f\rang</math> שניהם מצבים עצמיים של ההמילטוניאן הבלתי מופרע. ניתן לעשות זאת על-ידי הנוסחה:

<math>\ c_f^{(1)} (t) = -\frac{i}{\hbar} \lambda \int_{t_0}^{t}{ V_{fb}(t) e^{i \omega_{fb} t} dt}</math>

כאשר <math>V_{fb}(t)=\lang f|V(t)|b\rang</math> ו-<math>\omega_{fb} =\frac{E_f - E_b}{\hbar}</math>, וכן <math>t_0</math> הוא הזמן שממנו מתחילה ההפרעה. הסיכוי שחלקיק ימצא במצב קוונטי <math>|f\rang</math> בזמן <math>t</math> כלשהו הוא <math>|c_f^{(1)}(t)|^2</math>.

פיתוח זה תקף רק עבור <math>\lambda \ll  1</math> ו-<math>\ | c_f |^2 \ll 1</math>.

=== פיתוח מתמטי ===
בתורת ההפרעות התלויה בזמן יש לפתור את משוואת שרדינגר מהצורה:

<math>i\hbar\frac{\partial}{\partial t}\psi(t)=(H_0+\lambda V(t))|\psi(t)\rang</math>.

מנחשים פתרון מהצורה:

<math>\ | \psi (t) \rang = \sum_{n} c_n (t) e^{-i E_n t / \hbar} | n \rang </math>

כאשר אף כאן המצבים הקוונטיים <math>|n\rang</math> הם מצבים עצמיים של ההמילטוניאן הבלתי מופרע עם האנרגיה העצמית <math>E_n</math>. בגבול <math>\lambda \to 0</math> המקדמים <math>c_n</math> מתנוונים לקבועים בזמן. מציבים אותה במשוואה ומטילים על מצב עצמי כלשהו, ובכך מקבלים מערכת [[משוואה דיפרנציאלית|משוואות דיפרנציאליות]] מצומדות עבור המקדמים <math>c_n(t)</math>.
: <math>\ i \hbar \frac{\partial}{\partial t} c_n (t) = \lambda \sum_k c_k(t) e^{i(E_n - E_k)t / \hbar} \lang n | V(t) | k \rang </math>
על-ידי סימון <math>\ \omega_{nk} = ( E_n - E_k )/ \hbar</math> ו-<math>V_{nk} = \lang n | V(t) | k \rang </math>, מתקבלת הצורה הסופית:
: <math>\ i \hbar \frac{\partial}{\partial t} c_n (t) = \lambda \sum_k  e^{i\omega_{nk}t } V_{nk} c_k(t) </math>

אפשר לפתור את המשוואה על ידי לכסון מטריצת ההפרעה המוכפלת ב"תדירויות ההפרשים". לרוב, הפתרון הכללי מסובך לחישוב.

=== תמונת האינטראקציה של דיראק ===
בתמונת האינטראקציה של [[פול דיראק|דיראק]] מעבירים את התלות הטריוויאלית בזמן, זו הנובעת מ-H<sub>0</sub>, ה[[המילטוניאן]] הלא מופרע, מפונקציות הגל לאופרטורים. אזי:

: <math>\ | \psi (t) \rang _{Sc} \to | \psi (t) \rang _I = e^{i H_0 t / \hbar} | \psi \rang _{Sc} </math>
: <math>\ A_{Sc} \to A_I = e^{i H_0 t / \hbar} A_{Sc} e^{-i H_0 t / \hbar}</math>

ואז הנוסחה הכללית למצב המופרע נהפכת ל

: <math>\ | \psi (t) \rang _I = \sum_{n} c_n (t) | n \rang </math>

כאשר <math>\ \dot{c_n} \ne 0 \iff \lambda V \ne 0</math>. משוואת התנועה הופכת להיות

: <math>\  i\hbar \frac{\partial}{\partial t} | \psi \rang _I =  \lambda V_I(t) | \psi \rang _I </math>

ופתרונה הכללי הוא

: <math>\ | \psi (t) \rang _I = e^{-\frac{i}{\hbar} \lambda \int_{0}^{t}{ V_I(t') dt'}} | \psi (0) \rang _I</math>

כאשר אנו מניחים ש <math>\ \forall t_1,t_2 \ : \ \left[ V_I(t_1) , V_I(t_2) \right] = 0</math>.

נסמן <math>\ U_I = e^{-\frac{i}{\hbar} \lambda \int_{0}^{t}{ V_I(t') dt'}}</math> את האופרטור של הקידום בזמן של ההפרעה. זהו אופרטור [[אוניטריות|אוניטרי]].

אזי ההסתברות לעבור ממצב עצמי b למצב עצמי f נתונה על ידי

: <math>\ \mbox{Prob}_{b \to f} = | \lang f | U_I ( \lambda) | b \rang |^2</math>

וכדי לפתור זאת, מפתחים את [[אופרטור]] הקידום בזמן של ההפרעה לפי <math>\lambda</math> במה שנקרא [[סדרת דייסון]]

: <math>\ U_I(t) = e^{-\frac{i}{\hbar} \lambda \int_{0}^{t}{ V_I(t') dt'}} = 1 + \lambda U_I^{(1)} + \lambda^2 U_I^{(2)} + \cdots </math>.

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

: <math>\ i \hbar \frac{\partial}{\partial t} U_I(t) = \lambda V_I(t) U_I(t)</math>

מציבים ופותרים לכל סדר, וכך למשל מקבלים את התיקון הראשון

: <math>\ U_I^{(1)}(t) = -\frac{i}{\hbar} \int_0^t{ V_I(t') dt' }</math>

את התיקון השני

: <math>\ U_I^{(2)}(t) = \left( -\frac{i}{\hbar} \right) ^2  \int_0^t{dt' V_I(t') \int_0^{t'}{dt'' V_I(t'')}  }</math>

כאשר תמיד מתקיים <math>\ 0 \le t'' < t' < t</math>.

התיקון הראשון מייצג מעבר בין b ל-f בקפיצה אחת, התיקון השני מייצג מעבר בין b ל-f ב-2 קפיצות דרך רמת ביניים אחת, השני מייצג מעבר בין b ל-f ב-3 קפיצות דרך 2 רמות ביניים, וכן הלאה.

=== דוגמאות ===
==== תהליך אדיאבטי ====
ב[[תהליך אדיאבטי]] ההפרעה משתנה לאט מאוד כתלות בזמן. במקרה זה, המצב העצמי ה-<math>n</math>-י הבלתי מופרע משתנה בצורה [[רציפות|רציפה]] למצב ה-<math>n</math>-י המופרע המקביל לו והאנרגיה משתנה כך גם כן.

ניתן לתאר תהליך אדיאבטי באמצעות הפרעה מהצורה:

: <math>
V(t) = \left\{ \begin{matrix}
0 & \mbox{if }t = - \infty \\
\lambda V e^{t / \tau } & \mbox{if }t < 0 \\
\lambda V & \mbox{if }t \ge 0
\end{matrix} \right.
</math>

כאשר השינוי בהפרעה איטי מאוד, אזי:

: <math>\ | n(t \to + \infty) \rang \to | n' \rang</math>

כאשר <math>\ ( H + \lambda V ) | n' \rang = E_n' | n' \rang </math> ואותם אפשר למצוא באמצעות תורת ההפרעות שאינה תלויה בזמן.

הפתרון במקרה זה הוא

: <math>|\psi_n(t)\rang = e^{-\frac{i}{\hbar} \int_{t_0}^{t}{E_n(t') dt'}} | n(t) \rang</math>

וזה פתרון רק אם השינוי ב[[פונקציית גל|פונקציית הגל]] של המצב העצמי הוא מאוד איטי, כלומר, בגבול <math>\tau \to \infty</math>.

==== הפרעה מחזורית ====
נניח הפרעה מחזורית <math>\ V(t) = V \sin ( \omega t)</math> ב[[תדירות]] <math>\omega</math>. כמו בנוסחה הכללית ניתן לסמן <math>\ \omega_{fb} = \left( E_f - E_b \right) / \hbar</math>. במקרה זה, האמפליטודה למעבר ממצב b למצב f היא

: <math>\begin{align}c_f^{(1)} = -\frac{i}{\hbar}\lambda V_{fb} \int_{0}^{t}{ e^{i \omega_{fb}t'} \sin (\omega t') dt' }\\
= \frac{i \lambda V_{fb}}{2 \hbar} \left( \frac{1 - e^{i( \omega_{fb} + \omega )t}}{\omega_{fb} + \omega} + \frac{1 - e^{i( \omega_{fb} - \omega )t}}{ \omega_{fb} - \omega} \right) \\
= \frac{i \lambda V_{fb}}{2 \hbar} ( A_+ + A_- )\end{align} </math>

האיבר <math>\ A_-</math> מייצג בליעה של [[פוטון]] ועלייה ברמת האנרגיה, ואילו האיבר <math>\ A_+</math> מייצג פליטה של פוטון וירידה ברמת האנרגיה. כל איבר נהייה דומיננטי כאשר התדירות קרובות להפרש האנרגיות המתאים. אם מניחים כי <math>\ \omega \approx \omega_{fb}</math> מתקבל ש-<math>\ A_-</math> דומיננטי. במקרה זה, ה[[הסתברות]] למעבר רמות מ-b ל-f היא:

: <math> \begin{align}\mbox{Prob}_{b \to f} = | c_f |^2 \\
= \frac{\lambda^2 |V_fb|^2}{4 \hbar^2} | A_- |^2 \\
= \frac{\lambda^2 |V_fb|^2}{4 \hbar^2} \left| e^{i \frac{\omega_{fb} - \omega}{2} t} \right| \cdot \left| \frac{\sin \left( ( \omega_{fb} - \omega ) t/2 \right) }{ (\omega_{fb} - \omega )/2 } \right| ^2\end{align}  </math>

[[קובץ:Sinc-resonance01.png|ממוזער|250x250 פיקסלים|הסתברות למעבר רמה כתלות בתדירות ההפרעה. ההסתברות מתנהגת כמו פונקציית [[sinc]] סביב הפרש האנרגיות בין הרמות]]
ניתן לראות שההסתברות מתנהגת כמו פונקציית <math>\mbox{sinc}(\zeta) = \sin (\zeta) / \zeta</math> בריבוע כאשר <math>\zeta =( \omega_{fb} - \omega )t/2</math>. רק כאשר תדירות ההפרעה מקיימת <math>\ \omega_{fb} - \frac{2 \pi}{t} < \omega < \omega_{fb} + \frac{2 \pi}{t} </math> יש הסתברות משמעותית למעבר רמות. תופעה זאת נקראת [[תהודה|רזוננס]], שבה המערכת מגיבה רק לתדירות אחת. רוחב עקומת התהודה הוא <math>\ 4 \pi / t</math> ועבור <math>\ \omega \to \omega_{fb}</math> מקבלים שההסתברות למעבר רמות פרופורציונית למשך ההפרעה:

: <math>\mbox{Prob}_{b \to f} = \frac{\lambda^2 |V_{fb}|^2}{4 \hbar^2} |t|^2 </math>

הניתוח לעיל תקף עבור <math>\ \frac{4 \pi}{\omega_{fb}} \ll t \ll \frac{2 \hbar}{\lambda |V_{fb}|}</math>, כלומר:

# על המערכת לדגום מחזורים רבים של ההפרעה. תנאי זה נובע מרוחב עקומת התהודה.
# התנאי השני נובע מדרישת תקפות הפיתוח של תוה"פ, ש <math>\ |c_f|^2 \ll 1</math>. בשביל שיהיה t מספיק גדול שעבורו התנאי הראשון יתמלא, על ההפרעה להיות קטנה מאוד מהפרש האנרגיות, כלומר: <math>\ |V_{fb}| \ll \hbar \omega_{fb} = E_f - E_b</math>.

==קישורים חיצוניים==
* {{MathWorld}}
* {{דף שער בספרייה הלאומית|987007538619405171|הפרעות (דינמיקה קוונטית)}}

==הערות שוליים==
{{הערות שוליים|יישור=שמאל}}

{{בקרת זהויות}}

[[קטגוריה:ערכים שבהם תבנית בריטניקה אינה מתאימה]]
[[קטגוריה:אופטימיזציה רציפה]]
[[קטגוריה:מכניקת הקוונטים]]

[[de:Störungstheorie]]
[[ko:건드림이론]]