المذبذب التوافقي الكمي

بعض مسارات المذبذب التوافقي وفقًا لقوانين نيوتن للميكانيكا الكلاسيكية (أ-ب)، ووفقًا لمعادلة شرودنغر في ميكانيكا الكم (ج-ح). في أ-ب، يتذبذب الجسيم (المُمثَّل بكرة متصلة بنابض ) ذهابًا وإيابًا. في ج-ح، تُعرض بعض حلول معادلة شرودنغر، حيث يُمثل المحور الأفقي الموضع، والمحور الرأسي الجزء الحقيقي (باللون الأزرق) أو الجزء التخيلي (باللون الأحمر) من دالة الموجة . ج، د، هـ، و، وليس ز، ح، هي حالات طاقة ذاتية . ح هي حالة متماسكة - حالة كمومية تُقارب المسار الكلاسيكي.

المذبذب التوافقي الكمومي هو النظير الكمومي للمذبذب التوافقي الكلاسيكي . ولأن أي جهد أملس يمكن تقريبه عادةً بجهد توافقي بالقرب من نقطة اتزان مستقرة ، فهو يُعدّ من أهم الأنظمة النموذجية في ميكانيكا الكم. علاوة على ذلك، فهو من الأنظمة الكمومية القليلة التي يُعرف لها حل تحليلي دقيق. [ 1 ] [ 2 ] [ 3 ] [ 4 ]

مذبذب توافقي أحادي البعد

الهاميلتوني وحالات الطاقة الذاتية

تمثيلات دالة الموجة لأول ثمانية حالات ذاتية مقيدة، n = 0 إلى 7. يوضح المحور الأفقي الموضع x .
كثافات الاحتمال المقابلة.

هاميلتونيان الجسيم هو :ح^=ص^22م+12كx^2=ص^22م+12مω2x^2،{\displaystyle {\hat {H}}={\frac {{\hat {p}}^{2}}{2m}}+{\frac {1}{2}}k{\hat {x}}^{2}={\frac {{\hat {p}}^{2}}{2m}}+{\frac {1}{2}}m\omega ^{2}{\hat {x}}^{2}\,,} حيث m هي كتلة الجسيم، و k هو ثابت القوة،ω=ك/م{\textstyle \omega ={\sqrt {k/m}}}يمثل التردد الزاوي للمذبذب،x^{\displaystyle {\hat {x}}}هو عامل الموضع (المعطى بواسطة x في أساس الإحداثيات)، وص^{\displaystyle {\hat {p}}} هو مؤثر الزخم (المعطى بواسطةص^=-أنا/x{\displaystyle {\hat {p}}=-i\hbar \,\partial /\partial x}في أساس الإحداثيات). يمثل الحد الأول في الهاميلتوني الطاقة الحركية للجسيم، ويمثل الحد الثاني طاقته الكامنة، كما هو الحال في قانون هوك . [ 5 ]

معادلة شرودنغر المستقلة عن الزمن (TISE) هي، ح^|ψ=هـ|ψ ،{\displaystyle {\hat {H}}\left|\psi \right\rangle =E\left|\psi \right\rangle ~,} أينهـ{\displaystyle E}يشير إلى عدد حقيقي (يجب تحديده) يحدد مستوى طاقة مستقل عن الزمن ، أو قيمة ذاتية ، والحل |ψ{\displaystyle |\psi \rangle }يشير إلى الحالة الذاتية للطاقة لهذا المستوى . [ 6 ]

ثم قم بحل المعادلة التفاضلية التي تمثل مسألة القيم الذاتية هذه في أساس الإحداثيات، وذلك بالنسبة للدالة الموجية.x|ψ=ψ(x){\displaystyle \langle x|\psi \rangle =\psi (x)}باستخدام طريقة طيفية ، تبين وجود مجموعة من الحلول. في هذا الأساس، تُعتبر هذه الحلول دوال هيرميت ، [ 7 ] [ 8 ]ψن(x)=12نن!(مωπ)1/4هـ-مωx22حن(مωx)،ن=0،1،2،....{\displaystyle \psi _{n}(x)={\frac {1}{\sqrt {2^{n}\,n!}}}\left({\frac {m\omega }{\pi \hbar }}\right)^{1/4}e^{-{\frac {m\omega x^{2}}{2\hbar }}}H_{n}{\left({\sqrt {\frac {m\omega }{\hbar }}}x\right)},\qquad n=0,1,2,\ldots .}

الوظائفحن{\displaystyle H_{n}}هي كثيرات حدود هيرميت الفيزيائية ، حن(z)=(-1)ن هـz2دندzن(هـ-z2).{\displaystyle H_{n}(z)=(-1)^{n}~e^{z^{2}}{\frac {d^{n}}{dz^{n}}}\left(e^{-z^{2}}\right).}

مستويات الطاقة المقابلة هي [ 9 ]هـن=ω(ن+12).{\displaystyle E_{n}=\hbar \omega {\bigl (}n+{\tfrac {1}{2}}{\bigr )}.}يمكن استخلاص القيم المتوقعة للموقع والزخم، بالإضافة إلى تباين كل متغير، من دالة الموجة لفهم سلوك متجهات الطاقة الذاتية. وقد تبين أنهاx^=0{\textstyle \langle {\hat {x}}\rangle =0}وص^=0{\textstyle \langle {\hat {p}}\rangle =0}بسبب تناظر المسألة، بينما:

x2=(2ن+1)2مω=σx2ص2=(2ن+1)مω2=σص2{\displaystyle {\begin{aligned}\left\langle x^{2}\right\rangle &=(2n+1){\frac {\hbar }{2m\omega }}=\sigma _{x}^{2}\\[1ex]\left\langle p^{2}\right\rangle &=(2n+1){\frac {m\hbar \omega }{2}}=\sigma _{p}^{2}\end{aligned}}}

لوحظ ازدياد التباين في كل من الموضع والزخم مع ارتفاع مستويات الطاقة. ويبلغ أدنى مستوى طاقة قيمةً قدرهاσxσص=2{\textstyle \sigma _{x}\sigma _{p}={\frac {\hbar }{2}}}وهي قيمتها الدنيا بسبب علاقة عدم اليقين، وتتوافق أيضًا مع دالة موجية غاوسية. [ 10 ]

يُعدّ طيف الطاقة هذا جديرًا بالملاحظة لأربعة أسباب. أولًا، الطاقات مُكمّمة، أي أن قيم الطاقة الممكنة هي قيم منفصلة فقط (مضاعفات صحيحة زائد نصف لـ ħω )؛ وهذه سمة عامة للأنظمة الكمومية عندما يكون الجسيم محصورًا. ثانيًا، مستويات الطاقة المنفصلة هذه متساوية التباعد، على عكس نموذج بور للذرة، أو الجسيم في صندوق . ثالثًا، أدنى طاقة يُمكن بلوغها (طاقة الحالة n = 0 ، وتُسمى الحالة الأرضية ) لا تُساوي الحد الأدنى لبئر الجهد، بل تساوي ħω /2 فوقه؛ وتُسمى هذه طاقة النقطة الصفرية . وبسبب طاقة النقطة الصفرية، فإن موضع وزخم المذبذب في الحالة الأرضية ليسا ثابتين (كما هو الحال في المذبذب الكلاسيكي)، بل لهما نطاق تباين صغير، وفقًا لمبدأ هايزنبرغ للشك . رابعًا، مستويات الطاقة غير متماثلة، مما يعني أن كل قيمة ذاتية مرتبطة بحل واحد فقط (حالة واحدة). [ 11 ]

تتركز كثافة احتمالية الحالة الأرضية عند نقطة الأصل، مما يعني أن الجسيم يقضي معظم وقته في قاع بئر الجهد، كما هو متوقع لحالة ذات طاقة منخفضة. ومع ازدياد الطاقة، تبلغ كثافة الاحتمالية ذروتها عند "نقاط التحول" الكلاسيكية، حيث تتطابق طاقة الحالة مع طاقة الوضع. (انظر المناقشة أدناه حول الحالات المثارة بشدة). يتوافق هذا مع المذبذب التوافقي الكلاسيكي، حيث يقضي الجسيم معظم وقته (وبالتالي يكون وجوده أكثر احتمالاً) بالقرب من نقاط التحول، حيث يتحرك بأبطأ سرعة. وبذلك يتحقق مبدأ التوافق . علاوة على ذلك، تتذبذب حزم موجية خاصة غير مشتتة، ذات أدنى قدر من عدم اليقين، تُسمى الحالات المتماسكة ، بشكل مشابه جدًا للأجسام الكلاسيكية، كما هو موضح في الشكل؛ وهي ليست حالات ذاتية للهاميلتوني.

