مصفوفة الكثافة

في ميكانيكا الكم ، تُستخدم مصفوفة الكثافة (أو مؤثر الكثافة ) لحساب احتمالات نتائج القياسات التي تُجرى على الأنظمة الفيزيائية . [ 1 ] وهي تعميم لمتجهات الحالة أو الدوال الموجية : فبينما لا تستطيع هذه الأخيرة تمثيل سوى الحالات النقية ، تستطيع مصفوفات الكثافة تمثيل مجموعات مختلطة من الحالات. [ 2 ] : 73 [ 3 ] : 100 وتظهر هذه المصفوفات في ميكانيكا الكم في حالتين مختلفتين:

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

وبالتالي فإن مصفوفات الكثافة هي أدوات حاسمة في مجالات ميكانيكا الكم التي تتعامل مع الحالات المختلطة (لا ينبغي الخلط بينها وبين الحالات المتراكبة )، مثل ميكانيكا الكم الإحصائية ، وأنظمة الكم المفتوحة ، ومعلومات الكم .

التعريف والدافع

مصفوفة الكثافة هي تمثيل لمؤثر خطي يُسمى مؤثر الكثافة . تُستخلص مصفوفة الكثافة من مؤثر الكثافة باختيار أساس متعامد في الفضاء الأساسي. [ 4 ] عمليًا، يُستخدم مصطلحا مصفوفة الكثافة ومؤثر الكثافة غالبًا بشكل متبادل.

اختر أساسًا مع الولايات|0{\displaystyle |0\rangle }،|1{\displaystyle |1\rangle }في فضاء هيلبرت ثنائي الأبعاد ، يتم تمثيل عامل الكثافة بواسطة المصفوفة (ρأناج)=(ρ٠٠ρ01ρ10ρ11)=(ص0ρ01ρ01*ص1){\displaystyle (\rho _{ij})=\left({\begin{matrix}\rho _{00}&\rho _{01}\\\rho _{10}&\rho _{11}\end{matrix}}\right)=\left({\begin{matrix}p_{0}&\rho _{01}\\\rho _{01}^{*}&p_{1}\end{matrix}}\right)} حيث تكون العناصر القطرية أعدادًا حقيقية مجموعها يساوي واحدًا (وتسمى أيضًا عدد سكان الولايتين).|0{\displaystyle |0\rangle }،|1{\displaystyle |1\rangle }العناصر غير القطرية هي مرافقات معقدة لبعضها البعض (وتسمى أيضًا التماسك)؛ وهي مقيدة في قيمتها بالشرط التالي:(ρأناج){\displaystyle (\rho _{ij})}ليكن عاملًا شبه محدد موجبًا ، انظر أدناه.

مؤثر الكثافة هو مؤثر شبه موجب ، ذاتي الترافق، أثره يساوي واحدًا ، ويعمل على فضاء هيلبرت للنظام. [ 5 ] [ 6 ] [ 7 ] يمكن تبرير هذا التعريف بالنظر إلى حالة توجد فيها بعض الحالات النقية|ψج{\displaystyle |\psi _{j}\rangle }(التي ليست بالضرورة متعامدة) يتم إعدادها باحتماليةصج{\displaystyle p_{j}}كل منها. [ 8 ] يُعرف هذا باسم مجموعة من الحالات النقية. احتمال الحصول على نتيجة قياس إسقاطيم{\displaystyle m}عند استخدام أجهزة العرضΠم{\displaystyle \Pi _{m}}يُعطى بواسطة [ 3 ] : 99ص(م)=جصجψج|Πم|ψج،{\displaystyle p(m)=\sum _{j}p_{j}\left\langle \psi _{j}\right|\Pi _{m}\left|\psi _{j}\right\rangle ,} والتي يمكن إثبات أنها تساوي ص(م)=tr[Πم(جصج|ψجψج|)].{\displaystyle p(m)=\operatorname {tr} \left[\Pi _{m}\left(\sum _{j}p_{j}\left|\psi _{j}\right\rangle \left\langle \psi _{j}\right|\right)\right].} وبالتالي، فإن عامل الكثافة ، المعرّف على النحو التالي ρ=جصج|ψجψج|،{\displaystyle \rho =\sum _{j}p_{j}\left|\psi _{j}\right\rangle \left\langle \psi _{j}\right|,} يمثل هذا تمثيلاً مناسباً لحالة هذه المجموعة. هذا المؤثر شبه موجب، وذاتي الترافق، وله أثر يساوي واحد. وبالعكس، يترتب من نظرية الطيف أن كل مؤثر بهذه الخصائص يمكن كتابته على النحو التالي:جصج|ψجψج|{\textstyle \sum _{j}p_{j}\left|\psi _{j}\right\rangle \left\langle \psi _{j}\right|}بالنسبة لبعض الولايات|ψج{\displaystyle \left|\psi _{j}\right\rangle }والمعاملاتصج{\displaystyle p_{j}}الأعداد غير السالبة التي مجموعها يساوي واحدًا. [ 9 ] [ 3 ] : 102 ومع ذلك، لن يكون هذا التمثيل فريدًا، كما هو موضح في نظرية شرودنغر-HJW .

