مشكلة رابي

تتعلق مسألة رابي باستجابة الذرة لحقل كهربائي توافقي مطبق ، بتردد قريب جدًا من ترددها الطبيعي . وهي تقدم مثالًا بسيطًا وقابلًا للحل عمومًا لتفاعلات الضوء مع الذرة، وقد سُميت نسبةً إلى إيزيدور إسحاق رابي .

مشكلة رابي الكلاسيكية

في النهج الكلاسيكي، يمكن تمثيل مسألة رابي بحل المذبذب التوافقي المخمد المدفوع مع الجزء الكهربائي من قوة لورنتز كحد محرك:

x¨أ+2τ0x˙أ+ωأ2xأ=هـمهـ(ت،رأ)،{\displaystyle {\ddot {x}}_{a}+{\frac {2}{\tau _{0}}}{\dot {x}}_{a}+\omega _{a}^{2}x_{a}={\frac {e}{m}}E(t,\mathbf {r} _{a}),}

حيث يُفترض أن الذرة يمكن اعتبارها جسيمًا مشحونًا (بشحنة e ) يتذبذب حول موضع اتزانه حول ذرة متعادلة. هنا، x a هو مقدار تذبذبها اللحظي.ωأ{\displaystyle \omega _{a}}تردد تذبذبها الطبيعي، وτ0{\displaystyle \tau _{0}}عمره الطبيعي :

2τ0=2هـ2ωأ23مج3،{\displaystyle {\frac {2}{\tau _{0}}}={\frac {2e^{2}\omega _{a}^{2}}{3mc^{3}}},}

والتي تم حسابها بناءً على فقدان الطاقة للمذبذب ثنائي القطب من الإشعاع الكهرومغناطيسي.

ولتطبيق هذا على مسألة رابي، يفترض المرء أن المجال الكهربائي E متذبذب في الزمن وثابت في المكان:

هـ=هـ0[هـأناωت+هـ-أناωت]=2هـ0كوسωت،{\displaystyle E=E_{0}[e^{i\omega t}+e^{-i\omega t}]=2E_{0}\cos \omega t,}

ويتم تحليل x a إلى جزء u a متوافق في الطور مع المجال الكهربائي المحرك (المقابل للتشتت) وجزء v a غير متوافق في الطور (المقابل للامتصاص):

xأ=x0(uأكوسωت+vأالخطيئةωت).{\displaystyle x_{a}=x_{0}(u_{a}\cos \omega t+v_{a}\sin \omega t).}