طريقة عامل السلم

كثافات الاحتمال | ψn ( x ) | ² للحالات الذاتية المقيدة، بدءًا من الحالة الأرضية ( n = 0) في الأسفل، وتزداد الطاقة تدريجيًا باتجاه الأعلى. يُظهر المحور الأفقي الموضع x ، وتمثل الألوان الأكثر سطوعًا كثافات احتمال أعلى.

تتيح طريقة " المؤثر السلمي "، التي طورها بول ديراك ، استخراج القيم الذاتية للطاقة دون الحاجة إلى حل المعادلة التفاضلية مباشرةً. [ 12 ] وهي قابلة للتعميم على مسائل أكثر تعقيدًا، لا سيما في نظرية الحقل الكمومي . وباتباع هذا النهج، نُعرّف المؤثرات.أ^{\displaystyle {\hat {a}}}ومرافقهأ^{\displaystyle {\hat {a}}^{\dagger }}، أ^=مω2(x^+أنامωص^)أ^=مω2(x^-أنامωص^){\displaystyle {\begin{aligned}{\hat {a}}&={\sqrt {m\omega \over 2\hbar }}\left({\hat {x}}+{i \over m\omega }{\hat {p}}\right)\\{\hat {a}}^{\dagger }&={\sqrt {m\omega \over 2\hbar }}\left({\hat {x}}-{i \over m\omega }{\hat {p}}\right)\end{aligned}}}لاحظ أن هذه المؤثرات، من الناحية الكلاسيكية، هي بالضبط مولدات الدوران المعياري في فضاء الطور لـx{\displaystyle x}ومدxدت{\displaystyle m{\frac {dx}{dt}}}أي أنها تصف التطور الأمامي والخلفي في الزمن لمذبذب توافقي كلاسيكي.

تؤدي هذه العوامل إلى التمثيل التالي لـx^{\displaystyle {\hat {x}}}و ص^{\displaystyle {\hat {p}}}، x^=2مω(أ^+أ^)ص^=أنامω2(أ^-أ^) .{\displaystyle {\begin{aligned}{\hat {x}}&={\sqrt {\frac {\hbar }{2m\omega }}}({\hat {a}}^{\dagger }+{\hat {a}})\\{\hat {p}}&=i{\sqrt {\frac {\hbar m\omega }{2}}}({\hat {a}}^{\dagger }-{\hat {a}})~.\end{aligned}}}

المؤثر a ليس هيرميتيًا ، لأنه هو ومرافقه a غير متساويين. وتعطي حالات الطاقة الذاتية | n⟩ ، عند تطبيق هذه المؤثرات السلمية عليها، أ^|ن=ن+1|ن+1أ^|ن=ن|ن-1.{\displaystyle {\begin{aligned}{\hat {a}}^{\dagger }|n\rangle &={\sqrt {n+1}}|n+1\rangle \\{\hat {a}}|n\rangle &={\sqrt {n}}|n-1\rangle .\end{aligned}}}

انطلاقاً من العلاقات المذكورة أعلاه، يمكننا أيضاً تعريف عامل عددي N ، والذي يتمتع بالخاصية التالية: شمال^=أ^أ^شمال^|ن=ن|ن.{\displaystyle {\begin{aligned}{\hat {N}}&={\hat {a}}^{\dagger }{\hat {a}}\\{\hat {N}}\left|n\right\rangle &=n\left|n\right\rangle .\end{aligned}}}

يمكن الحصول على المبدلات التالية بسهولة عن طريق استبدال علاقة التبديل الأساسية ، [أ^،أ^]=1،[شمال^،أ^]=أ^،[شمال^،أ^]=-أ^،{\displaystyle [{\hat {a}},{\hat {a}}^{\dagger }]=1,\qquad [{\hat {N}},{\hat {a}}^{\dagger }]={\hat {a}}^{\dagger },\qquad [{\hat {N}},{\hat {a}}]=-{\hat {a}},}

ويمكن التعبير عن مؤثر هاميلتون على النحو التالي: ح^=ω(شمال^+12)،{\displaystyle {\hat {H}}=\hbar \omega \left({\hat {N}}+{\frac {1}{2}}\right),}

إذن، الحالات الذاتية لـشمال^{\displaystyle {\hat {N}}} وهي أيضًا حالات الطاقة الذاتية. ولإثبات ذلك، يمكننا تطبيقح^{\displaystyle {\hat {H}}}إلى حالة رقمية|ن{\displaystyle |n\rangle }:

ح^|ن=ω(شمال^+12)|ن.{\displaystyle {\hat {H}}|n\rangle =\hbar \omega \left({\hat {N}}+{\frac {1}{2}}\right)|n\rangle .}

باستخدام خاصية عامل العددشمال^{\displaystyle {\hat {N}}}:

شمال^|ن=ن|ن،{\displaystyle {\hat {N}}|n\rangle =n|n\rangle ,}

نحصل على:

ح^|ن=ω(ن+12)|ن.{\displaystyle {\hat {H}}|n\rangle =\hbar \omega \left(n+{\frac {1}{2}}\right)|n\rangle .}

وبالتالي، بما أن|ن{\displaystyle |n\rangle }يحل هذا الأمر معادلة TISE لمؤثر هاميلتونيح^{\displaystyle {\hat {H}}}، وهي أيضاً إحدى حالاتها الذاتية ذات القيمة الذاتية المقابلة:

هـن=ω(ن+12).{\displaystyle E_{n}=\hbar \omega \left(n+{\frac {1}{2}}\right).}

وهو المطلوب إثباته.

خاصية التبديل تعطي شمال^أ^|ن=(أ^شمال^+[شمال^،أ^])|ن=(أ^شمال^+أ^)|ن=(ن+1)أ^|ن،{\displaystyle {\begin{aligned}{\hat {N}}{\hat {a}}^{\dagger }|n\rangle &=\left({\hat {a}}^{\dagger }{\hat {N}}+[{\hat {N}},{\hat {a}}^{\dagger }]\right)|n\rangle \\&=\left({\hat {a}}^{\dagger }{\hat {N}}+{\hat {a}}^{\dagger }\right)|n\rangle \\&=(n+1){\hat {a}}^{\dagger }|n\rangle ,\end{aligned}}}

وبالمثل، شمال^أ^|ن=(ن-1)أ^|ن.{\displaystyle {\hat {N}}{\hat {a}}|n\rangle =(n-1){\hat {a}}|n\rangle .}

هذا يعني أنأ^{\displaystyle {\hat {a}}}أعمال بشأن|ن{\displaystyle |n\rangle }لإنتاج، حتى ثابت ضربي،|ن-1{\displaystyle |n-1\rangle }، وأ^{\displaystyle {\hat {a}}^{\dagger }}أعمال بشأن|ن{\displaystyle |n\rangle }لإنتاج|ن+1{\displaystyle |n+1\rangle }ولهذا السبب،أ^{\displaystyle {\hat {a}}}يُطلق عليه اسم عامل الإفناء ("عامل التخفيض")، وأ^{\displaystyle {\hat {a}}^{\dagger }}عامل الإنشاء (عامل الرفع). يُطلق على العاملين معًا اسم عوامل السلم .

بفرض أي حالة طاقة ذاتية، يمكننا تطبيق عامل التخفيض عليها، a ، لإنتاج حالة طاقة ذاتية أخرى بطاقة أقل بمقدار ħω . وبتكرار تطبيق عامل التخفيض، يبدو أنه يمكننا إنتاج حالات طاقة ذاتية تصل إلى E = −∞ . ومع ذلك، بما أن ن=ن|شمال|ن=ن|أأ|ن=(أ|ن)أ|ن0،{\displaystyle n=\langle n|N|n\rangle =\langle n|a^{\dagger }a|n\rangle ={\Bigl (}a|n\rangle {\Bigr )}^{\dagger }a|n\rangle \geqslant 0,}

