|
20080511～
13と7と11の倍数の論理積は13と7と11の積の倍数である。
和ァ・・・
さっきな、静かな図書館で落ち着いて手計算やった結果を家に帰ってscilabにぶち込んで確認してたんだ。5次行列をがんばって解いてた。
ついうっかり間違えて P'-inv(P) と打つところを P-inv(P) と打ってみたら、両方とも行列の中身がゼロになって、偶然気づいたんだけど 角運動量の固有状態になるユニタリ行列のx軸成分って ユニタリでもありながらエルミートでもあったのな！ なんかn次の行列全部に言えそうな気がする。 つまりn状態系すべてのx軸の固有状態がな。 パウリ行列σxを対角化するためのアレ （|sx+>、|sx->） だけじゃなかったのか・・・！ まじでびっくりした。 （|sx1>、|sx2>、・・・、|sxn-1>、|sxn>、）←こいつがユニタリ＆エルミート nが奇数だけじゃなく、偶数でも行けそう。ぼちぼち確かめてみるわ。 もしかしたらy成分はy成分で、ユニタリかつ歪エルミートかもしれん そういや行列式がiになる感じがしたっけ 並べ方次第かなあ。ノルムを維持しつつ符号を並べ替えたら、一意ではない気がする なんか特徴ありそう x方向はユニタリかつエルミートときどき歪エルミート？？ 5次ユニタリの固有値方程式を解いてみたい。 5次とは言っても簡単に解ける形になるはずだ。 代数的にではなくとも、三角関数を使えば解けるような解の公式みたいな条件・・・
PR
ムシャクシャしてやった。
パラメータや定数として文字を定義するのはまだ諦めていない 5状態・7状態・多状態系の角運動量演算子の行列表現 a=poly(0,"a") b=poly(0,"b") a+b !--error 144 指定したオペランドに関する処理が定義されていません. 関数 %p_f_p を確認またはオーバーロード定義してください. な、なんだってー！？ ウルフラム先生は、ちょーっと複雑にしただけで金要求してくるし 先日の続きで、5状態系に拡張してみた結果 こう定義してしまうと、球面調和関数Lx,Ly,Lzの固有値が揃わないどころか、整数から外れるものが出てしまうので、前提が間違っていたようだ。(背理法) ということで、以下のように定義してみると Lzを算出するだけで Lzの対角要素が-2,-1,0,1,2(固有値)であることを利用して B^2=2*2 A^2-B^2=1*2 の連立方程式が出来上がり、解くことができるようになった。 同様に、7状態系でも かなり楽ができて 対角要素だけ3つ計算すれば、あとはゼロか符号反転なので計算しなくてもよく C^2=3*2 B^2-C^2=2*2 A^2-B^2=1*2 の連立方程式ができ、これを解けば、係数は定まる。 9状態系以降も、同様に計算可能。 ただ、パラメータが増えるので、A,B,C・・・などというよりはa1,a2,a3・・・などとしたほうがよいだろう
昨日の日記はこちらです。3状態系の量子力学、行列表現です
ユニタリ行列Pの固有値は、僕自身がネタバレを回避するために、まだやってません 手計算かなにかで3次方程式をまず解いてみたいのです。 うろ覚えですが、確かこの量子数の個数って奇数でしたよね。 3次の次が4次すっ飛ばして5次方程式なのはなんだか残念です。 物理法則をやや無視して、無理やり4次のを作って解くのもありでしょうかねえ 2次のパウリ行列でできたんだから4次もできるんじゃ？と期待してしまいます y 軸方向の球面調和関数 A=[0,-1,0;1,0,-1;0,1,0] A=A*%i/sqrt(2) p1=[1/2;%i/sqrt(2);-1/2] p2=[1/sqrt(2);0;1/sqrt(2)] p3=[1/2;-%i/sqrt(2);-1/2] P=[p1,p2,p3] det(P) P'*A*P x 方向 A=[0,1,0;1,0,1;0,1,0] A=A/sqrt(2) p1=[1/2;1/sqrt(2);1/2] p2=[1/sqrt(2);0;-1/sqrt(2)] p3=[1/2;-1/sqrt(2);1/2] P=[p1,p2,p3] det(P) P'*A*P ああそういえば、スピン3/2で思い出したんですが 3/2階やほかの半整数階のテンソルってあるんでしょうか。
どこからどこまで自分で再発見した車輪なのかおぼろげなのですが、放送大学「量子物理」の小形正男さんによる14話、
3状態系の調和関数の代数的表現(行列表現)を参考に、ある程度自力で2状態系に応用したら それがスピンやパウリ行列そのものだった といったような記憶があって 先日の日記で、余裕がなかったので定性的な話を多めにしたら 自分でしゃべってようやく自分で気づきました。 固有ベクトルがユニタリになるのは必然なんじゃないかなって。 こういうところに定性的な話の重要性があるようですね。 そこで、量子物理14話に今一度戻って、3状態系の球面調和関数を例に、固有状態を計算してみたのがこちらです。 行列表現は抽象化されていて、あまり具体例を見ない気がするので、よかったら参考にしてもらえればと。 まず、このように3状態を定義します。 一番高い状態 中くらいの状態 一番低い状態 これらの状態を1つ上に押しやるのが上昇演算子L+ 1つ下に下げるのが下降演算子L-とすると、 3状態横に並べて、まとめてこのような行列の積で書けるはずです。 ただし、一番上の状態をさらに上げようとしても上がりませんし 一番下の状態をさらに下げようとしても下げられません。 この式の数学的な意味は、単に単位行列Eに行列Lを掛け算したら、行列Lになるよってそれだけのことです。 L=L×E ここで、係数Aはまだ確定していないものとします。 「量子物理」14話では、理由はよくわからないのですが 上昇・下降演算子L+とL-は角運動量のxとy軸方向成分LxとLyを使って、 L+=Lx+iLy L-=Lx-iLy と定義されており、また角運動量のz軸方向成分Lzは、交換関係を用いて iLz=[Lx,Ly] と定義されていたので、それに従うと、以下のようになります。 ここで、Lx、Ly、Lzの固有値を求めてみましょう。 Lx,Lyの固有値はλ=0,±A/√2 Lzの固有値はλ=0、±A^2/2 になるので、3軸の固有値がすべて整数、特に0、±1に揃うようなAを探すと、A=√2となります。 したがって、Lx,Ly,Lzはそれぞれ であることが判明しました。 固有値が求められたので、今度は固有ベクトルを求めてみましょう。 固有値が重解ではないので、固有ベクトルを規格化することが可能です。 Lxの固有ベクトルは、固有値が1,0,-1の場合にそれぞれ になるため、列ベクトルのノルムが1となるように規格化すると このようになります。 この3つ並んだ列ベクトルを行列とみなすと このようになって、行列式を計算してみると-1となり、ユニタリ行列であることがわかるかと思います。 LzはLz自身がすでに対角化されているので、状態固有ベクトルは単位行列となり、行列式はもちろん1の、ユニタリ行列です。 Lyの状態固有ベクトルは・・・材料は揃っているので計算してみてください。 昨日の日記のようになるはずです。 余裕があったら、3軸それぞれの対角化までやってみてください。
y方向角運動量の、行列表現のやつです。このエルミート行列Aを対角化する際は、対角化のための行列Pが規格化できるため、きっちりユニタリ行列になってくれます。Pの行列式の値はiです
バレバレだよ！お前らが犯してなかったらこんなに宇宙豊かじゃないよってくらい恩恵被ってるよ！！！
忘れてましたが、テーマを移る前に、できることから片づけてしまいましょう。
スピノルが何者なのかはまだわからないとして、 先日、ディラック方程式をタキオン用にちょっと改造したやつがあってですね クラインゴルドン方程式がターディオン用だったら こうなるところを タキオン用に こうしたいので、 ディラック方程式のベータだけ 結局、ターディオンだとこうなるところを タキオン用にこうしようかなと こう、4つ目の行列だけ変更した、連立偏微分方程式にします。 この解であるΨにx,y,z方向に進行する平面波 を仮定すると それぞれ、 x方向 y方向 z方向 このような固有値問題になります。 このうち、z方向の固有値問題だけは、2次のペアに分解できますが、 x,y方向の固有値問題は4次のまま分解できません。 にもかかわらずエネルギー固有値は これを解いた2つしかないため、重複するかと思ってジョルダン標準形を身構えてしまったのですが そこはやはり量子力学、ジョルダン標準形が出てくるシーンがなかなかシャイなようで めったにお目にかかれないみたいです。 まだ僕はジョルダン標準形に慣れていないので、wolfram alphaでネタバレを見てしまったのですが 固有値の片割れE0を それぞれ以下のようなブラケットに囲まれて対角化されるようです。 ｘ方向 ｙ方向 ｚ方向 エネルギー固有値が、タキオン用に改造したクラインゴルドン方程式と一致するのはまあ、数学的に無矛盾なだけの必然なんでしょうたぶん 残念ながら今の僕の知識ではまだ、対角化された行列を囲むブラとケットに相当するスピノル が何を意味するのかがわかりません。 ターディオンにおけるp=0を非相対論的極限、光速付近の状態を超相対論的極限(僕は常相対論的極限と呼びたい)と呼ぶように タキオンにおけるE=0を超相対論的極限、光速付近の状態を常相対論的極限と呼べば、 タキオンとターディオンの常相対論的極限は、ルクソンの状態に統一的に近似されるはずだと思います。 また、「タキオンのほうの」超相対論的極限E=0にも何かが見いだせるのではないかと期待してます あとできることといえばアレですね、カイラリティだかヘリシティだかいうやつ たぶんテンソルほど難しくないと思いますので、近日中に仕上げたいと思います。 できればターディオンの非相対論的極限もあげておきたい。っていうか基本そっちのが先なんだけど。 追記 ２１：１１ 間違えました訂正 ただしくはこうです ｘ方向 ｙ方向 ｚ方向
特に、「その他の状態」の項目が素晴らしいです！
から が導き出せるよね って言ってるだけなんですが、これがかなり奥深くて 興味のある人にとっては、この角度は何に対する何の角度なんだ！？って疑問がふつう湧くじゃないですか！ で、その計算のための道具が、上に全部まとめられてるんですよ！ とりあえずθを、θxとθyとθzの3タイプ用意しましょう。 ケットベクトルはこのたった1つの例をもとに、このように組み立てられるはずです。 そして、ケットに対するブラも、ただエルミート共役を取るだけなので、以下のようにできるはずです。 そうして、期待値を求めますと <θx|sx|θx>、<θx|sy|θx>、<θx|sz|θx> <θy|sx|θy>、<θy|sy|θy>、<θy|sz|θy> <θz|sx|θz>、<θz|sy|θz>、<θz|sz|θz> の9通りが考えられるわけですが。興味のある人はぜひ計算してみてください。 sx、sy、szの定義もちゃんとwikiにそろってます。 そうすると、 となるはずで、これが何を意味しているかというと、こういうことです。 軸からどのくらい傾いているのかの角度によるコサイン成分を知りたければ、軸と同じsをθで挟めばいい たとえばcosθx=<θx|sx|θx>、cosθy=<θy|sy|θy>、cosθz=<θz|sz|θz>のような感じです。 ただし、サイン成分を必ずしも算出できる組み合わせとは限らない といったところが、z軸＞x軸＞y軸のような不平等さの表れとしてでているみたいです。 まあ仕方ないですよね。自由度が奇数個の3軸で、2項演算を考えようとすると必ずハブられがちな軸は出てきますよ。やっぱりy軸なんですね z軸優遇なのは球座標系が前提にあるからなのかなと思います ＝＝＝＝＝＝＝＝ こういう議論をですね、たった2行の式を書くだけで、興味のある人には自力で全部理解させる まったくクドくないwikipediaのこのページの書き方は教育者として素晴らしいなあ！って思ったんです
「実は半整数も許されている」っつって、半整数の代数を始めると
いつの間にか軌道角運動量ではなくスピン角運動量の話にすり替えられていたんだな。 ところで、「軌道角運動量」の代数として、3次の行列はよく見かけるが、5次以上の奇数次の行列の例はほとんど見かけない。 というか分子の話になってしまうらしい。そういうのはノーサンキューだというのに。 で、軌道角運動量が整数しか許されないのなら、行列は奇数次に限られて偶数次の行列は考えるだけ無意味ということになるんだな。 ところで、スピンはよく±1/2の量子数を見かけるが ±3/2もどうやら許容されるらしい。 ということは、半整数の「スピン角運動量」の量子数としては、対応する偶数次の行列の演算を考えることに意味はあるということか。 ところで、sxとsyの積がszに比例するのは2行2列のパウリ行列のみらしく 一般の次数においては、交換関係[sx,sy]がszに比例(線形従属)するらしい。危うく間違えるところだった。 こういう「代数」は何という呼び名なのだろう？これが行列力学というものなのか？ 色々と定義が紛らわしい聞く。 その上、この「代数(微分演算子としてではなく)」は主量子数や調和振動子の量子数にも当てはめられるもんなんだろうか。なんかすごく今更感があるが、興味が出てきた。 磁気量子数に対する代数もあるのだろうか ＝＝＝＝＝＝＝ そんなことを考えながら、散歩中にARIAの主題歌を聞いていたら ふと、灯里の友達3種類の業種の人たちについて 「そういえばウンディーネを合わせると3種類ではなく4種類になるんだな」と気づかされた。 ウンディーネ、サラマンダー、ノーム、シルフ 4か なんだろう、男3人に女1人、まるで四次元時空の擬人化のようだ。 ノームが上下方向は確定として シルフとサラマンダーは直交座標系だろうか いや、アクア自体が(円柱座標系ではなく球座標の方の)極座標系としてあらわされたとしても、アクアの住民にとってはほぼほぼ直交座標系だろう。 実は惑星ではなく灯里たち人間を基準にした極座標なのではないか。 だとすると仰角(ピッチ)がサラマンダーで、何時の方向とかいうアレ(ヨー)がシルフのような気がする。 ウンディーネが時間軸担当なのは当然だ。いやあるいは世界線担当でもいいかもしれない。 ヒルベルト空間になってしまうだろうか？
20160905日記の続きです。
エネルギー固有値まで求められました。 これをもとに固有ベクトルを計算して、3軸それぞれの平面波の波動関数を求めたいと思います。 まずはx方向から順に行きましょう。 固有エネルギーを求めたのですから、2対の行列の式は、2対の永年方程式になっているはずで 1つの行列方程式の1行目と2行目の式は、固有エネルギーにおいて、同じ式を意味することになります。 A1とA4の比を求めることが固有ベクトルということになりますが とりあえず文字が多いので運動量については エネルギーについては としますと、 の下のほうの式に着目して A1とA4の比はA4/A1=tanθ、 同様に A3/A2=tanθになります。 このことから、AとA'を独立の規格化定数として、2つの波動関数が計算でき となります。 Eと-Eを置き換えて この式と この式になります。 今度は上の方に着目して解きますと A1/A4=-tanθ=A2/A3になるので、波動関数は wtの符号が先ほどと逆になっていることに注意してください。 同様にしてy,z方向の波動関数も求めてしまいますと y方向は z方向は となります。 ここで、ψの添え字であるx,y,zは方向、↑↓矢印はスピン、＋は粒子、－は反粒子を意味しています。 また、スピン演算子として、4次に拡張した以下のようなパウリ行列を用いると このように自分自身へのブーメランになるので、検算になると思います。 数検よりもこっちのが楽しすぎて困る  にほんブログ村
これディラック方程式。 行列をほどいて、連立方程式にするとこう。 これに、x,y,zそれぞれの方向に進行している平面波の波動関数 をぶち込みますと Ψxとψyは「Ψ1とΨ4」「Ψ2とΨ3」の2元連立方程式に分離でき Ψzは「Ψ1とΨ3」と「Ψ2とΨ4」の2元連立方程式のペアに分離できます。 x方向進行 y方向進行 z進行方向 これら3つは、どれも固有値問題を彷彿とさせる形になっていて、4つの波動関数の振幅Aが任意であるように永年方程式になる条件を探ると、案の定固有値問題に帰着します。 固有エネルギーħwを固有値として求めると、3つとも この2パターン(4パターン？)のいずれかになって このようなエネルギー固有値が算出されます。 ħwがエネルギー固有値E(アインシュタインの関係式)であるように、ħkは進行方向の運動量p(ド・ブロイの関係式)なので上の式は このように変形でき、クラインゴルドン方程式の着想そのものです。 固有ベクトルは工事中なのでもうちょっと待って～＞＜ まあとりあえずx,y,zの進行方向で、同じ固有値を算出する行列が、可能なすべてのパターンを示してくれて一安心です。 波動関数がΨ1・3対2・4だったりΨ1・4対2・3だったりするのは、まあ特に気にすることはないと思います。 もともと同種の素粒子には顔がないので、属性がすべて同じなら区別できませんしねせんしね ところで、お気づきかと思いますが、クラインゴルドン方程式をディラック方程式に改良してもなお、マイナスのエネルギーの存在という問題は解決していません。 これが未発見の粒子「反粒子」の予言につながるのです。 世は動物園じだ～い 未知なる素粒子を求めて探求するじだ～い  にほんブログ村
これが こうなって こうなる。 波動関数だってベクトルになりたいんだ！夢ならーたくさんぱぱっぱー しかしなんで運動量まで1階微分になったし。 そんなんで波動方程式が維持できるのかよぉーヨーﾖｰ
とりあえず計算した
実数行列だけでも数値的に計算可能 どうも、筋が通っていればいろんな定義方法があるらしく 僕が計算しているのは「ディラック表現 」というものらしい
まずはガンマ行列の定義から。 内積±1(反交換関係=2g)っぽいやつ 0番目だけはγ自身がエルミート、それ以外のγ自身は歪(反)エルミートで、 0番目以外は自身の2乗がマイナスの単位行列になることに注意 そこんとこ、前回のディラック行列aとは違うよね 内積0(反交換関係=0)っぽいやつ 今ならパウリ行列も詰め合わせてごらんのお値段！(書いてない)やっすいでしょー ちょっと怖い話 エルミート行列同士の積が歪エルミート行列になるのは、 4次行列に始まったことじゃなく、2次のパウリ行列ですでに起きていたんやでぇ・・・ﾄﾞﾛﾛﾛﾛﾛ・・・  にほんブログ村
まずはディラック行列の定義から。 内積1(反交換関係=2)っぽいやつ 内積0(反交換関係=0)っぽいやつ パウリ行列も詰め合わせておきました。 ちょっと怖い話 エルミート行列同士の積が歪エルミート行列になることがあるんですよ・・・ﾄﾞﾛﾛﾛﾛﾛ・・・ にほんブログ村
|
カレンダー
カテゴリー
ブログランキング
ブログランキング参戦中
よかったらポチッとお願いします＾＾
最新CM
[12/30 buy steroids credit card]
[09/26 Rositawok]
[03/24 hydraTep]
[03/18 Thomaniveigo]
[03/17 Robertaverm]
最新記事
(01/01)
(03/19)
(03/18)
(03/18)
(02/23)
(02/14)
(02/12)
(01/03)
(09/23)
(09/23)
(02/11)
(05/30)
(05/28)
(05/28)
(05/27)
最新TB
プロフィール
ＨＮ：
量子きのこ
年齢：
44
ＨＰ：
性別：
男性
誕生日：
1981/04/04
職業：
WinDOS.N臣T
趣味：
妄想・計算・測定・アニメ
自己紹介：
日記タイトルの頭についてるアルファベットは日記の番号です
26進数を右から読みます 例：H→7番目、XP→15(P)×26＋23(X)=413番目。 A=0とする仕様につき一番右の桁はAにできませんのでご了承くださいｽﾞｺｰ
バーコード
ブログ内検索
アーカイブ
最古記事
(05/11)
(05/11)
(05/13)
(05/13)
(05/13)
(05/13)
(05/13)
(05/13)
(05/14)
(05/14)
(05/14)
(05/14)
(05/16)
(05/16)
(05/16)
アクセス解析
|