ثمة دافع آخر لتعريف مؤثرات الكثافة، وهو النظر في القياسات المحلية على الحالات المتشابكة.|Ψ{\displaystyle |\Psi \rangle }ليكن حالة متشابكة نقية في فضاء هيلبرت المركبح1ح2{\displaystyle {\mathcal {H}}_{1}\otimes {\mathcal {H}}_{2}}احتمالية الحصول على نتيجة القياسم{\displaystyle m}عند قياس أجهزة العرضΠم{\displaystyle \Pi _{m}}في مساحة هيلبرتح1{\displaystyle {\mathcal {H}}_{1}}وحده معطى بواسطة [ 3 ] : 107ص(م)=Ψ|(Πمأنا)|Ψ،{\displaystyle p(m)=\left\langle \Psi \right|\left(\Pi _{m}\otimes I\right)\left|\Psi \right\rangle ,} والتي تصبح بعد المعالجة الجبرية ص(م)=tr[Πم(tr2|ΨΨ|)]،{\displaystyle p(m)=\operatorname {tr} \left[\Pi _{m}\left(\operatorname {tr} _{2}\left|\Psi \right\rangle \left\langle \Psi \right|\right)\right],} أينtr2{\displaystyle \operatorname {tr} _{2}}يشير إلى الأثر الجزئي على فضاء هيلبرتح2{\displaystyle {\mathcal {H}}_{2}}وهذا يجعل المشغل ρ=tr2|ΨΨ|{\displaystyle \rho =\operatorname {tr} _{2}\left|\Psi \right\rangle \left\langle \Psi \right|} أداة ملائمة لحساب احتمالات هذه القياسات المحلية. يتمتع هذا المؤثر بجميع خصائص مؤثر الكثافة، ويُعرف باسم مصفوفة الكثافة المختزلة لـ|Ψ{\displaystyle |\Psi \rangle }في النظام الفرعي 1. وعلى العكس من ذلك، تنص نظرية شرودنغر-هاميلتون-جونز-ويليامز على أنه يمكن كتابة جميع معاملات الكثافة على النحو التالي:tr2|ΨΨ|{\displaystyle \operatorname {tr} _{2}\left|\Psi \right\rangle \left\langle \Psi \right|}بالنسبة لبعض الولايات|Ψ{\displaystyle \left|\Psi \right\rangle }.

الحالات النقية والمختلطة

الحالة الكمومية النقية هي حالة لا يمكن كتابتها كمزيج احتمالي، أو توليفة محدبة ، من حالات كمومية أخرى. [ 7 ] توجد عدة توصيفات مكافئة للحالات النقية في لغة مؤثرات الكثافة. [ 2 ] : 73 يمثل مؤثر الكثافة حالة نقية إذا وفقط إذا:

  • يمكن كتابتها كحاصل ضرب خارجي لمتجه حالة|ψ{\displaystyle |\psi \rangle }مع نفسه، أيρ=|ψψ|.{\displaystyle \rho =|\psi \rangle \langle \psi |.}
  • إنه إسقاط ، وخاصة من الرتبة الأولى.
  • إنها متطابقة ، أيρ=ρ2.{\displaystyle \rho =\rho ^{2}.}
  • أي أن نقاوته واحدة.tr(ρ2)=1.{\displaystyle \operatorname {tr} (\rho ^{2})=1.}

من المهم التأكيد على الفرق بين مزيج احتمالي (أي مجموعة) من الحالات الكمومية وتراكب حالتين . إذا تم تحضير مجموعة بحيث يكون نصف أنظمتها في الحالة|ψ1{\displaystyle |\psi _{1}\rangle }والنصف الآخر في|ψ2{\displaystyle |\psi _{2}\rangle }ويمكن وصفها بواسطة مصفوفة الكثافة:

ρ=12(1001)،{\displaystyle \rho ={\frac {1}{2}}{\begin{pmatrix}1&0\\0&1\end{pmatrix}},}

أين|ψ1{\displaystyle |\psi _{1}\rangle }و|ψ2{\displaystyle |\psi _{2}\rangle }يُفترض أن تكون متعامدة وذات بُعد 2، تبسيطًا للأمر. من ناحية أخرى، ينتج عن التراكب الكمومي لهاتين الحالتين بسعات احتمالية متساوية الحالة النقية.|ψ=(|ψ1+|ψ2)/2،{\displaystyle |\psi \rangle =(|\psi _{1}\rangle +|\psi _{2}\rangle )/{\sqrt {2}},}مع مصفوفة الكثافة

|ψψ|=12(1111).{\displaystyle |\psi \rangle \langle \psi |={\frac {1}{2}}{\begin{pmatrix}1&1\\1&1\end{pmatrix}}.}

بخلاف الخليط الاحتمالي، يمكن لهذا التراكب أن يُظهر تداخلًا كميًا . [ 3 ] : 81

في تمثيل كرة بلوخ للكيوبت ، تمثل كل نقطة على كرة الوحدة حالة نقية. أما جميع مصفوفات الكثافة الأخرى فتمثل نقاطًا في الداخل.

هندسيًا، تُعدّ مجموعة مؤثرات الكثافة مجموعة محدبة ، والحالات النقية هي النقاط القصوى لتلك المجموعة. أبسط حالة هي حالة فضاء هيلبرت ثنائي الأبعاد، المعروف باسم الكيوبت . يمكن كتابة أي حالة مختلطة للكيوبت كمزيج خطي من مصفوفات باولي ، والتي تُشكّل مع مصفوفة الوحدة أساسًا لـ2×2{\displaystyle 2\times 2}المصفوفات ذاتية الترافق : [ 10 ] : 126

ρ=12(أنا+رxσx+رyσy+رzσz)،{\displaystyle \rho ={\frac {1}{2}}\left(I+r_{x}\sigma _{x}+r_{y}\sigma _{y}+r_{z}\sigma _{z}\right),}

حيث الأعداد الحقيقية(رx،رy،رz){\displaystyle (r_{x},r_{y},r_{z})}هي إحداثيات نقطة داخل كرة الوحدة و

σx=(0110)،σy=(0-أناأنا0)،σz=(100-1).{\displaystyle \sigma _{x}={\begin{pmatrix}0&1\\1&0\end{pmatrix}},\quad \sigma _{y}={\begin{pmatrix}0&-i\\i&0\end{pmatrix}},\quad \sigma _{z}={\begin{pmatrix}1&0\\0&-1\end{pmatrix}}.}

نقاط معرx2+رy2+رz2=1{\displaystyle r_{x}^{2}+r_{y}^{2}+r_{z}^{2}=1}تمثل الحالات النقية، بينما تمثل الحالات المختلطة بنقاط في الداخل. يُعرف هذا باسم صورة كرة بلوخ لفضاء حالة الكيوبت.

مثال: استقطاب الضوء

يُصدر المصباح المتوهج  (1) فوتونات مستقطبة عشوائية تمامًا  (2) ذات مصفوفة كثافة حالة مختلطة:
[0.5000.5]{\displaystyle {\begin{bmatrix}0.5&0\\0&0.5\end{bmatrix}}}.
بعد مرورها عبر مستقطب المستوى الرأسي  (3)، تكون جميع الفوتونات المتبقية مستقطبة رأسياً  (4) ولها مصفوفة كثافة حالة نقية:
[1000]{\displaystyle {\begin{bmatrix}1&0\\0&0\end{bmatrix}}}.