أصغر قيمة ذاتية لمؤثر العدد هي صفر، و أ|0=0.{\displaystyle a\left|0\right\rangle =0.}

في هذه الحالة، ستؤدي التطبيقات اللاحقة لمؤثر التخفيض إلى إنتاج الصفر فقط، بدلاً من حالات طاقة ذاتية إضافية. علاوة على ذلك، فقد أوضحنا أعلاه أن ح^|0=ω2|0{\displaystyle {\hat {H}}\left|0\right\rangle ={\frac {\hbar \omega }{2}}\left|0\right\rangle }

وأخيراً، من خلال العمل على|0{\displaystyle |0\rangle }باستخدام عامل الرفع والضرب بمعاملات تطبيع مناسبة ، يمكننا إنتاج مجموعة لانهائية من حالات الطاقة الذاتية {|0،|1،|2،...،|ن،...}،{\displaystyle \left\{\left|0\right\rangle ,\left|1\right\rangle ,\left|2\right\rangle ,\ldots ,\left|n\right\rangle ,\ldots \right\},}

بحيث ح^|ن=ω(ن+12)|ن،{\displaystyle {\hat {H}}\left|n\right\rangle =\hbar \omega \left(n+{\frac {1}{2}}\right)\left|n\right\rangle ,} وهو ما يتطابق مع طيف الطاقة الوارد في القسم السابق.

يمكن التعبير عن الحالات الذاتية العشوائية بدلالة|0{\displaystyle |0\rangle }, [ 13 ]|ن=(أ)نن!|0.{\displaystyle |n\rangle ={\frac {(a^{\dagger })^{n}}{\sqrt {n!}}}|0\rangle .}

دليل

ن|أأ|ن=ن|([أ،أ]+أأ)|ن=ن|(شمال+1)|ن=ن+1أ|ن=ن+1|ن+1|ن=1نأ|ن-1=1ن(ن-1)(أ)2|ن-2==1ن!(أ)ن|0.{\displaystyle {\begin{aligned}\langle n|aa^{\dagger }|n\rangle &=\langle n|\left([a,a^{\dagger }]+a^{\dagger }a\right)\left|n\right\rangle =\langle n|\left(N+1\right)|n\rangle =n+1\\[1ex]\Rightarrow a^{\dagger }|n\rangle &={\sqrt {n+1}}|n+1\rangle \\[1ex]\Rightarrow |n\rangle &={\frac {1}{\sqrt {n}}}a^{\dagger }\left|n-1\right\rangle ={\frac {1}{\sqrt {n(n-1)}}}\left(a^{\dagger }\right)^{2}\left|n-2\right\rangle =\cdots ={\frac {1}{\sqrt {n!}}}\left(a^{\dagger }\right)^{n}\left|0\right\rangle .\end{aligned}}}

أسئلة تحليلية

التحليل السابق جبري، ويعتمد فقط على علاقات التبادل بين مؤثرات الرفع والخفض. بعد إتمام التحليل الجبري، ينبغي الانتقال إلى المسائل التحليلية. أولًا، يجب إيجاد الحالة الأرضية، أي حل المعادلةأψ0=0{\displaystyle a\psi _{0}=0}في تمثيل الموضع، هذه هي المعادلة التفاضلية من الدرجة الأولى (x+مωددx)ψ0=0،{\displaystyle \left(x+{\frac {\hbar }{m\omega }}{\frac {d}{dx}}\right)\psi _{0}=0,} والتي يمكن إيجاد حلها بسهولة على أنه التوزيع الغاوسي [ nb 1 ]ψ0(x)=جهـ-مωx22.{\displaystyle \psi _{0}(x)=Ce^{-{\frac {m\omega x^{2}}{2\hbar }}}.} من الناحية المفاهيمية، من المهم أن يكون هناك حل واحد فقط لهذه المعادلة؛ فلو كان هناك، على سبيل المثال، حالتان أرضيتان مستقلتان خطيًا، لحصلنا على سلسلتين مستقلتين من المتجهات الذاتية للمذبذب التوافقي. بمجرد حساب الحالة الأرضية، يمكن إثبات استقرائيًا أن الحالات المثارة هي متعددات حدود هيرميت مضروبة في الحالة الأرضية الغاوسية، باستخدام الصيغة الصريحة لمؤثر الرفع في تمثيل الموضع. كما يمكن إثبات أنه، كما هو متوقع من تفرد الحالة الأرضية، فإن دوال هيرميت هي حالات الطاقة الذاتية.ψن{\displaystyle \psi _{n}}تشكل الدوال التي تم إنشاؤها باستخدام طريقة السلم مجموعة كاملة من الدوال المتعامدة. [ 14 ]

بما أن دوال هيرميت إما زوجية أو فردية، فإنه يمكن إثبات أن متوسط ​​الإزاحة ومتوسط ​​الزخم يساوي صفرًا لجميع الحالات في QHO. [ 11 ]

بالربط الصريح مع القسم السابق، يتم تحديد الحالة الأرضية |0⟩ في تمثيل الموضع بواسطةأ|0=0{\displaystyle a|0\rangle =0}، x|أ|0=0(x+مωددx)x|0=0{\displaystyle \left\langle x\mid a\mid 0\right\rangle =0\qquad \Rightarrow \left(x+{\frac {\hbar }{m\omega }}{\frac {d}{dx}}\right)\left\langle x\mid 0\right\rangle =0\qquad \Rightarrow }x|0=(مωπ)14خبرة(-مω2x2)=ψ0 ،{\displaystyle \left\langle x\mid 0\right\rangle =\left({\frac {m\omega }{\pi \hbar }}\right)^{\frac {1}{4}}\exp \left(-{\frac {m\omega }{2\hbar }}x^{2}\right)=\psi _{0}~,} لذلك x|أ|0=ψ1(x) ،{\displaystyle \langle x\mid a^{\dagger }\mid 0\rangle =\psi _{1}(x)~,} لهذا السبب.ψ1(x،ت)=x|هـ-3أناωت/2أ|0{\displaystyle \psi _{1}(x,t)=\langle x\mid e^{-3i\omega t/2}a^{\dagger }\mid 0\rangle }وهكذا دواليك.

مقاييس الطول والطاقة الطبيعية

يمتلك المذبذب التوافقي الكمومي مقاييس طبيعية للطول والطاقة، والتي يمكن استخدامها لتبسيط المسألة. ويمكن إيجاد هذه المقاييس عن طريق تحويلها إلى صيغة لا بُعدية .

والنتيجة هي أنه إذا تم قياس الطاقة بوحدات ħω والمسافة بوحدات √ħ /() ، فإن الهاميلتوني يتبسط إلى ح=-12د2دx2+12x2،{\displaystyle H=-{\frac {1}{2}}{d^{2} \over dx^{2}}+{\frac {1}{2}}x^{2},} بينما تتبسط الدوال الذاتية للطاقة والقيم الذاتية إلى دوال هيرميت وأعداد صحيحة مزاحة بمقدار النصف، ψن(x)=x|ن=12نن! π-1/4خبرة(-x2/2) حن(x)،{\displaystyle \psi _{n}(x)=\left\langle x\mid n\right\rangle ={1 \over {\sqrt {2^{n}n!}}}~\pi ^{-1/4}\exp \left(-x^{2}/2\right)~H_{n}(x),}هـن=ن+12 ،{\displaystyle E_{n}=n+{\tfrac {1}{2}}~,} حيث أن H n ( x ) هي كثيرات حدود هيرميت .

لتجنب الالتباس، لن يتم اعتماد هذه "الوحدات الطبيعية" في هذا المقال في أغلب الأحيان. ومع ذلك، فهي مفيدة للغاية عند إجراء العمليات الحسابية، إذ تساعد على تجنب التعقيدات.