هنا يُفترض أن x₀ ثابت ، ولكن يُسمح لـ uₐ و vₐ بالتغير مع الزمن. ومع ذلك، إذا كان النظام قريبًا جدًا من الرنين (ωωأ{\displaystyle \omega \approx \omega _{a}}إذاً، ستتغير هذه القيم ببطء مع مرور الوقت، ويمكننا أن نفترض أنu˙أωuأ{\displaystyle {\dot {u}}_{a}\ll \omega u_{a}}،v˙أωvأ{\displaystyle {\dot {v}}_{a}\ll \omega v_{a}}وu¨أω2uأ{\displaystyle {\ddot {u}}_{a}\ll \omega ^{2}u_{a}}،v¨أω2vأ{\displaystyle {\ddot {v}}_{a}\ll \omega ^{2}v_{a}}.

بناءً على هذه الافتراضات، يمكن إعادة كتابة معادلات قوة لورنتز للأجزاء المتوافقة في الطور وغير المتوافقة في الطور على النحو التالي:

u˙=-دلتاv-uتي،{\displaystyle {\dot {u}}=-\delta v-{\frac {u}{T}},}
v˙=دلتاu-vتي+κهـ0،{\displaystyle {\dot {v}}=\delta u-{\frac {v}{T}}+\kappa E_{0},}

حيث استبدلنا العمر الطبيعيτ0{\displaystyle \tau _{0}}مع عمر فعال أكثر عمومية T (والذي قد يشمل تفاعلات أخرى مثل التصادمات)، وقد تم إسقاط الرمز السفلي a لصالح الانحراف المُعرَّف حديثًادلتا=ω-ωأ{\displaystyle \delta =\omega -\omega _{a}}وهو ما يُستخدم بنفس الكفاءة لتمييز الذرات ذات الترددات الرنينية المختلفة. وأخيرًا، الثابت

κ =تعريف هـمωx0{\displaystyle \kappa \ {\stackrel {\text{def}}{=}}\ {\frac {e}{m\omega x_{0}}}}

تم تعريفها.

يمكن حل هذه المعادلات على النحو التالي:

u(ت؛دلتا)=[u0كوسدلتات-v0الخطيئةدلتات]هـ-ت/تي+κهـ00تدتالخطيئةدلتا(ت-ت)هـ-(ت-ت)/تي،{\displaystyle u(t;\delta )=[u_{0}\cos \delta t-v_{0}\sin \delta t]e^{-t/T}+\kappa E_{0}\int _{0}^{t}dt'\sin \delta (tt')e^{-(tt')/T},}
v(ت؛دلتا)=[u0كوسدلتات+v0الخطيئةدلتات]هـ-ت/تي-κهـ00تدتكوسدلتا(ت-ت)هـ-(ت-ت)/تي.{\displaystyle v(t;\delta )=[u_{0}\cos \delta t+v_{0}\sin \delta t]e^{-t/T}-\kappa E_{0}\int _{0}^{t}dt'\cos \delta (tt')e^{-(tt')/T}.}

بعد أن تتلاشى جميع التأثيرات العابرة ، يأخذ حل الحالة المستقرة الشكل البسيط التالي

xأ(ت)=هـمهـ0(هـأناωتωأ2-ω2+2أناω/تي+نسخة)،{\displaystyle x_{a}(t)={\frac {e}{m}}E_{0}\left({\frac {e^{i\omega t}}{\omega _{a}^{2}-\omega ^{2}+2i\omega /T}}+{\text{cc}}\right),}

حيث يرمز "cc" إلى المرافق المركب للحد المقابل.

ذرة ثنائية المستوى

النهج شبه الكلاسيكي

تُقدم مسألة رابي الكلاسيكية بعض النتائج الأساسية وصورة بسيطة للفهم للمسألة، ولكن من أجل فهم ظواهر مثل الانعكاس والانبعاث التلقائي وانزياح بلوخ-سيجرت ، فإن المعالجة الميكانيكية الكمية الكاملة ضرورية.

أبسط الطرق هي استخدام تقريب الذرة ثنائية المستوى ، حيث يُعامل مستويان فقط من مستويات طاقة الذرة المعنية. لا توجد ذرة في الواقع بمستويين فقط من الطاقة، ولكن يمكن التعامل مع الانتقال بين حالتين من حالات التفاعل الفائق الدقيق في الذرة، كتقريب أولي، كما لو أن هذين المستويين فقط موجودان، بافتراض أن المحرك ليس بعيدًا جدًا عن الرنين.

تكمن ميزة الذرة ثنائية المستوى في أن أي نظام ثنائي المستوى يتطور بنفس الطريقة الأساسية التي يتطور بها نظام الدوران 1/2 ، وفقًا لمعادلات بلوخ ، التي تحدد ديناميكيات متجه الدوران الزائف في مجال كهربائي:

u˙=-دلتاv،{\displaystyle {\dot {u}}=-\delta v,}
v˙=دلتاu+κهـw،{\displaystyle {\dot {v}}=\delta u+\kappa Ew,}
w˙=-κهـv،{\displaystyle {\dot {w}}=-\kappa Ev,}

حيث قمنا بتطبيق تقريب الموجة الدوارة عن طريق استبعاد الحدود ذات السرعة الزاوية العالية (وبالتالي تأثيرها الصغير على ديناميكيات الدوران الكلية على مدى فترات زمنية طويلة) وتحويلها إلى مجموعة من الإحداثيات تدور بترددω{\displaystyle \omega }.

ثمة تشابه واضح هنا بين هذه المعادلات وتلك التي حددت تطور مكونات التذبذب المتوافقة وغير المتوافقة في الحالة الكلاسيكية. لكن الآن، يوجد حد ثالث w ، يمكن تفسيره على أنه فرق عدد الذرات بين الحالة المثارة والحالة الأرضية (يتراوح من -1 للدلالة على الحالة الأرضية تمامًا إلى +1 للدلالة على الحالة المثارة تمامًا ). تجدر الإشارة إلى أنه في الحالة الكلاسيكية، كان هناك طيف طاقة متصل يمكن أن يشغله المذبذب الذري، بينما في الحالة الكمومية (كما افترضنا) لا يوجد سوى حالتين (ذاتيتين) محتملتين للمسألة.

يمكن أيضاً صياغة هذه المعادلات في شكل مصفوفة:

ددت[uvw]=[0-دلتا0دلتا0κهـ0-κهـ0][uvw].{\displaystyle {\frac {d}{dt}}{\begin{bmatrix}u\\v\\w\\\end{bmatrix}}={\begin{bmatrix}0&-\delta &0\\\delta &0&\kappa E\\0&-\kappa E&0\end{bmatrix}}{\begin{bmatrix}u\\v\\w\\\end{bmatrix}}.}

ومن الجدير بالذكر أن هذه المعادلات يمكن كتابتها كمعادلة ترنح متجهي :

دρدت=Ω×ρ،{\displaystyle {\frac {d\mathbf {\rho } }{dt}}=\mathbf {\Omega } \times \mathbf {\rho } ,}

أينρ=(u،v،w){\displaystyle \mathbf {\rho } =(u,v,w)}هو متجه الدوران الزائف، وΩ=(-κهـ،0،دلتا){\displaystyle \mathbf {\Omega } =(-\kappa E,0,\delta )}يعمل كعزم دوران فعال.

كما في السابق، يتم حل مشكلة رابي بافتراض أن المجال الكهربائي E متذبذب ذو مقدار ثابت E 0 :هـ=هـ0(هـأناωت+نسخة){\displaystyle E=E_{0}(e^{i\omega t}+{\text{c.c.}})}في هذه الحالة، يمكن إيجاد الحل بتطبيق دورتين متتاليتين على معادلة المصفوفة أعلاه، على النحو التالي:

[uvw]=[كوسχ0الخطيئةχ010-الخطيئةχ0كوسχ][uvw]{\displaystyle {\begin{bmatrix}u\\v\\w\end{bmatrix}}={\begin{bmatrix}\cos \chi &0&\sin \chi \\0&1&0\\-\sin \chi &0&\cos \chi \end{bmatrix}}{\begin{bmatrix}u'\\v'\\w'\end{bmatrix}}}

و

[uvw]=[1000كوسΩتالخطيئةΩت0-الخطيئةΩتكوسΩت][u"v"w"]،{\displaystyle {\begin{bmatrix}u'\\v'\\w'\end{bmatrix}}={\begin{bmatrix}1&0&0\\0&\cos \Omega t&\sin \Omega t\\0&-\sin \Omega t&\cos \Omega t\end{bmatrix}}{\begin{bmatrix}u''\\v''\\w''\end{bmatrix}},}

أين

لون برونزيχ=دلتاκهـ0،{\displaystyle \tan \chi ={\frac {\delta }{\kappa E_{0}}},}
Ω(دلتا)=دلتا2+(κهـ0)2.{\displaystyle \Omega (\delta )={\sqrt {\delta ^{2}+(\kappa E_{0})^{2}}}.}

هنا الترددΩ(دلتا){\displaystyle \Omega (\delta )}يُعرف هذا بتردد رابي المعمم ، والذي يُعطي معدل ترنح متجه الدوران الزائف حول المحور u ' المُحوَّل (المُعطى بتحويل الإحداثيات الأول أعلاه). على سبيل المثال، إذا كان المجال الكهربائي (أو الليزر ) في حالة رنين تام (بحيثدلتا=0{\displaystyle \delta =0}ثم سيدور متجه الدوران الزائف حول المحور u بمعدلκهـ0{\displaystyle \kappa E_{0}}إذا تم تسليط هذه النبضة (عند الرنين) على مجموعة من الذرات التي كانت جميعها في الأصل في حالتها الأرضية ( w = −1) لفترة زمنيةΔت=π/κهـ0{\displaystyle \Delta t=\pi /\kappa E_{0}}ثم بعد النبضة، ستكون جميع الذرات في حالتها المثارة ( w = +1) بسببπ{\displaystyle \pi }(أو دوران 180 درجة) حول المحور u . يُعرف هذا باسمπ{\displaystyle \pi }-نبضة ويكون لها نتيجة انعكاس كامل.

والنتيجة العامة معطاة بواسطة

[uvw]=[(κهـ0)2+دلتا2كوسΩتΩ2-دلتاΩالخطيئةΩت-دلتاκهـ0Ω2(1-كوسΩت)دلتاΩالخطيئةΩتكوسΩتκهـ0ΩالخطيئةΩتدلتاκهـ0Ω2(1-كوسΩت)-κهـ0ΩالخطيئةΩتدلتا2+(κهـ0)2كوسΩتΩ2][u0v0w0].{\displaystyle {\begin{bmatrix}u\\v\\w\end{bmatrix}}={\begin{bmatrix}{\frac {(\kappa E_{0})^{2}+\delta ^{2}\cos \Omega t}{\Omega ^{2}}}&-{\frac {\delta }{\Omega }}\sin {\Omega t}&-{\frac {\delta \kappa E_{0}}{\Omega ^{2}}}(1-\cos \Omega t)\\{\frac {\delta }{\Omega }}\sin \Omega t&\cos \Omega t&{\frac {\kappa E_{0}}{\Omega }}\sin \Omega t\\{\frac {\delta \kappa E_{0}}{\Omega ^{2}}}(1-\cos \Omega t)&-{\frac {\kappa E_{0}}{\Omega }}\sin {\Omega t}&{\frac {\delta ^{2}+(\kappa E_{0})^{2}\cos \Omega t}{\Omega ^{2}}}\end{bmatrix}}{\begin{bmatrix}u_{0}\\v_{0}\\w_{0}\end{bmatrix}}.}

يمكن تبسيط تعبير الانعكاس w بشكل كبير إذا افترضنا أن الذرة في البداية في حالتها الأرضية ( w₀ = -1 ) مع u₀ = v₀ = 0 ، وفي هذه الحالة

w(ت؛دلتا)=-1+2(κهـ0)2(κهـ0)2+دلتا2الخطيئة2(Ωت2).{\displaystyle w(t;\delta )=-1+{\frac {2(\kappa E_{0})^{2}}{(\kappa E_{0})^{2}+\delta ^{2}}}\sin ^{2}\left({\frac {\Omega t}{2}}\right).}

مسألة رابي في نظرية الاضطراب المعتمدة على الزمن

في النهج الكمومي، يمكن اعتبار القوة الدافعة الدورية بمثابة اضطراب دوري، وبالتالي، يمكن حل المشكلة باستخدام نظرية الاضطراب المعتمدة على الزمن، مع

ح(ت)=ح0+ح1(ت)،{\displaystyle H(t)=H^{0}+H^{1}(t),}

أينح0{\displaystyle H^{0}}هي دالة هاميلتونية مستقلة عن الزمن تُعطي الحالات الذاتية الأصلية، وح1(ت){\displaystyle H^{1}(t)}هو الاضطراب المعتمد على الزمن. افترض عند الزمنت{\displaystyle t}يمكننا توسيع نطاق الولاية كما

ϕ(ت)=ندن(ت)هـ-أناهـنت/|ن،{\displaystyle \phi (t)=\sum _{n}d_{n}(t)e^{-iE_{n}t/\hbar }|n\rangle ,}

أين|ن{\displaystyle |n\rangle }تمثل هذه القيم الحالات الذاتية للحالات غير المضطربة. بالنسبة لنظام غير مضطرب،دن(ت)=دن(0){\displaystyle d_{n}(t)=d_{n}(0)}هو ثابت. الآن، لنحسبدن(ت){\displaystyle d_{n}(t)}تحت تأثير اضطراب دوريح1(ت)=ح1هـ-أناωت{\displaystyle H^{1}(t)=H^{1}e^{-i\omega t}}تطبيق عامل التشغيلأنا/ت-ح0-ح1{\displaystyle i\hbar \partial /\partial t-H^{0}-H^{1}}على كلا جانبي المعادلة السابقة، يمكننا الحصول على

0=ن[أناد˙ن-ح1هـ-أناωتدن]هـ-أناهـن0ت/|ن،{\displaystyle 0=\sum _{n}[i\hbar {\dot {d}}_{n}-H^{1}e^{-i\omega t}d_{n}]e^{-iE_{n}^{0}t/\hbar }|n\rangle ,}

ثم اضرب طرفي المعادلة فيم|هـأناهـم0ت/{\displaystyle \langle m|e^{iE_{m}^{0}t/\hbar }}:

أناد˙م=نم|ح1|نهـأنا(ωمن-ω)تدن.{\displaystyle i\hbar {\dot {d}}_{m}=\sum _{n}\langle m|H^{1}|n\rangle e^{i(\omega _{mn}-\omega )t}d_{n}.}

عندما يكون تردد الإثارة عند الرنين بين حالتين|م{\displaystyle |m\rangle }و|ن{\displaystyle |n\rangle }، أيω=ωمن{\displaystyle \omega =\omega _{mn}}وبذلك، تصبح المسألة مسألة نمط طبيعي لنظام ثنائي المستوى، ومن السهل إيجاد ذلك.

دم،ن(ت)=دم،ن،+(0)هـأناΩت+دم،ن،-(0)هـ-أناΩت،{\displaystyle d_{m,n}(t)=d_{m,n,+}(0)e^{i\Omega t}+d_{m,n,-}(0)e^{-i\Omega t},}

أينΩ=م|ح1|ن.{\displaystyle \Omega ={\frac {\langle m|H^{1}|n\rangle }{\hbar }}.}

احتمالية التواجد في الحالة m عند الزمن t هي

Pم=دم(ت)*دم(ت)=دم،-2(0)+دم،+2(0)+2دم،-(0)دم،+(0)كوس(2Ωت).{\displaystyle P_{m}=d_{m}(t)^{*}d_{m}(t)=d_{m,-}^{2}(0)+d_{m,+}^{2}(0)+2d_{m,-}(0)d_{m,+}(0)\cos(2\Omega t).}

قيمةدم،±(0){\displaystyle d_{m,\pm }(0)}يعتمد ذلك على الحالة الأولية للنظام.

توصل رابي (1937) إلى حل دقيق لنظام ذي عزم مغزلي 1/2 في مجال مغناطيسي متذبذب. ويتضح من عمله أن تردد تذبذب رابي يتناسب طرديًا مع شدة المجال المغناطيسي المتذبذب.

نهج نظرية الحقل الكمومي

في منهج بلوخ، لا يتم تكميم المجال، ولا يتم تفسير التماسك الناتج ولا الرنين بشكل جيد.

هناك حاجة إلى عمل لنهج نظرية الحقول الكمومية، وخاصة نموذج جاينز-كامينغز .

انظر أيضاً

مراجع

  • ألين، ل؛ إيبرلي، جيه إتش (1987). الرنين البصري وذرات المستويين . نيويورك: دوفر. ISBN 978-0-486-65533-8. OCLC 17233252 .