من الأمثلة على الحالات النقية والمختلطة استقطاب الضوء . يمكن وصف الفوتون الواحد بأنه ذو استقطاب دائري يميني أو يساري ، ويتم وصف ذلك بواسطة الحالات الكمومية المتعامدة.|R{\displaystyle |\mathrm {R} \rangle }و|ل{\displaystyle |\mathrm {L} \rangle }أو تراكب الاثنين: يمكن أن يكون في أي حالةα|R+β|ل{\displaystyle \alpha |\mathrm {R} \rangle +\beta |\mathrm {L} \rangle }(مع|α|2+|β|2=1{\displaystyle |\alpha |^{2}+|\beta |^{2}=1}) المقابلة للاستقطاب الخطي أو الدائري أو الإهليلجي . لنفترض الآن فوتونًا مستقطبًا رأسيًا، موصوفًا بالحالة|V=(|R+|ل)/2{\displaystyle |\mathrm {V} \rangle =(|\mathrm {R} \rangle +|\mathrm {L} \rangle )/{\sqrt {2}}}إذا مررناه عبر مستقطب دائري يسمح فقط بأحد الخيارين|R{\displaystyle |\mathrm {R} \rangle }الضوء المستقطب، أو فقط|ل{\displaystyle |\mathrm {L} \rangle }في الضوء المستقطب، يتم امتصاص نصف الفوتونات في كلتا الحالتين. قد يُوحي هذا بأن نصف الفوتونات في حالة استقطاب.|R{\displaystyle |\mathrm {R} \rangle }والنصف الآخر في الولاية|ل{\displaystyle |\mathrm {L} \rangle }لكن هذا غير صحيح: إذا مررنا(|R+|ل)/2{\displaystyle (|\mathrm {R} \rangle +|\mathrm {L} \rangle )/{\sqrt {2}}}لا يحدث امتصاص على الإطلاق عند تمرير الضوء عبر مستقطب خطي ، ولكن إذا مررنا بأي من الحالتين|R{\displaystyle |\mathrm {R} \rangle }أو|ل{\displaystyle |\mathrm {L} \rangle }يتم امتصاص نصف الفوتونات.

لا يمكن وصف الضوء غير المستقطب (مثل الضوء المنبعث من مصباح متوهج ) بأنه أي حالة من حالات الشكلα|R+β|ل{\displaystyle \alpha |\mathrm {R} \rangle +\beta |\mathrm {L} \rangle }(استقطاب خطي، أو دائري، أو بيضاوي). على عكس الضوء المستقطب، يمر الضوء غير المستقطب عبر مستقطب مع فقدان 50% من شدته بغض النظر عن اتجاه المستقطب؛ ولا يمكن استقطابه بتمريره عبر أي صفيحة موجية . مع ذلك، يمكن وصف الضوء غير المستقطب بأنه مجموعة إحصائية، على سبيل المثال، كل فوتون له إما|R{\displaystyle |\mathrm {R} \rangle }الاستقطاب أو|ل{\displaystyle |\mathrm {L} \rangle }الاستقطاب باحتمالية 1/2. سيحدث السلوك نفسه إذا كان لكل فوتون استقطاب رأسي.|V{\displaystyle |\mathrm {V} \rangle }أو الاستقطاب الأفقي|ح{\displaystyle |\mathrm {H} \rangle }باحتمالية 1/2. لا يمكن التمييز بين هاتين المجموعتين تجريبياً، ولذلك تُعتبران حالة مختلطة واحدة. في هذا المثال للضوء غير المستقطب، يساوي معامل الكثافة [ 2 ] : 75

ρ=12|RR|+12|لل|=12|حح|+12|VV|=12(1001).{\displaystyle \rho ={\frac {1}{2}}|\mathrm {R} \rangle \langle \mathrm {R} |+{\frac {1}{2}}|\mathrm {L} \rangle \langle \mathrm {L} |={\frac {1}{2}}|\mathrm {H} \rangle \langle \mathrm {H} |+{\frac {1}{2}}|\mathrm {V} \rangle \langle \mathrm {V} |={\frac {1}{2}}{\begin{pmatrix}1&0\\0&1\end{pmatrix}}.}

توجد أيضًا طرق أخرى لتوليد ضوء غير مستقطب: إحدى هذه الطرق هي إدخال عنصر عدم اليقين في تحضير الفوتون، على سبيل المثال، تمريره عبر بلورة ثنائية الانكسار ذات سطح خشن، بحيث تكتسب أجزاء مختلفة قليلاً من شعاع الضوء استقطابات مختلفة. وثمة طريقة أخرى تتمثل في استخدام الحالات المتشابكة: إذ يمكن أن يُصدر الاضمحلال الإشعاعي فوتونين يسيران في اتجاهين متعاكسين، في الحالة الكمومية.(|R،ل+|ل،R)/2{\displaystyle (|\mathrm {R} ,\mathrm {L} \rangle +|\mathrm {L} ,\mathrm {R} \rangle )/{\sqrt {2}}}تكون الحالة المشتركة للفوتونين نقية، لكن مصفوفة الكثافة لكل فوتون على حدة، والتي تُحسب بأخذ الأثر الجزئي لمصفوفة الكثافة المشتركة، تكون مختلطة تمامًا. [ 3 ] : 106

مجموعات وتنقيات مكافئة

لا يُحدد مُؤثر الكثافة المُعطى بشكلٍ فريد أي مجموعة من الحالات النقية تُنتجه؛ فعمومًا، هناك عدد لا نهائي من المجموعات المختلفة التي تُولد مصفوفة الكثافة نفسها. [ 11 ] ولا يُمكن التمييز بينها بأي قياس. [ 12 ] يُمكن توصيف المجموعات المُكافئة توصيفًا كاملًا: ليكن{صج،|ψج}{\displaystyle \{p_{j},|\psi _{j}\rangle \}}ليكن مجموعة. إذن لأي مصفوفة مركبةيو{\displaystyle U}بحيثيويو=أنا{\displaystyle U^{\dagger }U=I}( تماثل جزئي )، المجموعة{qأنا،|φأنا}{\displaystyle \{q_{i},|\varphi _{i}\rangle \}}محدد بواسطة

