ميكانيكا الكم العشوائية

ميكانيكا الكم العشوائية هي إطار لوصف ديناميكيات الجسيمات الخاضعة لعمليات عشوائية داخلية، بالإضافة إلى قوى خارجية متنوعة. يوفر هذا الإطار اشتقاقًا لمعادلات الانتشار المرتبطة بهذه الجسيمات العشوائية. وتشتهر ميكانيكا الكم العشوائية باشتقاقها لمعادلة شرودنغر، المعروفة أيضًا بمعادلة كولموغوروف، لنوع معين من الانتشار المحافظ (أو الوحدوي). [ 1 ] [ 2 ]

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

في ثلاثينيات القرن العشرين، أدرك كلٌ من إرفين شرودنغر وراينهولد فورث وجود تشابه بين معادلات الانتشار الكلاسيكي وصيغة نظرية الكم، إلا أن أول نظرية عشوائية متماسكة نسبيًا لميكانيكا الكم طُرحت عام ١٩٤٦ على يد الفيزيائي المجري إيمري فينيس . [ ٤ ] [ ٥ ] [ ٦ ] [ ٧ ] شعر لويس دي بروي [ ٨ ] بضرورة دمج عملية عشوائية كامنة في ميكانيكا الكم لجعل الجسيمات تنتقل من موجة مرجعية إلى أخرى. [ ٧ ] تُنسب نظرية ميكانيكا الكم العشوائية إلى إدوارد نيلسون ، الذي اكتشف بشكل مستقل اشتقاقًا لمعادلة شرودنغر ضمن هذا الإطار. [ ١ ] [ ٢ ] كما طوّر هذه النظرية كلٌ من ديفيدسون، وغيرا ، وروجييرو، وبافون، وآخرون. [ ٧ ]

التفسير العشوائي لميكانيكا الكم

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

هذا التفسير معروف جيدًا في سياق الميكانيكا الإحصائية ، [ 9 ] وحركة براون على وجه الخصوص. لذا، وفقًا للتفسير العشوائي، ينبغي تفسير ميكانيكا الكم بطريقة مشابهة لحركة براون. [ 1 ] مع ذلك، في حالة حركة براون، فإن وجود مقياس احتمالي (يُسمى مقياس وينر [ 10 ] ) يُحدد التكامل المساري الإحصائي أمرٌ راسخ، ويمكن توليد هذا المقياس بواسطة عملية عشوائية تُسمى عملية وينر . [ 11 ] من ناحية أخرى، يواجه إثبات وجود مقياس احتمالي يُحدد التكامل المساري في ميكانيكا الكم صعوبات، [ 12 ] [ 13 ] وليس من المؤكد إمكانية توليد مثل هذا المقياس الاحتمالي بواسطة عملية عشوائية. الميكانيكا العشوائية هي الإطار المعني ببناء مثل هذه العمليات العشوائية التي تُولد مقياسًا احتماليًا لميكانيكا الكم.

في الحركة البراونية ، من المعروف أن التقلبات الإحصائية للجسيم البراوني غالبًا ما تنتج عن تفاعله مع عدد كبير من الجسيمات المجهرية. [ 14 ] [ 15 ] [ 16 ] [ 17 ] في هذه الحالة، يُستخدم وصف الحركة البراونية باستخدام عملية وينر كتقريب فقط، متجاهلاً ديناميكيات الجسيمات الفردية في الخلفية. وبدلاً من ذلك، يصف تأثير هذه الجسيمات الخلفية من خلال سلوكها الإحصائي.

لا يُحدد التفسير العشوائي لميكانيكا الكم أصل التقلبات الكمومية للجسيم الكمومي. فهو يُقدم هذه التقلبات كنتيجة لقانون عشوائي جديد في الطبيعة يُسمى فرضية الخلفية. [ 2 ] يمكن تفسير هذه الفرضية على أنها تطبيق دقيق لمقولة "الله يلعب النرد"، لكنها تُبقي الباب مفتوحًا أمام إمكانية استبدال هذه اللعبة بنظرية المتغيرات الخفية، كما هو الحال في نظرية الحركة البراونية. [ 18 ]

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

التكميم العشوائي

يمكن تلخيص مسلمات النظرية في شرط التكميم العشوائي الذي صاغه نيلسون. [ 2 ] بالنسبة لنظرية غير نسبية حولRن{\displaystyle \mathbb {R} ^{n}}ينص هذا الشرط على ما يلي: [ 19 ]

  • يتم وصف مسار الجسيم الكمومي بواسطة الإسقاط الحقيقي لشبه مارتينجال معقد :X(ت)=Rهـ[Z(ت)]{\displaystyle X(t)={\rm {Re}}[Z(t)]}معZ(ت)=ج(ت)+م(ت){\displaystyle Z(t)=C(t)+M(t)}، أينج(ت){\displaystyle C(t)}هي عملية تغير محدود مستمر وم(ت){\displaystyle M(t)}هو مارتينجال معقد ،
  • يُحسّن المسار بشكل عشوائي الفعلS=هـ[لدت]{\displaystyle S=\mathbb {E} \left[\int L\,dt\right]}،
  • المارتينجالم(ت){\displaystyle M(t)}هي عملية مستمرة ذات زيادات مستقلة وعزوم محدودة . علاوة على ذلك، فإن تغيرها التربيعي محدد بعلاقة البنية.د[مأنا،مج]=αمدلتاأناجدت{\displaystyle d[M^{i},M^{j}]={\frac {\alpha \,\hbar }{m}}\,\delta ^{ij}\,dt}أينم{\displaystyle m}كتلة الجسيم،{\displaystyle \hbar }ثابت بلانك المختزل ،α=|α|هـأناϕج{\displaystyle \alpha =|\alpha |e^{{\rm {i}}\phi }\in \mathbb {C} }ثابت بلا أبعاد، ودلتاأناج{\displaystyle \delta ^{ij}}دلتا كرونكر ، [ 20 ]
  • توجد عملية معكوسة زمنياً وتخضع لنفس القوانين الديناميكية.