على سبيل المثال، فإن الحل الأساسي ( المولد ) لـ Hi∂ t ، وهو عامل شرودنغر المعتمد على الوقت لهذا المذبذب ، يختزل ببساطة إلى نواة ميهلر ، [ 15 ] [ 16 ]x|خبرة(-أناتح)|yك(x،y؛ت)=12πأناالخطيئةتخبرة(أنا2الخطيئةت((x2+y2)كوست-2xy)) ،{\displaystyle \left\langle x\right|\exp(-itH)\left|y\right\rangle \equiv K(x,y;t)={\frac {1}{\sqrt {2\pi i\sin t}}}\exp \left({\frac {i}{2\sin t}}\left(\left(x^{2}+y^{2}\right)\cos t-2xy\right)\right)~,} حيث K ( x , y ; 0) = δ ( xy ) . وبالتالي، فإن الحل الأكثر عمومية لتكوين ابتدائي مُعطى ψ ( x , 0) هو ببساطة ψ(x،ت)=دy ك(x،y؛ت)ψ(y،0).{\displaystyle \psi (x,t)=\int dy~K(x,y;t)\psi (y,0)\,.}

حالات متماسكة

ديناميكيات الحالة المتماسكة لـα=10{\displaystyle \alpha ={\sqrt {10}}}، بوحدات طول المذبذب التوافقيx0=/مω{\displaystyle x_{0}={\sqrt {\hbar /m\omega }}}، مما يوضح كثافة الاحتمال|ψ(x،ت)|2{\displaystyle |\psi (x,t)|^{2}}والطور الكمي (اللون).

تُعدّ الحالات المتماسكة (المعروفة أيضًا بحالات غلاوبر) للمذبذب التوافقي حزمًا موجية خاصة غير تشتتية ، ذات حد أدنى من عدم اليقين σx σp = 2 ، وتتطور قيمها المتوقعة للكميات القابلة للرصد كما في النظام الكلاسيكي. وهي متجهات ذاتية لمؤثر الإفناء، وليست للهاميلتوني، وتشكل أساسًا مكتملًا بشكل زائد، وبالتالي يفتقر إلى التعامد. [ 17 ]

يتم فهرسة الحالات المتماسكة بواسطةαج{\displaystyle \alpha \in \mathbb {C} }ويُعبَّر عنها في أساس | n كما يلي

|α=ن=0|نن|α=هـ-12|α|2ن=0αنن!|ن=هـ-12|α|2هـαأهـ-α*أ|0.{\displaystyle {\begin{aligned}|\alpha \rangle &=\sum _{n=0}^{\infty }|n\rangle \langle n|\alpha \rangle \\&=e^{-{\frac {1}{2}}|\alpha |^{2}}\sum _{n=0}^{\infty }{\frac {\alpha ^{n}}{\sqrt {n!}}}|n\rangle \\&=e^{-{\frac {1}{2}}|\alpha |^{2}}e^{\alpha a^{\dagger }}e^{-{\alpha ^{*}a}}|0\rangle .\end{aligned}}}

بما أن الحالات المتماسكة ليست حالات طاقة ذاتية، فإن تطورها الزمني ليس مجرد تحول بسيط في طور دالة الموجة. ومع ذلك، فإن الحالات المتطورة زمنيًا هي أيضًا حالات متماسكة، ولكن بمعامل إزاحة طور α بدلاً من ذلك.α(ت)=α(0)هـ-أناωت=α0هـ-أناωت{\displaystyle \alpha (t)=\alpha (0)e^{-i\omega t}=\alpha _{0}e^{-i\omega t}}.

|α(ت)=ن=0هـ-أنا(ن+12)ωت|نن|α=هـ-أناωت2هـ-12|α|2ن=0(αهـ-أناωت)نن!|ن=هـ-أناωت2|αهـ-أناωت{\displaystyle {\begin{aligned}|\alpha (t)\rangle &=\sum _{n=0}^{\infty }e^{-i\left(n+{\frac {1}{2}}\right)\omega t}|n\rangle \langle n|\alpha \rangle \\&=e^{\frac {-i\omega t}{2}}e^{-{\frac {1}{2}}|\alpha |^{2}}\sum _{n=0}^{\infty }{\frac {(\alpha e^{-i\omega t})^{n}}{\sqrt {n!}}}|n\rangle \\&=e^{-{\frac {i\omega t}{2}}}|\alpha e^{-i\omega t}\rangle \end{aligned}}}

لأنأ|0=0{\displaystyle a\left|0\right\rangle =0}ومن خلال متطابقة كيرماك-مكراي، فإن الشكل الأخير يعادل عامل إزاحة وحدوي يعمل على الحالة الأرضية:|α=هـαأ^-α*أ^|0=د(α)|0{\displaystyle |\alpha \rangle =e^{\alpha {\hat {a}}^{\dagger }-\alpha ^{*}{\hat {a}}}|0\rangle =D(\alpha )|0\rangle }حساب القيم المتوقعة:

x^α(ت)=2مω|α0|كوس(ωت-ϕ)ص^α(ت)=-2مω|α0|الخطيئة(ωت-ϕ){\displaystyle {\begin{aligned}\langle {\hat {x}}\rangle _{\alpha (t)}&={\sqrt {\frac {2\hbar }{m\omega }}}\left|\alpha _{0}\right|\cos {(\omega t-\phi )}\\[1ex]\langle {\hat {p}}\rangle _{\alpha (t)}&=-{\sqrt {2m\hbar \omega }}\left|\alpha _{0}\right|\sin {(\omega t-\phi )}\end{aligned}}}

أينϕ{\displaystyle \phi }تمثل هذه الطور المساهمة من المركب α . وتؤكد هذه المعادلات السلوك التذبذبي للجسيم.

الشكوك المحسوبة باستخدام الطريقة العددية هي:

σx(ت)=2مωσص(ت)=مω2{\displaystyle {\begin{aligned}\sigma _{x}(t)&={\sqrt {\frac {\hbar }{2m\omega }}}\\\sigma _{p}(t)&={\sqrt {\frac {m\hbar \omega }{2}}}\end{aligned}}}