qأنا|φأنا=جيوأناجصج|ψج{\displaystyle {\sqrt {q_{i}}}\left|\varphi _{i}\right\rangle =\sum _{j}U_{ij}{\sqrt {p_{j}}}\left|\psi _{j}\right\rangle }

سيؤدي ذلك إلى نفس عامل الكثافة، وجميع المجموعات المتكافئة تكون من هذا الشكل.

ومن الحقائق ذات الصلة الوثيقة أن لمؤثر الكثافة المعطى عددًا لا نهائيًا من عمليات التنقية المختلفة ، وهي حالات نقية تولد مؤثر الكثافة عند أخذ أثر جزئي.

ρ=جصج|ψجψج|{\displaystyle \rho =\sum _{j}p_{j}|\psi _{j}\rangle \langle \psi _{j}|}

ليكن عامل الكثافة الناتج عن المجموعة{صج،|ψج}{\displaystyle \{p_{j},|\psi _{j}\rangle \}}، مع الولايات|ψج{\displaystyle |\psi _{j}\rangle }ليس بالضرورة متعامدًا. إذن، بالنسبة لجميع التماثلات الجزئيةيو{\displaystyle U}لدينا ذلك

|Ψ=جصج|ψجيو|أج{\displaystyle |\Psi \rangle =\sum _{j}{\sqrt {p_{j}}}|\psi _{j}\rangle U|a_{j}\rangle }

هو تنقية لـρ{\displaystyle \rho }، أين|أج{\displaystyle |a_{j}\rangle }هو أساس متعامد، وعلاوة على ذلك جميع عمليات تنقيةρ{\displaystyle \rho }وهي من هذا الشكل.

قياس

يتركأ{\displaystyle A}ليكن متغيرًا قابلًا للملاحظة للنظام، ولنفترض أن المجموعة في حالة مختلطة بحيث تكون كل حالة من الحالات النقية|ψج{\displaystyle \textstyle |\psi _{j}\rangle }يحدث باحتماليةصج{\displaystyle p_{j}}ثم يساوي عامل الكثافة المقابل

ρ=جصج|ψجψج|.{\displaystyle \rho =\sum _{j}p_{j}|\psi _{j}\rangle \langle \psi _{j}|.}

يمكن حساب القيمة المتوقعة للقياس عن طريق التوسع من حالة الحالات النقية:

أ=جصجψج|أ|ψج=جصجtr(|ψجψج|أ)=tr(جصج|ψجψج|أ)=tr(ρأ)،{\displaystyle \langle A\rangle =\sum _{j}p_{j}\langle \psi _{j}|A|\psi _{j}\rangle =\sum _{j}p_{j}\operatorname {tr} \left(|\psi _{j}\rangle \langle \psi _{j}|A\right)=\operatorname {tr} \left(\sum _{j}p_{j}|\psi _{j}\rangle \langle \psi _{j}|A\right)=\operatorname {tr} (\rho A),}

أينtr{\displaystyle \operatorname {tr} }يشير إلى الأثر . وبالتالي، فإن التعبير المألوفأ=ψ|أ|ψ{\displaystyle \langle A\rangle =\langle \psi |A|\psi \rangle }يتم استبدال الحالات النقية بـ

أ=tr(ρأ){\displaystyle \langle A\rangle =\operatorname {tr} (\rho A)}

بالنسبة للولايات المختلطة. [ 2 ] : 73

علاوة على ذلك، إذاأ{\displaystyle A}يتمتع بدقة طيفية

أ=أناأأناPأنا،{\displaystyle A=\sum _{i}a_{i}P_{i},}

أينPأنا{\displaystyle P_{i}}هو عامل الإسقاط في الفضاء الذاتي المقابل للقيمة الذاتيةأأنا{\displaystyle a_{i}}، يتم إعطاء عامل الكثافة بعد القياس بواسطة [ 13 ] [ 14 ]