باستخدام التفكيكZ=X+أناY{\displaystyle Z=X+{\rm {i}}\,Y}وحقيقة أنج{\displaystyle C}إذا كان للدالة تباين محدود، فإن التباين التربيعي لـX{\displaystyle X}وY{\displaystyle Y}يُعطى بواسطة

(د[Xأنا،Xج]د[Xأنا،Yج]د[Yأنا،Xج]د[Yأنا،Yج])=|α|2مدلتاأناج(1+كوسϕالخطيئةϕالخطيئةϕ1-كوسϕ)دت{\displaystyle {\begin{pmatrix}d[X^{i},X^{j}]&d[X^{i},Y^{j}]\\d[Y^{i},X^{j}]&d[Y^{i},Y^{j}]\end{pmatrix}}={\frac {|\alpha |\,\hbar }{2\,m}}\,\delta ^{ij}\,{\begin{pmatrix}1+\cos \phi &\sin \phi \\\sin \phi &1-\cos \phi \end{pmatrix}}\,dt}

وبالتالي، من خلال توصيف ليفي للحركة البراونية،X(ت){\displaystyle X(t)}وY(ت){\displaystyle Y(t)}صف عمليتين مترابطتين من عمليات وينر مع انجراف موصوف بعملية التباين المحدودج(ت){\displaystyle C(t)}ثابت انتشار يتناسب مع|α|[0،){\displaystyle |\alpha |\in [0,\infty )}وارتباط يعتمد على الزاويةϕ{\displaystyle \phi }تكون العمليات مترابطة إلى أقصى حد في الحد الكمي، وترتبط بـϕ{-π2،π2}{\displaystyle \phi \in \left\{-{\frac {\pi }{2}},{\frac {\pi }{2}}\right\}}ومقابل ذلكαأنا×R{\displaystyle \alpha \in {\rm {i}}\times \mathbb {R} }بينما لا توجد علاقة بينهما في حالة الحد البراوني، المرتبطة بـϕ{0،π}{\displaystyle \phi \in \{0,\pi \}}ومقابل ذلكαR{\displaystyle \alpha \in \mathbb {R} }،

استُخدم مصطلح التكميم العشوائي لوصف عملية التكميم هذه في سبعينيات القرن العشرين. [ 21 ] أما اليوم، فيُشير مصطلح التكميم العشوائي بشكل أكثر شيوعًا إلى إطار عمل طوره باريسي وو في عام 1981. ونتيجةً لذلك، يُشار أحيانًا إلى عملية التكميم المُطورة في الميكانيكا العشوائية باسم تكميم نيلسون العشوائي أو التكميم العشوائي. [ 3 ]

سرعة العملية

العملية العشوائيةZ(ت){\displaystyle Z(t)}من شبه المؤكد أنها غير قابلة للتفاضل في أي مكان، بحيث تكون السرعةZ˙(ت)=دZ(ت)دت{\displaystyle {\dot {Z}}(t)={\frac {dZ(t)}{dt}}}إن مسار العملية غير محدد بدقة. ومع ذلك، توجد حقول سرعة، تُعرَّف باستخدام التوقعات الشرطية. ويتم تحديدها بواسطة

w+(x،ت)=ليمدت0هـ[Z(ت+دت)-Z(ت)دت|X(ت)=x]،{\displaystyle w_{+}(x,t)=\lim _{dt\rightarrow 0}\mathbb {E} \left[{\frac {Z(t+dt)-Z(t)}{dt}}\,{\Big |}\,X(t)=x\right],}
w-(x،ت)=ليمدت0هـ[Z(ت)-Z(ت-دت)دت|X(ت)=x]،{\displaystyle w_{-}(x,t)=\lim _{dt\rightarrow 0}\mathbb {E} \left[{\frac {Z(t)-Z(t-dt)}{dt}}\,{\Big |}\,X(t)=x\right],}ويمكن ربطها بتكامل إيتو على طول العمليةZ(ت){\displaystyle Z(t)}بما أن العملية غير قابلة للتفاضل، فإن هذه السرعات، بشكل عام، لا تتساوى. والتفسير الفيزيائي لهذه الحقيقة هو كما يلي: في أي لحظةت{\displaystyle t}يتعرض الجسيم لقوة عشوائية تغير سرعته بشكل فوري منw-{\displaystyle w_{-}}لw+{\displaystyle w_{+}}بما أن حقلي السرعة غير متساويين، فلا يوجد مفهوم فريد للسرعة لهذه العملية.Z(ت){\displaystyle Z(t)}في الواقع، أي سرعة معطاة بواسطة
wأ=أ w++(1-أ) w-{\displaystyle w_{a}=a\ w_{+}+(1-a)\ w_{-}}

معأ[0،1]{\displaystyle a\in [0,1]}يمثل هذا خيارًا صحيحًا لسرعة العمليةZ(ت){\displaystyle Z(t)}وينطبق هذا بشكل خاص على الحالة الخاصة.أ=12{\displaystyle a={\frac {1}{2}}}يرمز إليه بـw=w++w-2{\displaystyle w_{\circ }={\frac {w_{+}+w_{-}}{2}}}، والذي يمكن ربطه بتكامل ستراتونوفيتش على طولZ(ت){\displaystyle Z(t)}.

منذZ(ت){\displaystyle Z(t)}إذا كان له تغير تربيعي غير معدوم ، فيمكن أيضًا تعريف حقول سرعة من الدرجة الثانية معطاة بواسطة