مما يعطيσx(ت)σص(ت)=2{\textstyle \sigma _{x}(t)\sigma _{p}(t)={\frac {\hbar }{2}}}بما أن الدالة الموجية الوحيدة التي يمكن أن يكون لها أقل قدر من عدم اليقين في الموضع والزخم،2{\textstyle {\frac {\hbar }{2}}}، هي دالة موجية غاوسية، وبما أن دالة الموجة للحالة المتماسكة لها أقل قدر من عدم اليقين في الموضع والزخم، فإننا نلاحظ أن دالة الموجة الغاوسية العامة في ميكانيكا الكم لها الشكل التالي:ψα(x)=(مωπ)14هـأناص^α(x-x^α2)-مω2(x-x^α)2.{\displaystyle \psi _{\alpha }(x')=\left({\frac {m\omega }{\pi \hbar }}\right)^{\frac {1}{4}}e^{{\frac {i}{\hbar }}\langle {\hat {p}}\rangle _{\alpha }(x'-{\frac {\langle {\hat {x}}\rangle _{\alpha }}{2}})-{\frac {m\omega }{2\hbar }}(x'-\langle {\hat {x}}\rangle _{\alpha })^{2}}.}وباستبدال القيم المتوقعة كدالة للزمن، نحصل على دالة الموجة المتغيرة مع الزمن المطلوبة.

يمكن حساب احتمالية كل حالة من حالات الطاقة الذاتية لإيجاد توزيع الطاقة للدالة الموجية:

P(هـن)=|ن|α|2=هـ-|α|2|α|2نن!{\displaystyle P(E_{n})=\left|\langle n|\alpha \rangle \right|^{2}={\frac {e^{-|\alpha |^{2}}|\alpha |^{2n}}{n!}}}

وهو ما يتوافق مع توزيع بواسون .

حالات الإثارة الشديدة

دالة الموجة (أعلى) وكثافة الاحتمال (أسفل) للحالة المثارة n = 30 للمذبذب التوافقي الكمومي. تشير الخطوط الرأسية المتقطعة إلى نقاط التحول الكلاسيكية، بينما يمثل الخط المنقط كثافة الاحتمال الكلاسيكية.

عندما تكون قيمة n كبيرة، تتركز الحالات الذاتية في المنطقة المسموحة كلاسيكيًا، أي المنطقة التي يمكن لجسيم كلاسيكي ذي طاقة E<sub> n</sub> أن يتحرك فيها. وتكون الحالات الذاتية في ذروتها بالقرب من نقاط الانعطاف: وهي النقاط الواقعة على طرفي المنطقة المسموحة كلاسيكيًا حيث يغير الجسيم الكلاسيكي اتجاهه. ويمكن التحقق من هذه الظاهرة من خلال السلوك التقاربي لكثيرات حدود هيرميت ، وكذلك من خلال تقريب WKB .

يتناسب تردد التذبذب عند النقطة x طرديًا مع زخم p(x) لجسيم كلاسيكي طاقته En وموقعه x . علاوة على ذلك ، يتناسب مربع السعة (الذي يحدد كثافة الاحتمال) عكسيًا مع p ( x ) ، مما يعكس المدة الزمنية التي يقضيها الجسيم الكلاسيكي بالقرب من x . لا يوجد تفسير كلاسيكي بسيط لسلوك النظام في جوار صغير لنقطة التحول، ولكن يمكن نمذجته باستخدام دالة آيري . باستخدام خصائص دالة آيري، يمكن تقدير احتمال وجود الجسيم خارج المنطقة المسموح بها كلاسيكيًا، ليكون تقريبًا 2ن1/332/3Γ2(13)=1ن1/37.46408092658...{\displaystyle {\frac {2}{n^{1/3}3^{2/3}\Gamma ^{2}({\tfrac {1}{3}})}}={\frac {1}{n^{1/3}\cdot 7.46408092658...}}} وهذا يُعطى أيضًا، بشكل تقاربي، بواسطة التكامل 12π0هـ(2ن+1)(x-12سينه(2x))دx .{\displaystyle {\frac {1}{2\pi }}\int _{0}^{\infty }e^{(2n+1)\left(x-{\tfrac {1}{2}}\sinh(2x)\right)}dx~.}

حلول فضاء الطور

في صياغة فضاء الطور لميكانيكا الكم، يمكن كتابة الحالات الذاتية للمذبذب التوافقي الكمومي في عدة تمثيلات مختلفة لتوزيع شبه الاحتمال بصيغة مغلقة. وأكثر هذه التمثيلات استخداماً هو توزيع شبه الاحتمال لويغنر .

توزيع شبه الاحتمال لويغنر لحالة الطاقة الذاتية | n هو، بالوحدات الطبيعية الموضحة أعلاه، [ 18 ]Fن(x،ص)=(-1)نπلن(2(x2+ص2))هـ-(x2+ص2)،{\displaystyle F_{n}(x,p)={\frac {(-1)^{n}}{\pi \hbar }}L_{n}{\left(2(x^{2}+p^{2})\right)}e^{-(x^{2}+p^{2})}\,,} حيث تمثل L n كثيرات حدود لاغير . يوضح هذا المثال كيفية ارتباط كثيرات حدود هيرميت ولاغير من خلال خريطة ويغنر .

في الوقت نفسه، تتخذ دالة هوسيمي Q لحالات الطاقة الذاتية للمذبذب التوافقي شكلاً أبسط. إذا عملنا بالوحدات الطبيعية الموضحة أعلاه، فسنحصل على سؤالن(x،ص)=(x2+ص2)نن!هـ-(x2+ص2)π{\displaystyle Q_{n}(x,p)={\frac {(x^{2}+p^{2})^{n}}{n!}}{\frac {e^{-(x^{2}+p^{2})}}{\pi }}} يمكن التحقق من هذا الادعاء باستخدام تحويل سيغال-بارغمان . تحديدًا، بما أن عامل الرفع في تمثيل سيغال-بارغمان هو ببساطة الضرب فيz=x+أناص{\displaystyle z=x+ip}والحالة الأرضية هي الدالة الثابتة 1، وحالات المذبذب التوافقي المعياري في هذا التمثيل هي ببساطةzن/ن!{\displaystyle z^{n}/{\sqrt {n!}}}عند هذه النقطة، يمكننا اللجوء إلى صيغة دالة Husimi Q من حيث تحويل Segal–Bargmann.

المذبذبات التوافقية ثنائية الأبعاد

تمت معالجة المذبذب التوافقي ثنائي الأبعاد في الإحداثيات الديكارتية والمذبذب التوافقي المتساوي الخواص ثنائي الأبعاد في الإحداثيات الأسطوانية بالتفصيل في كتاب مولر-كيرستن. [ 19 ]

مذبذب توافقي متساوي الخواص ذو أبعاد N

يمكن تعميم المذبذب التوافقي أحادي البعد بسهولة إلى N بُعدًا، حيث N = 1، 2، 3، ... . في بُعد واحد، يُحدد موضع الجسيم بإحداثي واحد ، x . في N بُعدًا، يُستبدل هذا الإحداثي بـ N إحداثيات موضع، نُرمز لها بـ x₁ ، ...، xₙ . لكل إحداثي موضع زخم ، نُرمز له بـ p₁ ، ... ، pₙ . علاقات التبادل الأساسية بين هذه المؤثرات هي [xأنا،صج]=أنادلتاأنا،ج[xأنا،xج]=0[صأنا،صج]=0{\displaystyle {\begin{aligned}{[}x_{i},p_{j}{]}&=i\hbar \delta _{i,j}\\{[}x_{i},x_{j}{]}&=0\\{[}p_{i},p_{j}{]}&=0\end{aligned}}}

الهاميلتوني لهذا النظام هو ح=أنا=1شمال(صأنا22م+12مω2xأنا2).{\displaystyle H=\sum _{i=1}^{N}\left({p_{i}^{2} \over 2m}+{1 \over 2}m\omega ^{2}x_{i}^{2}\right).}

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

هذه الملاحظة تجعل الحل واضحاً. بالنسبة لمجموعة معينة من الأعداد الكمومية{ن}{ن1،ن2،...،نشمال}{\displaystyle \{n\}\equiv \{n_{1},n_{2},\dots ,n_{N}\}}يتم التعبير عن الدوال الذاتية للطاقة للمذبذب ذي الأبعاد N بدلالة الدوال الذاتية أحادية البعد على النحو التالي:

x|ψ{ن}=أنا=1شمالxأنا|ψنأنا{\displaystyle \langle \mathbf {x} |\psi _{\{n\}}\rangle =\prod _{i=1}^{N}\langle x_{i}\mid \psi _{n_{i}}\rangle }

في طريقة عامل السلم، نقوم بتعريف N مجموعة من عوامل السلم،

أأنا=مω2(xأنا+أنامωصأنا)،أأنا=مω2(xأنا-أنامωصأنا).{\displaystyle {\begin{aligned}a_{i}&={\sqrt {m\omega \over 2\hbar }}\left(x_{i}+{i \over m\omega }p_{i}\right),\\a_{i}^{\dagger }&={\sqrt {m\omega \over 2\hbar }}\left(x_{i}-{i \over m\omega }p_{i}\right).\end{aligned}}}

وباتباع إجراء مماثل للحالة أحادية البعد، يمكننا أن نُبين أن كلًا من المؤثرين a i و a i يُخفض الطاقة ويرفعها بمقدار ℏω على التوالي. الهاميلتوني هو ح=ωأنا=1شمال(أأناأأنا+12).{\displaystyle H=\hbar \omega \,\sum _{i=1}^{N}\left(a_{i}^{\dagger }\,a_{i}+{\frac {1}{2}}\right).} هذا الهاميلتوني ثابت تحت مجموعة التناظر الديناميكي U ( N ) (المجموعة الوحدوية في N بُعدًا)، المعرفة بواسطة يوأأنايو=ج=1شمالأجيوجأناللجميعيويو(شمال)،{\displaystyle U\,a_{i}^{\dagger }\,U^{\dagger }=\sum _{j=1}^{N}a_{j}^{\dagger }\,U_{ji}\quad {\text{for all}}\quad U\in U(N),} أينيوجأنا{\displaystyle U_{ji}}هو عنصر في التمثيل المصفوفي المحدد لـ U ( N ) .

مستويات الطاقة في النظام هي هـ=ω[(ن1++نشمال)+شمال2].{\displaystyle E=\hbar \omega \left[(n_{1}+\cdots +n_{N})+{N \over 2}\right].}نأنا=0،1،2،...(مستوى الطاقة في البعد أنا).{\displaystyle n_{i}=0,1,2,\dots \quad ({\text{the energy level in dimension }}i).}

كما هو الحال في الحالة أحادية البعد، تكون الطاقة مُكمّمة. طاقة الحالة الأرضية تساوي N ضعف طاقة الحالة الأرضية أحادية البعد، كما هو متوقع عند استخدام القياس على N من المذبذبات أحادية البعد المستقلة. ثمة فرق آخر: في الحالة أحادية البعد، يتوافق كل مستوى طاقة مع حالة كمومية فريدة. أما في N بُعد، فباستثناء الحالة الأرضية، تكون مستويات الطاقة متساوية الطاقة ، أي أن هناك عدة حالات لها نفس الطاقة.

يمكن حساب درجة الانحلال بسهولة نسبية. على سبيل المثال، لنفترض الحالة ثلاثية الأبعاد: لنُعرّف n = n₁ + n₂ + n₃ . جميع الحالات التي لها نفس قيمة n ستكون لها نفس الطاقة. بالنسبة لقيمة n معينة ، نختار قيمة n₁ محددة . عندئذٍ ، n₂ + n₃ = n - n₁ . يوجد n - n₁ + 1 زوجًا ممكنًا { n₂ , n₃ } . يمكن أن تأخذ n₂ القيم من 0 إلى n - n₁ ، ولكل قيمة n₂ ، تكون قيمة n₃ ثابتة. وبالتالي ، تكون درجة الانحلال كما يلي :زن=ن1=0نن-ن1+1=(ن+1)(ن+2)2{\displaystyle g_{n}=\sum _{n_{1}=0}^{n}n-n_{1}+1={\frac {(n+1)(n+2)}{2}}} صيغة عامة لـ N و n [ حيث n هو بُعد التمثيل المتناظر غير القابل للاختزال من الرتبة n للمجموعة الوحدوية U ( N ) ]: زن=(شمال+ن-1ن)=(شمال+ن-1)!ن!(شمال-1)!،ن=ن1+ن2+ن3+.{\displaystyle g_{n}={\binom {N+n-1}{n}}={\frac {(N+n-1)!}{n!(N-1)!}},\;\;\;n=n_{1}+n_{2}+n_{3}+\cdots .} تنتج الحالة الخاصة N = 3، المذكورة أعلاه، مباشرةً من هذه المعادلة العامة. إلا أن هذا لا يصح إلا للجسيمات القابلة للتمييز، أي في إحصاءات ماكسويل-بولتزمان (وليس في الإحصاءات الكمومية) أو جسيم واحد في N بُعدًا (حيث تكون الأبعاد قابلة للتمييز). في حالة N بوزون في مصيدة توافقية أحادية البُعد، يتناسب الانحلال طرديًا مع عدد طرق تقسيم عدد صحيح n باستخدام أعداد صحيحة أقل من أو تساوي N. ويمكن إثبات أن العدد الكبير-هـ{\displaystyle E}السلوك التقاربي للانحلالزن{\displaystyle g_{n}}وهو مستقل عمليًا عن الطاقةهـ{\displaystyle E}- يختلف عن الحالة الكلاسيكية التي ينحرف فيها هذا الأمر. [ 20 ] هذا الانحطاط هو

زن=ص(شمال-،ن).{\displaystyle g_{n}=p(N_{-},n).}

ينشأ هذا بسبب قيد وضع N من الكمات في حالة كيت حيثك=0كنك=ن{\textstyle \sum _{k=0}^{\infty }kn_{k}=n} وك=0نك=شمال{\textstyle \sum _{k=0}^{\infty }n_{k}=N}وهي نفس القيود الموجودة في تقسيم الأعداد الصحيحة.

مثال: مذبذب توافقي متناحٍ ثلاثي الأبعاد

حلول مدارية توافقية كروية ثلاثية الأبعاد لشرودنغر في مخططات كثافة ثنائية الأبعاد؛ يوجد كود المصدر لبرنامج Mathematica المستخدم لإنشاء المخططات في الأعلى

يمكن حل معادلة شرودنغر لجسيم في مذبذب توافقي ثلاثي الأبعاد متناظر كرويًا بشكل صريح عن طريق فصل المتغيرات. هذه العملية مماثلة للفصل الذي تم إجراؤه في مسألة الذرة الشبيهة بالهيدروجين ، ولكن مع جهد متناظر كرويًا مختلف.V(ر)=12μω2ر2،{\displaystyle V(r)={1 \over 2}\mu \omega ^{2}r^{2},} حيث μ هي كتلة الجسيم. ولأن m ستُستخدم لاحقًا للدلالة على العدد الكمي المغناطيسي، فسيتم الإشارة إلى الكتلة بـ μ بدلاً من m كما ذُكر سابقًا في هذه المقالة.

حل المعادلة هو: [ 21 ]ψكلم(ر،θ،ϕ)=شمالكلرلهـ-νر2لك(ل+12)(2νر2)Yلم(θ،ϕ){\displaystyle \psi _{klm}(r,\theta ,\phi )=N_{kl}r^{l}e^{-\nu r^{2}}L_{k}^{\left(l+{1 \over 2}\right)}(2\nu r^{2})Y_{lm}(\theta ,\phi )} أين

شمالكل=2ν3π2ك+2ل+3ك!νل(2ك+2ل+1)!!  {\displaystyle N_{kl}={\sqrt {{\sqrt {\frac {2\nu ^{3}}{\pi }}}{\frac {2^{k+2l+3}\;k!\;\nu ^{l}}{(2k+2l+1)!!}}}}~~}هو ثابت تطبيع؛νμω2 {\displaystyle \nu \equiv {\mu \omega \over 2\hbar }~}؛
لك(ل+12)(2νر2){\displaystyle {L_{k}}^{(l+{1 \over 2})}(2\nu r^{2})}

هي كثيرات حدود لاغير المعممة ؛ رتبة k لكثير الحدود هي عدد صحيح غير سالب، يتطابق مع عدد عقد الجزء القطري من الدالة الموجية؛

قيمة الطاقة الذاتية هي هـ=ω(2ك+ل+32).{\displaystyle E=\hbar \omega \left(2k+l+{\frac {3}{2}}\right).} عادةً ما يتم وصف الطاقة بواسطة عدد كمي واحدن2ك+ل.{\displaystyle n\equiv 2k+l\,.}

بما أن k عدد صحيح غير سالب، فلكل عدد زوجي n لدينا = 0، 2، ...، n − 2، ولكل عدد فردي n لدينا = 1، 3، ...، n − 2، n . العدد الكمي المغناطيسي m عدد صحيح يحقق m ، لذا لكل n و يوجد 2   + 1 حالة كمية  مختلفة ، مُصنَّفة بـ m . وبالتالي، فإن الانحلال عند المستوى n هو ل=...،ن-2،ن(2ل+1)=(ن+1)(ن+2)2،{\displaystyle \sum _{l=\ldots ,n-2,n}(2l+1)={(n+1)(n+2) \over 2}\,,} حيث يبدأ المجموع من 0 أو 1، وفقًا لما إذا كان n زوجيًا أم فرديًا. تتوافق هذه النتيجة مع صيغة الأبعاد المذكورة أعلاه، وتُعادل بُعد التمثيل المتناظر لـ SU(3) ، [ 22 ] وهي مجموعة الانحلال ذات الصلة.

التطبيقات

شبكة المذبذبات التوافقية: الفونونات

يمكن توسيع مفهوم المذبذب التوافقي ليشمل شبكة أحادية البعد تتكون من العديد من الجسيمات. لنفترض سلسلة توافقية كمومية أحادية البعد مكونة من N ذرة متطابقة. هذا هو أبسط نموذج كمومي للشبكة، وسنرى كيف تنشأ الفونونات منه. يمكن تعميم الصيغة التي سنطورها لهذا النموذج بسهولة إلى بعدين وثلاثة أبعاد. كما في القسم السابق، نرمز إلى مواضع الكتل بـ x₁ , x₂ , ... ، مقاسة من مواضع اتزانها (أي xᵢ = 0 إذا كان الجسيم i في موضع اتزانه). في بعدين أو أكثر، تكون xᵢ كميات متجهة. هاميلتوني هذا النظام هو

ح=أنا=1شمالصأنا22م+12مω2{أناج}(نن)(xأنا-xج)2،{\displaystyle \mathbf {H} =\sum _{i=1}^{N}{p_{i}^{2} \over 2m}+{1 \over 2}m\omega ^{2}\sum _{\{ij\}(nn)}(x_{i}-x_{j})^{2}\,,} حيث m هي كتلة كل ذرة (بافتراض أنها متجانسة)، و xᵢ و pᵢ هما مؤثرا الموضع والزخم للذرة رقم i ، ويتم الجمع على أقرب الجيران (nn). مع ذلك، جرت العادة على إعادة كتابة الهاميلتوني بدلالة الأنماط الطبيعية لمتجه الموجة بدلاً من إحداثيات الجسيم، وذلك لتسهيل العمل في فضاء فورييه الأكثر ملاءمة .

تراكب ثلاثة ثنائيات أقطاب متذبذبة - يوضح الانتشار الزمني للدالة الموجية المشتركة لقيم مختلفة لـ n و l و m

نقدم إذن مجموعة من N "الإحداثيات الطبيعية" Q k ، المعرفة على أنها تحويلات فورييه المنفصلة لـ x s ، و N "الزخم المترافق" Π المعرفة على أنها تحويلات فورييه لـ p s ، سؤالك=1شماللهـأناكألxل{\displaystyle Q_{k}={1 \over {\sqrt {N}}}\sum _{l}e^{ikal}x_{l}}Πك=1شماللهـ-أناكألصل.{\displaystyle \Pi _{k}={1 \over {\sqrt {N}}}\sum _{l}e^{-ikal}p_{l}\,.}

ستكون الكمية k n هي العدد الموجي للفونون، أي مقسومًا على الطول الموجي . وهي تأخذ قيمًا كمية، لأن عدد الذرات محدود.

يحافظ هذا على علاقات التبادل المطلوبة إما في الفضاء الحقيقي أو في فضاء متجه الموجة

يوضح مثال آخر لانتشار الدالة الموجية المشتركة مع الزمن لثلاث ذرات مختلفة تأثير الزخم الزاوي على سلوك التوزيع

[xل،صم]=أنادلتال،م[سؤالك،Πك]=1شمالل،مهـأناكألهـ-أناكأم[xل،صم]=أناشمالمهـأناأم(ك-ك)=أنادلتاك،ك[سؤالك،سؤالك]=[Πك،Πك]=0 .{\displaystyle {\begin{aligned}\left[x_{l},p_{m}\right]&=i\hbar \delta _{l,m}\\\left[Q_{k},\Pi _{k'}\right]&={1 \over N}\sum _{l,m}e^{ikal}e^{-ik'am}[x_{l},p_{m}]\\&={i\hbar \over N}\sum _{m}e^{iam(k-k')}=i\hbar \delta _{k,k'}\\\left[Q_{k},Q_{k'}\right]&=\left[\Pi _{k},\Pi _{k'}\right]=0~.\end{aligned}}}

من النتيجة العامة لxلxل+م=1شمالككسؤالكسؤالكلهـأناأل(ك+ك)هـأناأمك=كسؤالكسؤال-كهـأناأمكلصل2=كΠكΠ-ك ،{\displaystyle {\begin{aligned}\sum _{l}x_{l}x_{l+m}&={1 \over N}\sum _{kk'}Q_{k}Q_{k'}\sum _{l}e^{ial\left(k+k'\right)}e^{iamk'}=\sum _{k}Q_{k}Q_{-k}e^{iamk}\\\sum _{l}{p_{l}}^{2}&=\sum _{k}\Pi _{k}\Pi _{-k}~,\end{aligned}}} من السهل إثبات، من خلال حساب المثلثات الأساسي، أن حد طاقة الوضع هو 12مω2ج(xج-xج+1)2=12مω2كسؤالكسؤال-ك(2-هـأناكأ-هـ-أناكأ)=12مكωك2سؤالكسؤال-ك ،{\displaystyle {\begin{aligned}{1 \over 2}m\omega ^{2}\sum _{j}\left(x_{j}-x_{j+1}\right)^{2}&={1 \over 2}m\omega ^{2}\sum _{k}Q_{k}Q_{-k}\left(2-e^{ika}-e^{-ika}\right)\\&={1 \over 2}m\sum _{k}{\omega _{k}}^{2}Q_{k}Q_{-k}~,\end{aligned}}} أين ωك=2ω2(1-كوس(كأ)) .{\displaystyle \omega _{k}={\sqrt {2\omega ^{2}(1-\cos(ka))}}~.}

يمكن كتابة الهاميلتوني في فضاء متجه الموجة على النحو التالي: ح=12مك(ΠكΠ-ك+م2ωك2سؤالكسؤال-ك) .{\displaystyle \mathbf {H} ={1 \over {2m}}\sum _{k}\left({\Pi _{k}\Pi _{-k}}+m^{2}\omega _{k}^{2}Q_{k}Q_{-k}\right)~.}

لاحظ أن اقترانات متغيرات الموضع قد تم تحويلها بعيدًا؛ إذا كانت Q s و Π s هيرميتية (وهي ليست كذلك)، فإن الهاميلتوني المحول سيصف N من المذبذبات التوافقية غير المقترنة .

يعتمد شكل التكميم على اختيار الشروط الحدية؛ ولتبسيط الأمر، نفرض شروطًا حدية دورية ، مع تعريف الذرة ( N + 1) على أنها مكافئة للذرة الأولى. فيزيائيًا، يتوافق هذا مع وصل السلسلة من طرفيها. والتكميم الناتج هو

ك=كن=2نπشمالأل ن=0،±1،±2،...،±شمال2.{\displaystyle k=k_{n}={2n\pi \over Na}\quad {\hbox{for}}\ n=0,\pm 1,\pm 2,\ldots ,\pm {N \over 2}.}

يأتي الحد الأعلى لـ n من الحد الأدنى للطول الموجي، وهو ضعف المسافة بين الشبكات a ، كما نوقش أعلاه.

القيم الذاتية للمذبذب التوافقي أو مستويات الطاقة للنمط ω k هي هـن=(12+ن)ωكلن=0،1،2،3،...{\displaystyle E_{n}=\left({1 \over 2}+n\right)\hbar \omega _{k}\quad {\hbox{for}}\quad n=0,1,2,3,\ldots }

إذا تجاهلنا طاقة النقطة الصفرية ، فإن المستويات تكون متباعدة بانتظام عند 0، ω، 2ω، 3ω، {\displaystyle 0,\ \hbar \omega ,\ 2\hbar \omega ,\ 3\hbar \omega ,\ \cdots }

لذا ، يجب تزويد شبكة المذبذب التوافقي بكمية محددة من الطاقة ħω لدفعها إلى مستوى الطاقة التالي. قياسًا على حالة الفوتون عندما يكون المجال الكهرومغناطيسي مُكمّمًا، يُطلق على كمّ الطاقة الاهتزازية اسم فونون .

تُظهر جميع الأنظمة الكمومية خصائص موجية وخصائص جسيمية. ويمكن فهم الخصائص الجسيمية للفونون على أفضل وجه باستخدام طرق التكميم الثاني وتقنيات المؤثرات الموصوفة في مكان آخر. [ 23 ]

في حالة التقارب المستمر ، a → 0 و N → ∞ ، بينما تبقى Na ثابتة. تتحول الإحداثيات الكنسية Q k إلى أنماط الزخم المنفصلة لحقل قياسي.ϕك{\displaystyle \phi _{k}}بينما يصبح مؤشر الموقع i ( وليس متغير الإزاحة الديناميكي ) هو وسيط المعلمة x للحقل القياسي،ϕ(x،ت){\displaystyle \phi (x,t)}.

الاهتزازات الجزيئية

  • تُعدّ اهتزازات الجزيء ثنائي الذرة مثالاً على نموذج ثنائي الجسم للمذبذب التوافقي الكمومي. في هذه الحالة، يُعطى التردد الزاوي بالعلاقة التالية:ω=كμ{\displaystyle \omega ={\sqrt {\frac {k}{\mu }}}}أينμ=م1م2م1+م2{\displaystyle \mu ={\frac {m_{1}m_{2}}{m_{1}+m_{2}}}}الكتلة المختزلة وم1{\displaystyle m_{1}}وم2{\displaystyle m_{2}}[ 24 ] هما كتل الذرتين.
  • نمذجة الفونونات، كما نوقش أعلاه.
  • رسومq{\displaystyle q}مع الكتلةم{\displaystyle m}في مجال مغناطيسي منتظمب{\displaystyle \mathbf {B} }يُعد مثالاً على المذبذب التوافقي الكمي أحادي البعد: تكميم لاندو .
  • يُقارب نموذج المذبذب التوافقي الجهد بين النوى لجزيء ثنائي الذرات، حيث تشبه الحالات الاهتزازية المنخفضة النموذج بشكل كبير، بينما تنحرف الحالات الأعلى بسبب اللا توافقية. [ 11 ]

قانون هوك

  • ذرة هوك هي نموذج بسيط لذرة الهيليوم باستخدام المذبذب التوافقي الكمي.
  • يُمثل قانون هوك نموذجاً لكتلة تتحرك على نابض حيث تتناسب القوة المؤثرة على الكتلة طردياً مع إزاحتها. [ 11 ]
  • يمكن اشتقاق الحل العام لكتلة على نابض من هذا الافتراض. [ 11 ]
  • تصل إزاحة الكتلة إلى أقصى حد وأدنى حد عند A و -A حيث تسمى A السعة.
  • يُعتبر هذا النظام نظامًا محافظًا حيث تظل الطاقة الكلية دون تغيير ويتم إعادة توزيعها باستمرار بين الطاقة الحركية والطاقة الكامنة.

المذبذب التوافقي المعكوس

تمت دراسة المذبذب التوافقي المعكوس بالتفصيل من قبل جي. بارتون. [ 25 ] انظر أيضًا إتش جيه دبليو مولر-كيرستن [ 26 ] وسي. يوس، إيه. كيلين وإيه. كوروه. [ 27 ]

مذبذب ديراك

دراسة المذبذب التوافقي، على سبيل المثال، من منظور الطاقةهـ=ص2/2م+مω2q2/2{\displaystyle E={\mathbf {p}}^{2}/2m+m\omega ^{2}q^{2}/2}قياسًا على اشتقاق معادلة ديراك - إن صح التعبير من "الجذر التربيعي" للمعادلةصμصμ+م2=0{\displaystyle p_{\mu }p^{\mu }+m^{2}=0}– تم استكشافها بواسطة لوريلا م. جونز. [ 28 ]

انظر أيضاً

ملحوظات

  1. ثابت التطبيع هوج=(مωπ)1/4{\displaystyle C=\left({\frac {m\omega }{\pi \hbar }}\right)^{{1}/{4}}}، ويحقق شرط التطبيع-ψ0(x)*ψ0(x)دx=1{\displaystyle \int _{-\infty }^{\infty }\psi _{0}(x)^{*}\psi _{0}(x)dx=1}.

مراجع

  1. غريفيثس 2004 .
  2. ليبوف 2002 .
  3. رشيد، منير أ. (2006). "سعة الانتقال للمذبذب التوافقي الخطي المتغير مع الزمن مع إضافة حدود خطية متغيرة مع الزمن إلى الهاميلتوني" (ملف PDF) . مركز محمد علي رشيد للرياضيات والفيزياء المتقدمة . المركز الوطني للفيزياء . مؤرشف من الأصل ( ملف PDF - مايكروسوفت باوربوينت ) في 3 مارس 2016. تم الاطلاع عليه في 19 أكتوبر 2010 .
  4. هارالد جيه دبليو مولر-كيرستن، مقدمة في ميكانيكا الكم: معادلة شرودنغر وتكامل المسار، الطبعة الثانية (2012)، دار النشر العالمية، ISBN 978-9810-4397-5. يتضمن هذا المرجع شرحًا مفصلًا لطريقة المؤثرات، ودوال هيرميت، والتكامل الكفافي، الصفحات 130-141.
  5. ^ زويباخ (2022) ، ص 233-234.
  6. Zwiebach (2022) ، ص 234.
  7. Zwiebach (2022) ، ص 241.
  8. غبور، غريغوري ج. (2011). الأساليب الرياضية للفيزياء البصرية والهندسة . مطبعة جامعة كامبريدج. ص 631-633 . ISBN  978-0-521-51610-5.
  9. Zwiebach (2022) ، ص 240.
  10. ^ زويباخ (2022) ، ص 249-250.
  11. 1 2 3 4 5 ماكواري، دونالد (7 أغسطس 2023). الكيمياء الفيزيائية: مدخل جزيئي ( الطبعة الأولى). منشورات جامعة العلوم. الصفحات 157-173 . ISBN   9781891389962.
  12. ^ زويباخ (2022) ، ص 246-249.
  13. Zwiebach (2022) ، ص 248.
  14. هول، برايان سي. (2013)، نظرية الكم للرياضيين ، نصوص الدراسات العليا في الرياضيات، المجلد 267، سبرينغر، النظرية 11.4، رمز Bibcode : 2013qtm..book.....H ، ISBN  978-1461471158
  15. باولي، و. (2000)، ميكانيكا الموجات: المجلد 5 من محاضرات باولي في الفيزياء (كتب دوفر في الفيزياء). ISBN 978-0486414621 القسم 44.
  16. كوندون، الاتحاد الأوروبي (1937). "غمر تحويل فورييه في مجموعة متصلة من التحويلات الوظيفية"، وقائع الأكاديمية الوطنية للعلوم في الولايات المتحدة الأمريكية 23 ، 158-164. متاح على الإنترنت
  17. ^ زويباخ (2022) ، ص 481-492.
  18. سيريسولا، فيديريكو؛ مايو، فرانكو؛ رونكاليا، أوغوستو ج. (11 أكتوبر 2023). "توزيع ويغنر شبه الاحتمالي للعمل" . إنتروبي . 25 ( 10): 1439. arXiv : 2303.08755 . Bibcode : 2023Entrp..25.1439C . doi : 10.3390/e25101439 . ISSN 1099-4300 . PMC 10606729. PMID 37895560 .   
  19. هارالد جيه دبليو مولر-كيرستن، مقدمة في ميكانيكا الكم: معادلة شرودنغر وتكامل المسار، الطبعة الثانية (2012)، وورلد ساينتيفيك، ص 111-114.
  20. هارالد جيه دبليو مولر-كيرستن، أساسيات الفيزياء الإحصائية، الطبعة الثالثة (2022)، وورلد ساينتيفيك، ISBN 978-9810-125-609-7، الصفحات 203-204.
  21. ألبرت ميسيا ، ميكانيكا الكم ، 1967، نورث هولاند، الفصل الثاني عشر، القسم 15، صفحة 456. متوفر على الإنترنت
  22. فرادكين، د.م. (1965). "المذبذب التوافقي المتناحي ثلاثي الأبعاد وSU3". المجلة الأمريكية للفيزياء . 33 (3): 207-211 . doi : 10.1119/1.1971373 .
  23. ماهان، جي دي (1981). فيزياء الجسيمات المتعددة . نيويورك: سبرينغر. ISBN 978-0306463389.
  24. "المذبذب التوافقي الكمي" . هايبرفيزيكس . تم الاطلاع عليه بتاريخ 24 سبتمبر 2009 .
  25. جي. بارتون. المذبذب التوافقي المعكوس، حوليات الفيزياء (نيويورك) 166 (1986) 322.
  26. هارالد ج. مولر-كيرستن، مقدمة في ميكانيكا الكم: معادلة شرودنغر وتكامل المسار، الطبعة الثانية (2012) وورلد ساينتيفيك، ص 130-141.
  27. C. Yuce, A. Killen and A. Coruh, Invert Oscillator, quant-ph/0703234.
  28. لوريلا م. جونز، مذبذب ديراك آخر، تقرير جامعة إلينوي في أوربانا-شامبين ILLL-(TH)-91-24 (1991).

فهرس