ρأنا=PأناρPأناtr[ρPأنا]{\displaystyle \rho _{i}'={\frac {P_{i}\rho P_{i}}{\operatorname {tr} \left[\rho P_{i}\right]}}}

عند الحصول على النتيجة i . في حالة عدم معرفة نتيجة القياس، يتم وصف المجموعة بدلاً من ذلك بواسطة

ρ=أناPأناρPأنا.{\displaystyle \;\rho '=\sum _{i}P_{i}\rho P_{i}.}

إذا افترضنا أن احتمالات نتائج القياس هي دوال خطية للمسقطاتPأنا{\displaystyle P_{i}}إذاً، يجب أن تُعطى هذه القيم بأثر المسقط مع مؤثر الكثافة. تُبين نظرية غليسون أنه في فضاءات هيلبرت ذات البعد 3 أو أكبر، يمكن استبدال افتراض الخطية بافتراض عدم السياقية . [ 15 ] يمكن إزالة هذا القيد على البعد بافتراض عدم السياقية لـ POVMs أيضًا، [ 16 ] [ 17 ] ولكن هذا الأمر وُوجه بانتقادات لعدم وجود أساس فيزيائي له. [ 18 ]

إنتروبيا

إنتروبيا فون نيومانS{\displaystyle S}يمكن التعبير عن خصائص الخليط بدلالة القيم الذاتية لـρ{\displaystyle \rho }أو بدلالة أثر ولوغاريتم عامل الكثافةρ{\displaystyle \rho }. منذρ{\displaystyle \rho }هو مؤثر شبه موجب، وله تحليل طيفي بحيثρ=أناλأنا|φأناφأنا|{\displaystyle \rho =\textstyle \sum _{i}\lambda _{i}|\varphi _{i}\rangle \langle \varphi _{i}|}، أين|φأنا{\displaystyle |\varphi _{i}\rangle }هي متجهات متعامدة،λأنا0{\displaystyle \lambda _{i}\geq 0}، وλأنا=1{\displaystyle \textstyle \sum \lambda _{i}=1}ثم إنتروبيا نظام كمومي ذي مصفوفة كثافةρ{\displaystyle \rho }يكون

S=-أناλأناlnλأنا=-tr(ρlnρ).{\displaystyle S=-\sum _{i}\lambda _{i}\ln \lambda _{i}=-\operatorname {tr} (\rho \ln \rho ).}

يشير هذا التعريف إلى أن إنتروبيا فون نيومان لأي حالة نقية تساوي صفرًا. [ 19 ] : 217 إذاρأنا{\displaystyle \rho _{i}}إذا كانت الحالات التي لها دعم على الفضاءات الفرعية المتعامدة، فإن إنتروبيا فون نيومان لمجموعة محدبة من هذه الحالات،

ρ=أناصأناρأنا،{\displaystyle \rho =\sum _{i}p_{i}\rho _{i},}

يتم تحديدها بواسطة إنتروبيا فون نيومان للحالاتρأنا{\displaystyle \rho _{i}}وإنتروبيا شانون لتوزيع الاحتماليةصأنا{\displaystyle p_{i}}:

S(ρ)=ح(صأنا)+أناصأناS(ρأنا).{\displaystyle S(\rho )=H(p_{i})+\sum _{i}p_{i}S(\rho _{i}).}

عندما تكون الولاياتρأنا{\displaystyle \rho _{i}}إذا لم يكن للمصفوفات دعامات متعامدة، فإن المجموع على الجانب الأيمن يكون أكبر تمامًا من إنتروبيا فون نيومان للتركيبة المحدبة.ρ{\displaystyle \rho }[ 3 ] : 518

بالنظر إلى عامل الكثافةρ{\displaystyle \rho }وقياس إسقاطي كما في القسم السابق، الدولةρ{\displaystyle \rho '}مُعرَّف بواسطة التركيبة المحدبة

ρ=أناPأناρPأنا،{\displaystyle \rho '=\sum _{i}P_{i}\rho P_{i},}

والتي يمكن تفسيرها على أنها الحالة الناتجة عن إجراء القياس دون تسجيل النتيجة التي حدثت، [ 10 ] : 159 لها إنتروبيا فون نيومان أكبر من إنتروبيا فون نيومان لـρ{\displaystyle \rho }إلا إذاρ=ρ{\displaystyle \rho =\rho '}ومع ذلك، من الممكن أنρ{\displaystyle \rho '}ينتج عن القياس المعمم ، أو POVM ، أن يكون له إنتروبيا فون نيومان أقل منρ{\displaystyle \rho }[ 3 ] : 514

معادلة فون نيومان للتطور الزمني

كما تصف معادلة شرودنغر كيفية تطور الحالات النقية مع الزمن، تصف معادلة فون نيومان (المعروفة أيضًا باسم معادلة ليوفيل-فون نيومان ) كيفية تطور مؤثر الكثافة مع الزمن. وتنص معادلة فون نيومان على ما يلي [ 20 ] [ 21 ] [ 22 ]

أناددتρ=[ح،ρ] ،{\displaystyle i\hbar {\frac {d}{dt}}\rho =[H,\rho ]~,}

حيث تشير الأقواس إلى المبدل .

لا تصح هذه المعادلة إلا عندما يتم اعتبار عامل الكثافة في صورة شرودنغر ، على الرغم من أن هذه المعادلة تبدو للوهلة الأولى وكأنها تحاكي معادلة هايزنبرغ للحركة في صورة هايزنبرغ ، مع اختلاف جوهري في الإشارة:

أناددتأح=-[حح،أح] ،{\displaystyle i\hbar {\frac {d}{dt}}A_{\text{H}}=-[H_{\text{H}},A_{\text{H}}]~,}

أينأح(ت){\displaystyle A_{\text{H}}(t)}هو أحد عوامل صورة هايزنبرغ ؛ ولكن في هذه الصورة، لا تعتمد مصفوفة الكثافة على الزمن ، وتضمن الإشارة النسبية أن المشتق الزمني للقيمة المتوقعةأ{\displaystyle \langle A\rangle }[ 7 ] النتيجة هي نفسها كما في صورة شرودنغر .

إذا كان الهاميلتوني مستقلاً عن الزمن، فيمكن حل معادلة فون نيومان بسهولة للحصول على

ρ(ت)=هـ-أناحت/ρ(0)هـأناحت/.{\displaystyle \rho (t)=e^{-iHt/\hbar }\rho (0)e^{iHt/\hbar }.}

بالنسبة لهاملتونيان أكثر عمومية، إذاجي(ت){\displaystyle G(t)}إذا كان لدينا مُوَصِّل دالة الموجة على فترة زمنية معينة، فإن التطور الزمني لمصفوفة الكثافة على نفس الفترة الزمنية يُعطى بالعلاقة التالية:

ρ(ت)=جي(ت)ρ(0)جي(ت).{\displaystyle \rho (t)=G(t)\rho (0)G(t)^{\dagger }.}

إذا دخل المرء في صورة التفاعل ، يختار التركيز على أحد المكوناتح1{\displaystyle H_{1}}من الهاميلتونيح=ح0+ح1{\displaystyle H=H_{0}+H_{1}}، معادلة تطور عامل كثافة صورة التفاعلρأنا(ت){\displaystyle \rho _{\,\mathrm {I} }(t)}يمتلك بنية مطابقة لمعادلة فون نيومان، باستثناء أنه يجب أيضًا تحويل الهاميلتوني إلى الصورة الجديدة:

أناددتρأنا(ت)=[ح1،أنا(ت)،ρأنا(ت)]،{\displaystyle {\displaystyle i\hbar {\frac {d}{dt}}\rho _{\text{I}}(t)=[H_{1,{\text{I}}}(t),\rho _{\text{I}}(t)],}}

أينح1،أنا(ت)=هـأناح0ت/ح1هـ-أناح0ت/{\displaystyle {\displaystyle H_{1,{\text{I}}}(t)=e^{iH_{0}t/\hbar }H_{1}e^{-iH_{0}t/\hbar }}}.

دوال ويغنر والتشابهات الكلاسيكية

يمكن أيضًا تمثيل مؤثر مصفوفة الكثافة في فضاء الطور . في ظل تحويل ويغنر ، تتحول مصفوفة الكثافة إلى دالة ويغنر المكافئة .

دبليو(x،ص) =دهـو 1π-ψ*(x+y)ψ(x-y)هـ2أناصy/دy.{\displaystyle W(x,p)\,\ {\stackrel {\mathrm {def} }{=}}\ \,{\frac {1}{\pi \hbar }}\int _{-\infty }^{\infty }\psi ^{*}(x+y)\psi (x-y)e^{2ipy/\hbar }\,dy.}

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

دبليو(x،ص،ت)ت=-{{دبليو(x،ص،ت)،ح(x،ص)}}،{\displaystyle {\frac {\partial W(x,p,t)}{\partial t}}=-\{\{W(x,p,t),H(x,p)\}\},}

أينح(x،ص){\displaystyle H(x,p)}هو الهاميلتوني، و{{،}}{\displaystyle \{\{\cdot ,\cdot \}\}}هو قوس مويال ، وهو تحويل المبدل الكمي .

تكون معادلة تطور دالة ويغنر مماثلة لمعادلة حدها الكلاسيكي، وهي معادلة ليوفيل في الفيزياء الكلاسيكية . في حالة انعدام ثابت بلانك{\displaystyle \hbar }،دبليو(x،ص،ت){\displaystyle W(x,p,t)}يختزل إلى دالة كثافة الاحتمال الكلاسيكية لليوفيل في فضاء الطور .

تطبيقات نموذجية

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

  • تستخدم الميكانيكا الإحصائية مصفوفات الكثافة، وخاصة للتعبير عن فكرة تحضير النظام عند درجة حرارة غير صفرية. ويؤدي بناء مصفوفة الكثافة باستخدام مجموعة متعارف عليها إلى نتيجة على الشكل التالي:ρ=خبرة(-βح)/Z(β){\displaystyle \rho =\exp(-\beta H)/Z(\beta )}، أينβ{\displaystyle \beta }هي درجة الحرارة المعكوسة(كبتي)-1{\displaystyle (k_{\rm {B}}T)^{-1}}وح{\displaystyle H}هو هاميلتونيان النظام. شرط التطبيع الذي يكون فيه أثرρ{\displaystyle \rho }إن كون القيمة تساوي 1 يحدد دالة التقسيم على النحو التالي:Z(β)=ترخبرة(-βح){\displaystyle Z(\beta )=\mathrm {tr} \exp(-\beta H)}إذا لم يكن عدد الجسيمات المشاركة في النظام مؤكدًا، فيمكن تطبيق مجموعة إحصائية كبرى ، حيث تُستمد الحالات التي جُمعت لتكوين مصفوفة الكثافة من فضاء فوك . [ 23 ] : 174
  • تتضمن نظرية التفكك الكمومي عادةً أنظمة كمومية غير معزولة تُطوّر تشابكًا مع أنظمة أخرى، بما في ذلك أجهزة القياس. تُسهّل مصفوفات الكثافة وصف العملية وحساب نتائجها. يُفسّر التفكك الكمومي سبب انتقال نظام يتفاعل مع بيئة من حالة نقية، تُظهر تراكبات، إلى حالة مختلطة، وهي مزيج غير متماسك من البدائل الكلاسيكية. هذا الانتقال قابل للانعكاس جوهريًا، لأن الحالة المُدمجة للنظام والبيئة تظل نقية، ولكنه عمليًا غير قابل للانعكاس، لأن البيئة نظام كمومي كبير ومعقد للغاية، وليس من الممكن عكس تفاعلهما. لذا، يُعدّ التفكك الكمومي بالغ الأهمية لتفسير الحد الكلاسيكي لميكانيكا الكم، ولكنه لا يُفسّر انهيار الدالة الموجية، لأن جميع البدائل الكلاسيكية لا تزال موجودة في الحالة المختلطة، ولا يختار انهيار الدالة الموجية سوى واحد منها. [ 24 ]
  • وبالمثل، في الحوسبة الكمومية ، ونظرية المعلومات الكمومية ، والأنظمة الكمومية المفتوحة ، وغيرها من المجالات التي تتسم فيها عملية تحضير الحالة بالتشويش وإمكانية حدوث فقدان الترابط الكمومي، تُستخدم مصفوفات الكثافة بشكل متكرر. غالبًا ما يتم نمذجة التشويش عبر قناة إزالة الاستقطاب أو قناة تخميد السعة . التصوير المقطعي الكمومي هو عملية يتم من خلالها، بالنظر إلى مجموعة من البيانات التي تمثل نتائج القياسات الكمومية، حساب مصفوفة كثافة تتوافق مع نتائج تلك القياسات. [ 25 ] [ 26 ]
  • عند تحليل نظام يحتوي على العديد من الإلكترونات، مثل الذرة أو الجزيء ، فإن التقريب الأولي غير الكامل ولكنه مفيد هو اعتبار الإلكترونات غير مترابطة أو أن لكل منها دالة موجية مستقلة. هذه هي نقطة البداية المعتادة عند بناء محدد سلاتر في طريقة هارتري-فوك .شمال{\displaystyle N}الإلكترونات تملأشمال{\displaystyle N}دوال الموجة للجسيم الواحد|ψأنا{\displaystyle |\psi _{i}\rangle }وإذا اقتصر الأمر على رصد الجسيمات المفردة فقط، فإن قيمها المتوقعة لـشمال{\displaystyle N}يمكن حساب نظام الإلكترونات باستخدام مصفوفة الكثافةأنا=1شمال|ψأناψأنا|{\textstyle \sum _{i=1}^{N}|\psi _{i}\rangle \langle \psi _{i}|}( مصفوفة كثافة الجسيمات المفردة لـشمال{\displaystyle N}نظام الإلكترونات). [ 27 ]

صياغة الحالات باستخدام الجبر C*

من المقبول عمومًا الآن أن وصف ميكانيكا الكم الذي تُمثل فيه جميع المؤثرات ذاتية الترافق الكميات القابلة للرصد غير مقبول. [ 28 ] [ 29 ] لهذا السبب، تُعرَّف الكميات القابلة للرصد بعناصر جبر C* مجرد A (أي جبر بدون تمثيل مميز كجبر للمؤثرات)، وتُعرَّف الحالات بأنها دوال خطية موجبة على A. مع ذلك، باستخدام بناء GNS ، يمكننا استعادة فضاءات هيلبرت التي تُجسد A كجبر فرعي للمؤثرات.

هندسياً، الحالة النقية على جبر C* -A هي حالة تمثل نقطة قصوى لمجموعة جميع الحالات على A. وبحسب خصائص بناء GNS، فإن هذه الحالات تتوافق مع التمثيلات غير القابلة للاختزال لـ A.

تتوافق حالات الجبر C* للمؤثرات المدمجة K ( H ) تمامًا مع مؤثرات الكثافة، وبالتالي فإن الحالات النقية لـ K ( H ) هي بالضبط الحالات النقية بمعنى ميكانيكا الكم.

يمكن ملاحظة أن الصيغة الجبرية C* تشمل الأنظمة الكلاسيكية والكمومية على حد سواء. عندما يكون النظام كلاسيكيًا، يصبح جبر الكميات القابلة للرصد جبرًا تبديليًا C*. في هذه الحالة، تصبح الحالات مقاييس احتمالية.

تاريخ

طُوِّر هذا الشكل الرسمي للمؤثرات والمصفوفات عام 1927 على يد جون فون نيومان [ 30 ] ، وبشكل مستقل، ولكن بشكل أقل منهجية، على يد ليف لانداو [ 31 ] ، ولاحقًا عام 1946 على يد فيليكس بلوخ [ 32 ] . قدّم فون نيومان المصفوفة لتطوير كلٍّ من الميكانيكا الإحصائية الكمومية ونظرية القياسات الكمومية. أُدخل مصطلح الكثافة على يد ديراك عام 1931 عندما استخدم مؤثر فون نيومان لحساب سحب كثافة الإلكترونات [ 33 ] [ 34 ] .

في الوقت الحاضر، اكتسب مصطلح "مصفوفة الكثافة" أهمية خاصة به، ويتوافق مع مقياس احتمالية فضاء الطور الكلاسيكي (توزيع احتمالية الموضع والزخم) في الميكانيكا الإحصائية الكلاسيكية ، والذي قدمه يوجين ويجنر في عام 1932. [ 5 ]

في المقابل، كان الدافع الذي ألهم لاندو هو استحالة وصف نظام فرعي من نظام كمي مركب بواسطة متجه حالة. [ 31 ]

انظر أيضاً

ملاحظات ومراجع

  1. شانكار، رامامورتي (2014). مبادئ ميكانيكا الكم (الطبعة الثانية، [الطبعة التاسعة عشرة المصححة]  ). نيويورك، نيويورك: سبرينغر. ISBN 978-0-306-44790-7.
  2. 1 2 3 4 بيريز، آشر (1995). نظرية الكم: المفاهيم والأساليب . كلوير. ISBN 978-0-7923-3632-7. OCLC 901395752 . 
  3. 1 2 3 4 5 6 7 8 نيلسن، مايكل؛ تشوانغ، إسحاق (2000)، الحوسبة الكمومية والمعلومات الكمومية ، مطبعة جامعة كامبريدج ، ISBN 978-0-521-63503-5.
  4. بالنتين، ليزلي (2009). "مصفوفة الكثافة". موسوعة الفيزياء الكمية . برلين، هايدلبرغ: سبرينغر برلين هايدلبرغ. ص 166. doi : 10.1007/978-3-540-70626-7_51 . ISBN  978-3-540-70622-9.
  5. 1 2 فانو، يو. (1957). "وصف الحالات في ميكانيكا الكم باستخدام مصفوفة الكثافة وتقنيات المؤثرات". مراجعات الفيزياء الحديثة . 29 (1): 74-93 . Bibcode : 1957RvMP...29...74F . doi : 10.1103/RevModPhys.29.74 .
  6. هوليفو، ألكسندر س. (2001). البنية الإحصائية لنظرية الكم . سلسلة محاضرات في الفيزياء. سبرينغر. ISBN 3-540-42082-7. OCLC 318268606 . 
  7. 1 2 3 هول، برايان سي. (2013). "الأنظمة والأنظمة الفرعية، الجسيمات المتعددة". نظرية الكم للرياضيين . نصوص الدراسات العليا في الرياضيات. المجلد 267. الصفحات 419-440 . doi : 10.1007/978-1-4614-7116-5_19 . ISBN   978-1-4614-7115-8.
  8. ^ كوهين تنودجي، كلود. ديو، برنارد. لالوي ، فرانك (2019). ميكانيكا الكم، المجلد الأول . فاينهايم، ألمانيا: جون وايلي وأولاده. ص 301 – 303. ISBN  978-3-527-34553-3..
  9. ديفيدسون، إرنست روي (1976). مصفوفات الكثافة المختزلة في الكيمياء الكمية . دار النشر الأكاديمية ، لندن.
  10. 1 2 وايلد، مارك م. (2017). نظرية المعلومات الكمومية ( الطبعة الثانية). مطبعة جامعة كامبريدج. arXiv : 1106.1445 . doi : 10.1017/9781316809976.001 . ISBN  978-1-107-17616-4. OCLC 973404322 . S2CID 2515538 .  
  11. كيركباتريك، ك. أ. (فبراير 2006). "نظرية شرودنغر-هـ. ج. و.". رسائل أسس الفيزياء . 19 (1): 95-102 . arXiv : quant-ph/0305068 . Bibcode : 2006FoPhL..19...95K . doi : 10.1007/s10702-006-1852-1 . ISSN: 0894-9875 . S2CID : 15995449 .  
  12. أوكس، فيلهلم (1981-11-01). "بعض التعليقات على مفهوم الحالة في ميكانيكا الكم". إركنتنيس . 16 (3): 339-356 . doi : 10.1007/BF00211375 . ISSN 1572-8420 . S2CID 119980948 .  
  13. ^ لودرز غيرهارت (1950). "Über die Zustandsänderung durch den Messprozeß". أنالين دير فيزيك . 443 ( 5– 8): 322. بيب كود : 1950AnP...443..322L . دوى : 10.1002/andp.19504430510 .ترجمة ك. أ. كيركباتريك ، لودرز، جيرهارت (2006-04-03). "حول تغير الحالة الناتج عن عملية القياس". حوليات الفيزياء . 15 (9): 663-670 . arXiv : quant-ph/0403007 . Bibcode : 2006AnP...518..663L . doi : 10.1002/andp.200610207 . S2CID 119103479 . 
  14. ^ بوش، بول ؛ لاهتي، بيكا (2009)، “قاعدة لودرز”، في دانيال غرينبرغر؛ هنتشل، كلاوس؛ وينرت ، فريدل (محرران)، خلاصة وافية لفيزياء الكم ، سبرينغر برلين هايدلبرغ، الصفحات من 356 إلى 358، دوى : 10.1007 / 978-3-540-70626-7_110 ، ISBN  978-3-540-70622-9
  15. غليسون، أندرو م. (1957). "القياسات على الفضاءات الجزئية المغلقة لفضاء هيلبرت" . مجلة الرياضيات بجامعة إنديانا . 6 (4): 885-893 . doi : 10.1512/iumj.1957.6.56050 . MR 0096113 . 
  16. بوش، بول (2003). "الحالات الكمومية والكميات القابلة للرصد المعممة: برهان بسيط لنظرية غليسون". رسائل المراجعة الفيزيائية . 91 (12) 120403. arXiv : quant-ph/9909073 . Bibcode : 2003PhRvL..91l0403B . doi : 10.1103 /PhysRevLett.91.120403 . PMID 14525351. S2CID 2168715 .  
  17. كيفز، كارلتون م .؛ فوكس، كريستوفر أ.؛ مان، كيران ك.؛ رينيس، جوزيف م. (2004). "اشتقاقات من نوع غليسون لقاعدة الاحتمال الكمي للقياسات المعممة". أسس الفيزياء . 34 (2): 193-209 . arXiv : quant-ph/0306179 . Bibcode : 2004FoPh...34..193C . doi : 10.1023/B:FOOP.0000019581.00318.a5 . S2CID 18132256 . 
  18. أندريه غرودكا؛ باويل كورزينسكي (2008). "هل يوجد سياق لكيوبت واحد؟". رسائل المراجعة الفيزيائية . 100 (16) 160401. arXiv : 0705.0181 . Bibcode : 2008PhRvL.100p0401G . doi : 10.1103/PhysRevLett.100.160401 . PMID 18518167. S2CID 13251108 .  
  19. ريفيل، إليانور ج .؛ بولاك، وولفغانغ هـ. (4 مارس 2011). الحوسبة الكمومية: مقدمة مبسطة . مطبعة معهد ماساتشوستس للتكنولوجيا. ISBN 978-0-262-01506-6.
  20. بروير، هاينز؛ بيتروتشيوني، فرانشيسكو (2002)، نظرية الأنظمة الكمومية المفتوحة ، مطبعة جامعة أكسفورد، ص 110، ISBN  978-0-19-852063-4
  21. شوابل، فرانز (2002)، الميكانيكا الإحصائية ، سبرينغر، ص 16، ISBN  978-3-540-43163-3
  22. مولر-كيرستن، هارالد جيه دبليو (2008)، الميكانيكا الكلاسيكية والنسبية ، وورلد ساينتيفيك، ص 175-179 ، ISBN  978-981-283-251-1
  23. كاردار، مهران (2007). الفيزياء الإحصائية للجسيمات . مطبعة جامعة كامبريدج . ISBN 978-0-521-87342-0. OCLC 860391091 . 
  24. شلوسهاور، م. (2019). "فك الترابط الكمي". تقارير الفيزياء . 831 : 1-57 . arXiv : 1911.06282 . Bibcode : 2019PhR...831....1S . doi : 10.1016/j.physrep.2019.10.001 . S2CID 208006050 . 
  25. غراناد، كريستوفر؛ كومبس، جوشوا؛ كوري، دي جي (2016-01-01). "التصوير المقطعي البايزي العملي". المجلة الجديدة للفيزياء . 18 (3) 033024. arXiv : 1509.03770 . Bibcode : 2016NJPh...18c3024G . doi : 10.1088/1367-2630/18/3/033024 . ISSN 1367-2630 . S2CID 88521187 .  
  26. أرديلا، لويس؛ هايل، ماركوس؛ إيكاردت، أندريه (28 ديسمبر 2018). "قياس مصفوفة كثافة الجسيمات المفردة للفيرميونات والبوزونات ذات النواة الصلبة في شبكة بصرية". رسائل المراجعة الفيزيائية . 121 (260401): 6. arXiv : 1806.08171 . Bibcode : 2018PhRvL.121z0401P . doi : 10.1103/PhysRevLett.121.260401 . PMID 30636128. S2CID 51684413 .  
  27. كيتل، تشارلز (1963). نظرية الكم للمواد الصلبة . نيويورك: وايلي. ص 101. 
  28. انظر الملحق، ماكي، جورج وايتلو (1963)، الأسس الرياضية لميكانيكا الكم ، كتب دوفر في الرياضيات، نيويورك: منشورات دوفر ، رقم ISBN 978-0-486-43517-6{{citation}}عدم توافق رقم ISBN / التاريخ ( مساعدة )
  29. إمخ، جيرارد ج. (1972)، الأساليب الجبرية في الميكانيكا الإحصائية ونظرية الحقل الكمومي ، وايلي-إنترساينس ، رقم ISBN 978-0-471-23900-0
  30. ^ فون نيومان ، جون (1927) ، “Wahrscheinlichkeitstheoretischer Aufbau der Quantenmechanik” ، Göttinger Nachrichten ، 1 : 245–272
  31. 1 2 "مشكلة التخميد في ميكانيكا الموجات (1927)". أوراق بحثية مجمعة لـ ل. د. لانداو . 1965. ص 8-18 . doi : 10.1016/B978-0-08-010586-4.50007-9 . ISBN  978-0-08-010586-4.
  32. فانو، أوجو (1995). "مصفوفات الكثافة كمتجهات استقطاب". رينديكونتي لينسي . 6 (2): 123-130 . doi : 10.1007/BF03001661 . S2CID 128081459 . 
  33. ديراك، PAM (يوليو 1930). "ملاحظة حول ظواهر التبادل في ذرة توماس" . وقائع الجمعية الفلسفية في كامبريدج . 26 (3): 376-385 . Bibcode : 1930PCPS...26..376D . doi : 10.1017/S0305004100016108 . ISSN 0305-0041 . 
  34. ديراك، PAM (أبريل 1931). "ملاحظة حول تفسير مصفوفة الكثافة في مسألة الإلكترونات المتعددة" . وقائع الجمعية الفلسفية في كامبريدج . 27 (2): 240-243 . Bibcode : 1931PCPS...27..240D . doi : 10.1017/S0305004100010343 . ISSN 0305-0041 .