w2،+(x،ت)=ليمدت0هـ[[Z(ت+دت)-Z(ت)][Z(ت+دت)-Z(ت)]دت | X(ت)=x]،{\displaystyle w_{2,+}(x,t)=\lim _{dt\rightarrow 0}\mathbb {E} \left[{\frac {[Z(t+dt)-Z(t)]\otimes [Z(t+dt)-Z(t)]}{dt}}\ {\Big |}\ X(t)=x\right],}
w2،-(x،ت)=ليمدت0هـ[[Z(ت)-Z(ت-دت)][Z(ت)-Z(ت-دت)]دت | X(ت)=x].{\displaystyle w_{2,-}(x,t)=\lim _{dt\rightarrow 0}\mathbb {E} \left[{\frac {[Z(t)-Z(t-dt)]\otimes [Z(t)-Z(t-dt)]}{dt}}\ {\Big |}\ X(t)=x\right].}

تفرض فرضية انعكاس الزمن علاقة بين هذين المجالين بحيثw2،±=±w2{\displaystyle w_{2,\pm }=\pm \,w_{2}}علاوة على ذلك، باستخدام علاقة البنية التي يتم من خلالها تحديد التغير التربيعي، نجد أنw2أناج(x،ت)=αمدلتاأناج{\displaystyle w_{2}^{ij}(x,t)={\frac {\alpha \,\hbar }{m}}\,\delta ^{ij}}ويترتب على ذلك أنه في صيغة ستراتونوفيتش، يتلاشى الجزء من الدرجة الثانية للسرعة، أيw2،=0{\displaystyle w_{2,\circ }=0}.

يُرمز إلى الجزء الحقيقي والجزء التخيلي من السرعات بـ

v=Rهـ(w)أندu=أنام(w).{\displaystyle v={\rm {Re}}(w)\qquad {\rm {and}}\qquad u={\rm {Im}}(w)\,.}

باستخدام وجود حقول السرعة هذه، يمكن تعريف عمليات السرعة بشكل رسميدبليو±(ت){\displaystyle W_{\pm }(t)}بواسطة إيتو المتكاملةدبليو±(ت) دت=د±Z(ت){\displaystyle \int W_{\pm }(t)\ dt=\int d_{\pm }Z(t)}وبالمثل، يمكن تعريف عملية ما بشكل رسميدبليو(ت){\displaystyle W_{\circ }(t)}بواسطة تكامل ستراتونوفيتشدبليو(ت) دت=دZ(ت){\displaystyle \int W_{\circ }(t)\ dt=\int d_{\circ }Z(t)}وعملية سرعة من الدرجة الثانيةدبليو2(ت){\displaystyle W_{2}(t)}بواسطة تكامل ستيلتجسدبليو2(ت) دت=د[Z،Z](ت){\displaystyle \int W_{2}(t)\ dt=\int d[Z,Z](t)}باستخدام علاقة البنية، نجد أن عملية السرعة من الدرجة الثانية معطاة بـدبليو2أناج(ت)=α م دلتاأناج{\displaystyle W_{2}^{ij}(t)={\frac {\alpha \ \hbar }{m}}\ \delta ^{ij}}ومع ذلك، فإن العملياتدبليو±(ت){\displaystyle W_{\pm }(t)}ودبليو(ت){\displaystyle W_{\circ }(t)}ليست محددة جيدًا: اللحظات الأولى موجودة ويتم تحديدها بواسطةهـ[دبليو±(ت) | X(ت)]=w±(X(ت)،ت){\displaystyle \mathbb {E} [W_{\pm }(t)\ |\ X(t)]=w_{\pm }(X(t),t)}لكن العزوم التربيعية تتباعد، أيهـ[||دبليو±(ت)||2 | X(ت)]={\displaystyle \mathbb {E} [||W_{\pm }(t)||^{2}\ |\ X(t)]=\infty }التفسير الفيزيائي لهذا التباين هو أنه في تمثيل الموضع يكون الموضع معروفًا بدقة، لكن السرعة لديها عدم يقين لا نهائي.

الفعل العشوائي

ينص شرط التكميم العشوائي على أن المسار العشوائي يجب أن يحقق أقصى قيمة لفعل عشوائيS=هـ[لدت]{\displaystyle S=\mathbb {E} \left[\int L\,dt\right]}، لكنه لا يحدد لاغرانجيان العشوائيل{\displaystyle L}يمكن الحصول على هذه الدالة اللاغرانجية من دالة لاغرانجية كلاسيكية باستخدام إجراء قياسي. هنا، نعتبر دالة لاغرانجية كلاسيكية على الشكل التالي:

ل(x،v،ت)=م2دلتاأناجvأناvج+qأأنا(x،ت)vأنا-يو(x،ت).{\displaystyle L(x,v,t)={\frac {m}{2}}\delta _{ij}v^{i}v^{j}+q\,A_{i}(x,t)\,v^{i}-{\mathfrak {U}}(x,t).}

هنا،(x،v){\displaystyle (x,v)}هي إحداثيات في فضاء الطور ( حزمة المماسدلتاأناج{\displaystyle \delta _{ij}}دالة كرونكر دلتا تصف المقياس علىRن{\displaystyle \mathbb {R} ^{n}}،م{\displaystyle m}يرمز إلى كتلة الجسيم،q{\displaystyle q}الشحنة تحت الجهد الشعاعيأ{\displaystyle A}، ويو{\displaystyle {\mathfrak {U}}}هو جهد قياسي. علاوة على ذلك، يُفترض استخدام اصطلاح جمع أينشتاين .

من الخصائص المهمة لهذا اللاغرانجي مبدأ ثبات القياس . ويمكن توضيح ذلك بتعريف فعل جديد.S~{\displaystyle {\tilde {S}}}من خلال إضافة حد مشتق كلي إلى الفعل الأصلي، بحيث

S~[x(ت)]=S[x(ت)]+دF(x،ت)=[م2دلتاأناجvأناvج+qأأناvأنا-يو+تF+vأناأناF]دت=[م2دلتاأناجvأناvج+qأ~أناvأنا-يو~]دت،{\displaystyle {\tilde {S}}[x(t)]=S[x(t)]+\int dF(x,t)=\int \left[{\frac {m}{2}}\delta _{ij}v^{i}v^{j}+qA_{i}v^{i}-{\mathfrak {U}}+\partial _{t}F+v^{i}\partial _{i}F\right]dt=\int \left[{\frac {m}{2}}\delta _{ij}v^{i}v^{j}+q{\tilde {A}}_{i}v^{i}-{\tilde {\mathfrak {U}}}\right]dt,}

أينأ~أنا=أأنا+q-1أناF{\displaystyle {\tilde {A}}_{i}=A_{i}+q^{-1}\partial _{i}F}ويو~=يو-تF{\displaystyle {\tilde {\mathfrak {U}}}={\mathfrak {U}}-\partial _{t}F}وبالتالي، بما أن الديناميكيات لا تتأثر بإضافة مشتقة كلية إلى الفعل، فإن الفعل يظل ثابتًا بالنسبة للمقياس في ظل إعادة تعريف الكمونات المذكورة أعلاه لدالة قابلة للتفاضل بشكل عشوائيF{\displaystyle F}.

لإنشاء لاغرانجي عشوائي يُطابق هذا اللاغرانجي الكلاسيكي، يجب البحث عن امتداد أدنى للاغرانجي المذكور أعلاه يحافظ على ثبات القياس هذا. [ 22 ] في صياغة ستراتونوفيتش للنظرية، يمكن القيام بذلك مباشرةً، لأن المؤثر التفاضلي في صياغة ستراتونوفيتش يُعطى بواسطة

دF(x،ت)=(تF+vأناأناF)دت.{\displaystyle \int d_{\circ }F(x,t)=\int \left(\partial _{t}F+v_{\circ }^{i}\partial _{i}F\right)dt\,.}

لذلك، يمكن الحصول على لاغرانجيان ستراتونوفيتش عن طريق استبدال السرعة الكلاسيكيةv{\displaystyle v}بواسطة السرعة المركبةw{\displaystyle w_{\circ }}بحيث

ل(x،w،ت)=م2دلتاأناجwأناwج+qأأناwأنا-يو.{\displaystyle L_{\circ }(x,w_{\circ },t)={\frac {m}{2}}\delta _{ij}w_{\circ }^{i}w_{\circ }^{j}+qA_{i}w_{\circ }^{i}-{\mathfrak {U}}\,.}

في صيغة إيتو، تصبح الأمور أكثر تعقيداً، حيث يتم إعطاء المشتقة الكلية بواسطة ليمّة إيتو :د±F(x،ت)=(تF+v±أناأناF±12v2أناججأناF)دت.{\displaystyle \int d_{\pm }F(x,t)=\int \left(\partial _{t}F+v_{\pm }^{i}\partial _{i}F\pm {\frac {1}{2}}v_{2}^{ij}\partial _{j}\partial _{i}F\right)dt\,.}

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

ل±(x،w±،w2،ت)=م2دلتاأناجw±أناw±ج+qأأناw±أنا±q2جأأناw2أناج-يو.{\displaystyle L_{\pm }(x,w_{\pm },w_{2},t)={\frac {m}{2}}\delta _{ij}w_{\pm }^{i}w_{\pm }^{j}+q\,A_{i}w_{\pm }^{i}\pm {\frac {q}{2}}\partial _{j}A_{i}w_{2}^{ij}-{\mathfrak {U}}\,.}

يمكن تعريف الفعل العشوائي باستخدام لاغرانجيان ستراتونوفيتش، وهو يساوي الفعل المحدد بواسطة لاغرانجيان إيتو حتى حد متباعد: [ 2 ]

S=هـ[لدت]=هـ[ل±دت]±هـ[لدت].{\displaystyle S=\mathbb {E} \left[\int L_{\circ }\,dt\right]=\mathbb {E} \left[\int L_{\pm }dt\right]\pm \mathbb {E} \left[\int L_{\infty }dt\right]\,.}

يمكن حساب الحد المتباعد ويُعطى بالصيغة التالية:

هـ[لدت]=م2γدلتاأناجw2أناجتدت=απأناأنا=1نكأنا،{\displaystyle \mathbb {E} \left[\int L_{\infty }dt\right]={\frac {m}{2}}\oint _{\gamma }{\frac {\delta _{ij}w_{2}^{ij}}{t}}dt=\alpha \,\hbar \,\pi \,{\rm {i}}\,\sum _{i=1}^{n}k_{i},}

أينكأناZ{\displaystyle k_{i}\in \mathbb {Z} }هي أرقام متعرجة تحسب مدى تعرج المسارγ(ت){\displaystyle \gamma (t)}حول العمود عندت=0{\displaystyle t=0}[ 23 ]

بما أن الحد المتباعد ثابت، فإنه لا يُسهم في معادلات الحركة. ولهذا السبب، تم إهمال هذا الحد في الدراسات المبكرة للميكانيكا العشوائية. [ 2 ] مع ذلك، عند إهمال هذا الحد، لا تستطيع الميكانيكا العشوائية تفسير ظهور الأطياف المنفصلة في ميكانيكا الكم. تُعرف هذه المشكلة بنقد والستروم، [ 24 ] [ 25 ] ويمكن حلها من خلال مراعاة الحد المتباعد بشكل صحيح.

يوجد أيضًا صياغة هاميلتونية للميكانيكا العشوائية. [ 22 ] [ 26 ] وهي تبدأ من تعريف الزخم الكنسي:

ص،أنا=لwأنا=مدلتاأناجwج+qأأنا،{\displaystyle p_{\circ ,i}={\frac {\partial L_{\circ }}{\partial w_{\circ }^{i}}}=m\,\delta _{ij}w_{\circ }^{j}+q\,A_{i}\,,}
ص±،أنا=ل±w±أنا=مدلتاأناجw±ج+qأأنا.{\displaystyle p_{\pm ,i}={\frac {\partial L_{\pm }}{\partial w_{\pm }^{i}}}=m\,\delta _{ij}w_{\pm }^{j}+q\,A_{i}\,.}

يمكن بعد ذلك الحصول على الهاميلتوني في صيغة ستراتونوفيتش عن طريق تحويل ليجندر من الدرجة الأولى :

ح(x،ص،ت)=ص،أناvأنا-ل(x،v،ت).{\displaystyle H_{\circ }(x,p_{\circ },t)=p_{\circ ,i}v_{\circ }^{i}-L(x,v_{\circ },t)\,.}

من ناحية أخرى، في صياغة إيتو، يتم الحصول على الهاميلتوني من خلال تحويل ليجندر من الدرجة الثانية: [ 27 ]

ح±(x،ص±،ص±،ت)=صأنا±w±أنا±12w2أناجصأناج±-ل(x،w±،w2،ت).{\displaystyle H_{\pm }(x,p^{\pm },\partial p^{\pm },t)=p_{i}^{\pm }w_{\pm }^{i}\pm {\frac {1}{2}}w_{2}^{ij}\partial p_{ij}^{\pm }-L(x,w_{\pm },w_{2},t)\,.}

معادلات أويلر-لاغرانج

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

د(لwأنا)=(لxأنا)دت.{\displaystyle \int d_{\circ }\left({\frac {\partial L_{\circ }}{\partial w_{\circ }^{i}}}\right)=\int \left({\frac {\partial L_{\circ }}{\partial x^{i}}}\right)dt\,.}

بالنسبة للدالة اللاغرانجية، التي نوقشت في القسم السابق، يؤدي هذا إلى معادلة تفاضلية عشوائية من الدرجة الثانية على النحو التالي وفقًا لمعنى ستراتونوفيتش :

مدلتاأناجد2Zج(ت)=qFأناج(X(ت)،ت)دZج(ت)دت-qتأأنا(X(ت)،ت)دت2-أنايو(X(ت)،ت)دت2،{\displaystyle m\,\delta _{ij}\,d_{\circ }^{2}Z^{j}(t)=q\,F_{ij}(X(t),t)\,d_{\circ }Z^{j}(t)\,dt-q\,\partial _{t}A_{i}(X(t),t)\,dt^{2}-\partial _{i}{\mathfrak {U}}(X(t),t)\,dt^{2}\,,}

حيث تُعطى شدة المجال بواسطةFأناج:=أناأج-جأأنا{\displaystyle F_{ij}:=\partial _{i}A_{j}-\partial _{j}A_{i}}. هذه المعادلة بمثابة نسخة عشوائية من قانون نيوتن الثاني .

في صيغة إيتو، تُعطى معادلات أويلر-لاغرانج العشوائية على النحو التالي

د±(ل±w±أنا)=(ل±xأنا)دت.{\displaystyle \int d_{\pm }\left({\frac {\partial L_{\pm }}{\partial w_{\pm }^{i}}}\right)=\int \left({\frac {\partial L_{\pm }}{\partial x^{i}}}\right)dt\,.}

يؤدي هذا إلى معادلة تفاضلية عشوائية من الدرجة الثانية بالمعنى الذي وضعه إيتو ، معطاة بصيغة عشوائية لقانون نيوتن الثاني على النحو التالي:

مدلتاأناجد±2Zج(ت)=qFأناج(X(ت)،ت)د±Zج(ت)دت±αq2مدلتاجككFأناج(X(ت)،ت)دت2-qتأأنا(X(ت)،ت)دت2-أنايو(X(ت)،ت)دت2.{\displaystyle m\,\delta _{ij}\,d_{\pm }^{2}Z^{j}(t)=q\,F_{ij}(X(t),t)\,d_{\pm }Z^{j}(t)\,dt\pm {\frac {\alpha \,\hbar \,q}{2\,m}}\,\delta ^{jk}\,\partial _{k}F_{ij}(X(t),t)\,dt^{2}-q\,\partial _{t}A_{i}(X(t),t)\,dt^{2}-\partial _{i}{\mathfrak {U}}(X(t),t)\,dt^{2}\,.}

معادلات هاميلتون-جاكوبي

يمكن أيضًا الحصول على معادلات الحركة في تعميم عشوائي لصياغة هاميلتون-جاكوبي للميكانيكا الكلاسيكية. في هذه الحالة، يبدأ المرء بتعريف دالة هاميلتون الرئيسية. بالنسبة لدالة لاغرانجل+{\displaystyle L_{+}}تُعرَّف هذه الدالة على النحو التالي:

S+(x،ت؛xو،تو):=-هـ[تتول+(X(s)،دبليو+(s)،دبليو2(s)،s)دs|X(ت)=x،X(تو)=xو]،{\displaystyle S_{+}(x,t;x_{f},t_{f}):=-\mathbb {E} \left[\int _{t}^{t_{f}}L_{+}(X(s),W_{+}(s),W_{2}(s),s)\,ds\,{\Big |}\,X(t)=x,X(t_{f})=x_{f}\right],}

حيث يُفترض أن العملية{X(s):s[ت،تو]}{\displaystyle \{X(s):s\in [t,t_{f}]\}}تخضع لمعادلات أويلر-لاغرانج العشوائية. وبالمثل، بالنسبة للدالة اللاغرانجيةل-{\displaystyle L_{-}}تُعرَّف دالة هاميلتون الرئيسية على النحو التالي:

S-(x،ت؛x0،ت0):=هـ[ت0تل-(X(s)،دبليو-(s)،دبليو2(s)،s)دs|X(ت)=x،X(ت0)=x0]،{\displaystyle S_{-}(x,t;x_{0},t_{0}):=\mathbb {E} \left[\int _{t_{0}}^{t}L_{-}(X(s),W_{-}(s),W_{2}(s),s)\,ds\,{\Big |}\,X(t)=x,X(t_{0})=x_{0}\right],}

حيث يُفترض أن العملية{X(s):s[ت0،ت]}{\displaystyle \{X(s):s\in [t_{0},t]\}}تخضع لمعادلات أويلر-لاغرانج العشوائية. ونظرًا للجزء المتباعد من الفعل، تخضع هذه الدوال الرئيسية لعلاقة التكافؤ.

S~±S±أناوكZن،suجحتحأتS~±=S±±απأناأنا=1نكأنا.{\displaystyle {\tilde {S}}_{\pm }\equiv S_{\pm }\;{\rm {if}}\;\exists k\in \mathbb {Z} ^{n},\;{\rm {such\,that}}\;{\tilde {S}}_{\pm }=S_{\pm }\pm \alpha \,\pi \,{\rm {i}}\,\hbar \,\sum _{i=1}^{n}k_{i}\,.}

عن طريق تغيير الدوال الرئيسية بالنسبة للنقطة(x،ت){\displaystyle (x,t)}يتم إيجاد معادلات هاميلتون-جاكوبي. وهي معطاة بواسطة

{xأناS±(x،ت)=صأنا±(x،ت)،تS±(x،ت)=-ح±(x،ص±(x،ت)،ص±(x،ت)،ت).{\displaystyle {\begin{cases}{\frac {\partial }{\partial x^{i}}}S_{\pm }(x,t)&=p_{i}^{\pm }(x,t)\,,\\{\frac {\partial }{\partial t}}S_{\pm }(x,t)&=-H_{\pm }(x,p_{\pm }(x,t),\partial p_{\pm }(x,t),t)\,.\end{cases}}}

لاحظ أن هذه المعادلات تبدو مماثلة للحالة الكلاسيكية. مع ذلك، يُحسب الهاميلتوني في معادلة هاميلتون-جاكوبي الثانية باستخدام تحويل ليجندر من الرتبة الثانية. علاوة على ذلك، ونظرًا للجزء المتباعد من الفعل، توجد معادلة هاميلتون-جاكوبي ثالثة، تأخذ شكل قيد تكاملي غير تافه.

(صأنا±v±أنا±12v2أناجأناصج)دت=±απأناأنا=1نكأنا.{\displaystyle \oint \left(p_{i}^{\pm }\,v_{\pm }^{i}\pm {\frac {1}{2}}v_{2}^{ij}\,\partial _{i}p_{j}\right)dt=\pm \alpha \,\hbar \,\pi \,{\rm {i}}\,\sum _{i=1}^{n}k_{i}\,.}

بالنسبة للدالة اللاغرانجية المعطاة، فإن أول معادلتين من معادلات هاميلتون-جاكوبي تعطيان

{أناS=مدلتاأناجw±ج+qأأنا،تS=-م2دلتاأناجw±أناw±جم2دلتاأناجw2أناككw±ج-يو.{\displaystyle {\begin{cases}\partial _{i}S&=m\,\delta _{ij}w_{\pm }^{j}+q\,A_{i}\,,\\\partial _{t}S&=-{\frac {m}{2}}\,\delta _{ij}w_{\pm }^{i}w_{\pm }^{j}\mp {\frac {m}{2}}\,\delta _{ij}w_{2}^{ik}\partial _{k}w_{\pm }^{j}-{\mathfrak {U}}\,.\end{cases}}}

يمكن دمج هاتين المعادلتين، مما ينتج عنه

[مدلتاأناج(ت+w±كك±12w2كللك)-qFأناج]w±ج=±q2w2جككFأناج-qتأأنا-أنايو.{\displaystyle \left[m\,\delta _{ij}\left(\partial _{t}+w_{\pm }^{k}\partial _{k}\pm {\frac {1}{2}}\,w_{2}^{kl}\partial _{l}\partial _{k}\right)-q\,F_{ij}\right]w_{\pm }^{j}=\pm {\frac {q}{2}}\,w_{2}^{jk}\partial _{k}F_{ij}-q\,\partial _{t}A_{i}-\partial _{i}{\mathfrak {U}}\,.}

باستخدام ذلكw2أناج=αمدلتاأناج{\displaystyle w_{2}^{ij}={\frac {\alpha \,\hbar }{m}}\,\delta ^{ij}}، هذه المعادلة، الخاضعة لشرط التكامل والشرط الابتدائيw+(x،ت0)=w0(x){\displaystyle w_{+}(x,t_{0})=w_{0}(x)}أو حالة نهائيةw-(x،تو)=wو(x){\displaystyle w_{-}(x,t_{f})=w_{f}(x)}يمكن حلها لـw±(x،ت){\displaystyle w_{\pm }(x,t)}ويمكن بعد ذلك إدخال الحل في معادلة إيتو.

{د±Zأنا(ت)=w±أنا(x،ت)دت+دمأنا(ت)،د[مأنا،مج](ت)=αمدلتاأناجدت،{\displaystyle {\begin{cases}d_{\pm }Z^{i}(t)&=w_{\pm }^{i}(x,t)\,dt+dM^{i}(t)\,,\\d[M^{i},M^{j}](t)&={\frac {\alpha \,\hbar }{m}}\delta ^{ij}\,dt\,,\end{cases}}}

والتي يمكن حلها في هذه العملية{Z(ت):ت[ت0،تو]}{\displaystyle \{Z(t):t\in [t_{0},t_{f}]\}}وبالتالي، عندما تكون الحالة الأوليةX(ت0)=x0{\displaystyle X(t_{0})=x_{0}}(للمعادلة الموجهة المستقبلية التي تحمل علامة+{\displaystyle +}) أو حالة نهائيةX(تو)=xو{\displaystyle X(t_{f})=x_{f}}(للمعادلة الموجهة الماضية التي تحمل علامة-{\displaystyle -}إذا تم تحديد ) ، فسيتم العثور على عملية عشوائية فريدة{X(ت):ت[ت0،تو]}{\displaystyle \{X(t):t\in [t_{0},t_{f}]\}}وهذا يصف مسار الجسيم.

معادلة الانتشار

تتمثل النتيجة الرئيسية للميكانيكا العشوائية لنيلسون في أنها تستنتج معادلة شرودنغر من العملية العشوائية المفترضة. وفي هذا الاشتقاق، تُستخدم معادلات هاميلتون-جاكوبي.

{xأناS±(x،ت)=صأنا±(x،ت)تS±(x،ت)=-ح±(x،ص±(x،ت)،ص±(x،ت)،ت){\displaystyle {\begin{cases}{\frac {\partial }{\partial x^{i}}}S_{\pm }(x,t)&=p_{i}^{\pm }(x,t)\\{\frac {\partial }{\partial t}}S_{\pm }(x,t)&=-H_{\pm }(x,p_{\pm }(x,t),\partial p_{\pm }(x,t),t)\end{cases}}}

يتم دمجها، بحيث نحصل على المعادلة

2م(تS±+يو)+دلتاأناج(أناS±جS±±αجأناS±-2qأأناجS±αqجأأنا+q2أأناأج)=0.{\displaystyle 2\,m\left(\partial _{t}S_{\pm }+{\mathfrak {U}}\right)+\delta ^{ij}\left(\partial _{i}S_{\pm }\partial _{j}S_{\pm }\pm \alpha \,\hbar \,\partial _{j}\partial _{i}S_{\pm }-2\,q\,A_{i}\partial _{j}S_{\pm }\mp \alpha \,\hbar \,q\,\partial _{j}A_{i}+q^{2}\,A_{i}A_{j}\right)=0\,.}

بعد ذلك، يتم تعريف الدالة الموجية

Ψ±(x،ت)=خبرة(±S±(x،ت)α).{\displaystyle \Psi _{\pm }(x,t)=\exp \left(\pm {\frac {S_{\pm }(x,t)}{\alpha \,\hbar }}\right).}

بما أن الدوال الرئيسية لهاملتون متعددة القيم، فإن الدوال الموجية تخضع لعلاقات التكافؤ.

Ψ~+Ψ+أناوΨ~+=±Ψ+أندΨ~-Ψ-أناوΨ~-=±Ψ-.{\displaystyle {\tilde {\Psi }}_{+}\equiv \Psi _{+}\quad {\rm {if}}\quad {\tilde {\Psi }}_{+}=\pm \Psi _{+}\qquad {\rm {and}}\qquad {\tilde {\Psi }}_{-}\equiv \Psi _{-}\quad {\rm {if}}\quad {\tilde {\Psi }}_{-}=\pm \Psi _{-}\,.}

علاوة على ذلك، تخضع الدوال الموجية لمعادلات الانتشار المعقدة

-αتΨ+=[دلتاأناج2م(αxأنا+qأأنا)(αxج+qأج)+يو]Ψ+،{\displaystyle -\alpha \,\hbar \,{\frac {\partial }{\partial t}}\Psi _{+}=\left[{\frac {\delta ^{ij}}{2\,m}}\left(\alpha \,\hbar \,{\frac {\partial }{\partial x^{i}}}+q\,A_{i}\right)\left(\alpha \,\hbar \,{\frac {\partial }{\partial x^{j}}}+q\,A_{j}\right)+{\mathfrak {U}}\right]\Psi _{+}\,,}
αتΨ-=[دلتاأناج2م(αxأنا+qأأنا)(αxج+qأج)+يو]Ψ-.{\displaystyle \alpha \,\hbar \,{\frac {\partial }{\partial t}}\Psi _{-}=\left[{\frac {\delta ^{ij}}{2\,m}}\left(\alpha \,\hbar \,{\frac {\partial }{\partial x^{i}}}+q\,A_{i}\right)\left(\alpha \,\hbar \,{\frac {\partial }{\partial x^{j}}}+q\,A_{j}\right)+{\mathfrak {U}}\right]\Psi _{-}\,.}

وبالتالي، لأي عملية تحقق مسلمات الميكانيكا العشوائية، يمكن بناء دالة موجية تخضع لمعادلات الانتشار هذه. وبفضل علاقات التكافؤ على دالة هاميلتون الرئيسية، فإن العكس صحيح أيضًا: لأي حل لمعادلات الانتشار المعقدة هذه، يمكن بناء عملية عشوائية.{X(ت):ت[ت0،تو]}{\displaystyle \{X(t):t\in [t_{0},t_{f}]\}}هذا حلٌّ لمسلمات الميكانيكا العشوائية. وقد تم التوصل إلى نتيجة مماثلة من خلال نظرية فاينمان-كاك .

وأخيرًا، يمكن للمرء أن يبني دالة كثافة احتمالية

ρ±(x،ت):=|Ψ±(x،ت)|2|Ψ±(y،ت)|2دنy،{\displaystyle \rho _{\pm }(x,t):={\frac {|\Psi _{\pm }(x,t)|^{2}}{\int |\Psi _{\pm }(y,t)|^{2}d^{n}y}}\,,}

والتي تصف احتمالات الانتقال للعملية{X(ت):ت[ت0،تو]}{\displaystyle \{X(t):t\in [t_{0},t_{f}]\}}وبشكل أدق،ρ+{\displaystyle \rho _{+}}يصف احتمالية التواجد في الحالة(x،ت){\displaystyle (x,t)}بافتراض أن النظام ينتهي في الحالة(xو،تو){\displaystyle (x_{f},t_{f})}لذلك، فإن معادلة الانتشار لـΨ+{\displaystyle \Psi _{+}}يمكن تفسيرها على أنها معادلة كولموغوروف العكسية للعملية{X(ت):ت[ت0،تو]}{\displaystyle \{X(t):t\in [t_{0},t_{f}]\}}. بصورة مماثلة، ρ-{\displaystyle \rho _{-}}يصف احتمالية التواجد في الحالة(x،ت){\displaystyle (x,t)}بافتراض أن النظام ينتهي في الحالة(x0،ت0){\displaystyle (x_{0},t_{0})}، عندما يتم تطويرها عكسيًا في الزمن. لذلك، فإن معادلة الانتشار لـΨ-{\displaystyle \Psi _{-}}يمكن تفسيرها على أنها معادلة كولموغوروف العكسية للعملية{X(ت):ت[ت0،تو]}{\displaystyle \{X(t):t\in [t_{0},t_{f}]\}}عندما يتطور نحو الماضي. بعكس اتجاه الزمن، يجد المرء أنρ-{\displaystyle \rho _{-}}يصف احتمالية التواجد في الحالة(x،ت){\displaystyle (x,t)}بافتراض أن النظام يبدأ في الحالة(x0،ت0){\displaystyle (x_{0},t_{0})}، عندما يتطور مع مرور الوقت. وبالتالي، فإن معادلة الانتشار لـΨ-{\displaystyle \Psi _{-}}ويمكن تفسيرها أيضًا على أنها معادلة كولموغوروف الأمامية للعملية{X(ت):ت[ت0،تو]}{\displaystyle \{X(t):t\in [t_{0},t_{f}]\}}عندما يتطور نحو المستقبل.

الجوانب الرياضية

الحالات المحددة

تتضمن النظرية العديد من الحدود الخاصة:

  • النهاية الكلاسيكية معα=0{\displaystyle \alpha =0}في هذه الحالة، العمليةX{\displaystyle X}والعملية المساعدةY{\displaystyle Y}يصف مسارين حتميين منفصلين.
  • حد براون معα(0،){\displaystyle \alpha \in (0,\infty )}في هذه الحالة، العمليةX{\displaystyle X}يصف عملية وينر (المعروفة أيضًا باسم الحركة البراونية ) التي تم إثبات النتيجة المذكورة أعلاه لها بواسطة نظرية فاينمان-كاك ، بينما العملية المساعدةY{\displaystyle Y}يصف عملية حتمية.
  • الحد الكمي معαأنا×(0،){\displaystyle \alpha \in {\rm {i}}\times (0,\infty )}في هذه الحالة، العمليةX{\displaystyle X}والعملية المساعدةY{\displaystyle Y}صف عمليتين من عمليات وينر مرتبطتين ارتباطًا إيجابيًا.
  • الحد البراوني المعكوس زمنيًا معα(-،0){\displaystyle \alpha \in (-\infty ,0)}في هذه الحالة، العمليةX{\displaystyle X}يصف عملية حتمية، بينما العملية المساعدةY{\displaystyle Y}يصف عملية وينر .
  • الحد الكمي المعكوس زمنيًا معαأنا×(-،0){\displaystyle \alpha \in {\rm {i}}\times (-\infty ,0)}في هذه الحالة، العمليةX{\displaystyle X}والعملية المساعدةY{\displaystyle Y}صف عمليتين من عمليات وينر ذات ارتباط سلبي.

في حالة الحد البراوني مع الشرط الابتدائيu(x،ت0)=0{\displaystyle u(x,t_{0})=0}أو حالة نهائيةu(x،تو)=0{\displaystyle u(x,t_{f})=0}(وهذا يعنيu(x،ت)=0ت{\displaystyle u(x,t)=0\quad \forall t}) العملياتX{\displaystyle X}وY{\displaystyle Y}يتم فصلها، بحيث تكون ديناميكيات العملية المساعدةY{\displaystyle Y}يمكن التخلص منها، وX{\displaystyle X}يتم وصفها بواسطة عملية وينر حقيقية . في جميع الحالات الأخرى معα0{\displaystyle \alpha \neq 0}، وتكون العمليات مرتبطة ببعضها البعض، بحيث تكون العملية المساعدةY{\displaystyle Y}يجب أخذ ذلك في الاعتبار عند استخلاص ديناميكياتX{\displaystyle X}.

تناظر انعكاس الزمن

تكون النظرية متناظرة تحت عملية عكس الزمن(ت،α،q)(-ت،-α،-q){\displaystyle (t,\alpha ,q)\leftrightarrow (-t,-\alpha ,-q)}.

في حدود الحركة البراونية، تكون النظرية مُبددة للطاقة إلى أقصى حد ، بينما تكون الحدود الكمومية وحدوية ، بحيث

ددت|Ψ±(y،ت)|2دنy|αأنا×R=0.{\displaystyle {\frac {d}{dt}}\int |\Psi _{\pm }(y,t)|^{2}d^{n}y{\Big |}_{\alpha \in {\rm {i}}\times \mathbb {R} }=0\,.}

علاقات التبادل الأساسية

يمكن إعادة كتابة معادلة الانتشار على النحو التالي

αتΨ±=ح^(x^،ص^±،ت)Ψ±،{\displaystyle \mp \alpha \,\hbar \,{\frac {\partial }{\partial t}}\Psi _{\pm }={\hat {H}}({\hat {x}},{\hat {p}}^{\pm },t)\,\Psi _{\pm }\,,}

أينح^{\displaystyle {\hat {H}}}هو مؤثر هاميلتوني . وهذا يسمح بإدخال مؤثرات الموضع والزخم كـ

x^أنا=xأناأندص^أنا±=±αxأنا،{\displaystyle {\hat {x}}^{i}=x^{i}\qquad {\rm {and}}\qquad {\hat {p}}_{i}^{\pm }=\pm \alpha \,\hbar \,{\frac {\partial }{\partial x^{i}}}\,,}

بحيث يكون للهاميلتوني شكله المألوف

ح(x^،ص^،ت)=دلتاأناج2م[ص^أنا-qأأنا(x^،ت)][ص^ج-qأج(x^،ت)]+يو(x^،ت).{\displaystyle H({\hat {x}},{\hat {p}},t)={\frac {\delta ^{ij}}{2\,m}}\,{\Big [}{\hat {p}}_{i}-q\,A_{i}({\hat {x}},t){\Big ]}{\Big [}{\hat {p}}_{j}-q\,A_{j}({\hat {x}},t){\Big ]}+{\mathfrak {U}}({\hat {x}},t)\,.}

تخضع هذه المؤثرات لعلاقة التبادل المتعارف عليها

[x^أنا،ص^ج±]=αدلتاجأنا.{\displaystyle [{\hat {x}}^{i},{\hat {p}}_{j}^{\pm }]=\mp \alpha \,\hbar \,\delta _{j}^{i}\,.}

انظر أيضاً

ملحوظات

مراجع

أوراق

الكتب