File failed to load: https://cdnjs.cloudflare.com/ajax/libs/mathjax/2.7.0/extensions/TeX/mathtools.js

2025-02-25

道志みちから大室山ピストン


野原吊橋から大室山まで。これまで3回トライして時間切れ、冬装備忘れなどで撤退していたので今回はチャンと準備してチャンと登頂しました! 行きは黙々と登り通しでしたが、帰りは山頂で出会った牛久の藤田さんと会話しながら下山できたので楽しかったです。懸案だった大室山に登頂できたのでこれで一安心。いろいろ勉強させてもらいました。ルートなど詳しくはこちらから。

2025-02-23

テニス:2024年度の活動報告

地元の市民大会の話ですが、団体戦は1部残留、個人戦はシニアで単複準優勝。シニアのカテゴリだけど漸くシングルスで結果が出ました。来年度は団体戦で勝てるよう頑張ります!

2025-02-22

高校の野球部OBOG会に参加

昨年末に2年上の先輩から案内があり、先週、本当に久しぶりに同期の仲間と再会、懐かしい先輩方とも歓談できました!


これを機会にンッ十年振りに連絡が取れた同期のメンバーもいるので久しぶりに同期会しようということになりました。もうみんないい年なので少しは落ち着いてきたかな。

2025-02-19

13. ウィグナーの D 関数とその応用 vol.5

前回に引き続き原子内電子の光吸収・発光過程の選択則について考える。選択則は放射過程の行列要素
α|Hint|β=emα|Ap|β=ieα|A(H0xxH0)|β
にウィグナー-エッカルトの定理を適用して導出できる。前回はゼロ次近似から電気双極子遷移の選択則を求めたが、今回は1次近似から磁気双極子遷移と電気四重極子遷移の選択則を求める。

磁気双極子遷移と電気四重極子遷移の選択則

  微小因子 kx について1次までのオーダーで行列要素(13.85)を表すと
α|Hint|β=emα|Ap|βeα|(1+ikx)Aωˆe˙x|β=α|H(0)int|β+α|H(1)int|β
と書ける。ただし、α|H(0)int|β前回求めたゼロ次のオーダーの項
α|H(0)int|β=ieˆeaAω(EαEβ)α|xa|β=e˙Aα|x|β=eEα|x|β
である。ただし、E=tA=˙A は外部電場を表す。ベクトル・ポテンシャル A
A=ˆeAei(ωtkx)=ˆeAωeikx
とパラメータ表示される。ただし、Aω=Aeiωt であり、角運動量 ω はエネルギー保存則から
ω=EαEβ=ωαωβ
と決まる。ここでは、放射ゲージを
ϕ=0,    A=0
採用していることに注意しよう。1次のオーダーの行列要素は
α|H(1)int|β=iekaˆebAωα|xa˙xb|β=iemkaˆebAωα|xapb|β=ie2mkaˆebAωα|(xapb+pbxa)|β
と表せる。ただし、xapb を反対称成分と対称成分に
xapb=12(xapbpbxa)+12(xapb+pbxa)=i2δab+12(xapb+pbxa)
と分離して、関係式 kaˆebδab=0 を用いた。対称成分はハイゼンベルク方程式
1mp=˙x=i[H0,x]
を用いて
xapb+pbxa=m2(xa˙xb+˙xbxa)=im2xa(H0xbxbH0)+m2˙xbxa=im2(H0xaxbxaxbH0)m2(˙xaxb˙xbxa)
と計算できる。よって、1次のオーダーの行列要素(13.97)は
α|H(1)int|β=e2kaˆebAωα|[H0,xaxb]|βie2mkaˆebAωα|(paxbpbxa)|β
と書ける。電気双極子近似(13.90)の場合と同様に、右辺の第1項は
e2kaˆebAωα|[H0,xaxb]|β=e2kaˆebAω(EαEβ)α|xaxb|β=e2ka˙Abα|xaxb|β=ie2aEbα|xaxb|β=ie2aEbα|Tab|β
と変形できる。ただし、Tab13.2節で定義した階数2の対称テンソル
Tab=xaxb13δabx2
である。外部電場に対して aEbδab=E=0 が成り立つことに注意しよう。式(13.100)右辺の第2項は
ie2mkaˆebAωα|(paxbpbxa)|β=ie2mkaˆebAωα|ϵabcϵklcpkxl|β=ie2mϵabckaˆebAωα|Lc|β=e2mϵabcaAbα|Lc|β=e2mBα|L|β
と変形できる。ただし、関係式 B=×AL=x×p を用いた。

 以上、まとめると
α|H(1)int|β=ie2aEbα|Tab|β+e2mBα|L|β
と求まる。右辺の第1項、第2項はそれぞれ電気四重極子遷移、磁気双極子遷移を記述する。電気双極子遷移の場合と同様に、第2項の選択則はクレブシュ-ゴルダン係数 C1llamm で決定される。ただし、磁気双極子遷移の場合は状態のパリティが保存される。つまり、Δl=ll に対して、(1)Δl=1 が課される。言い換えると、磁気双極子遷移の選択則は
Δl=0Δj=0,±1    ただし (j,j)(0,0)Δm=0,±1    ただし Δj=0 の場合は (m,m)(0,0)
で与えられる。

 同様に、ウィグナー-エッカルトの定理を適用すると電気四重極子遷移の選択則はクレブシュ-ゴルダン係数
C2llAmm=2Alm|lm=(1)ll2lm2A|lm=δm,m+A(1)ll2lm2A|lm
から導ける。ただし、A=0,±1,±2 である。上式より磁気量子数 m に関する選択則は Δm=0,±1,±2 であることが簡単に分かる。磁気双極子遷移の場合と同じく、始状態と終状態のパリティは同じである。よって、軌道角運動量量子数 l に関する選択則は例外的な場合を除いて Δl=0,±2 と書ける。例外となる場合は下表で示すクレブシュ-ゴルダン係数 lm2A|lm+A の具体的な形から判別できる。
lA=2A=1l+2(l+m+1)(l+m+2)(l+m+3)(l+m+4)(2l+1)(2l+2)(2l+3)(2l+4)(l+m+1)(l+m+2)(l+m+3)(lm+1)(l+1)(l+2)(2l+1)(2l+3)l+1(l+m+1)(l+m+2)(l+m+3)(lm)l(l+1)(2l+1)(2l+4)(l2m)(l+m+1)(l+m+2)l(l+1)(2l+1)(2l+4)l3(l+m+1)(l+m+2)(lm1)(lm)2l(l+1)(2l1)(2l+3)(2m+1)3(l+m+1)(lm)2l(l+1)(2l1)(2l+3)l1(l+m+1)(lm2)(lm1)(lm)(l1)l(2l+1)(2l+2)(l+2m+1)(lm1)(lm)(l1)l(2l+1)(2l+2)l2(lm3)(lm2)(lm1)(lm)(2l2)(2l1)2l(2l+1)(l+m)(lm2)(lm1)(lm)(l1)l(2l1)(2l+1)
lA=0l+23(l+m+1)(l+m+2)(lm+1)(lm+2)(l+1)(2l+1)(2l+3)(2l+4)l+1m3(l+m+1)(lm+1)l(l+1)(l+2)(2l+1)l3m2l(l+1)l(l+1)(2l1)(2l+3)l1m3(l+m)(lm)(l1)l(l+1)(2l+1)l23(l+m1)(l+m)(lm1)(lm)l(2l2)(2l1)(2l+1)
lA=1A=2l+2(lm+1)(lm+2)(lm+3)(l+m+1)(l+1)(l+2)(2l+1)(2l+3)(lm+1)(lm+2)(lm+3)(lm+4)(2l+1)(2l+2)(2l+3)(2l+4)l+1(l+2m)(lm+1)(lm+2)l(l+1)(2l+1)(2l+4)(lm+1)(lm+2)(lm+3)(l+m)l(l+1)(2l+1)(2l+4)l(2m1)3(lm+1)(l+m)2l(l+1)(2l1)(2l+3)3(lm+1)(lm+2)(l+m1)(l+m)2l(l+1)(2l1)(2l+3)l1(l2m+1)(l+m1)(l+m)(l1)l(2l+1)(2l+2)(lm+1)(l+m2)(l+m1)(l+m)(l1)l(2l+1)(2l+2)l2(lm)(l+m2)(l+m1)(l+m)(l1)l(2l1)(2l+1)(l+m3)(l+m2)(l+m1)(l+m)(2l2)(2l1)2l(2l+1)
これらの値はクレブシュ-ゴルダン係数の一般形 (Racah 公式) 
j1m1j2m2|JM=δM,m1+m2(2J+1)(J+j1j2)!(Jj1+j2)!(j1+j2J)!(j1+j2+J+1)!×(J+M)!(JM)!(j1+m1)!(j1m1)!(j2+m2)!(j2m2)!×k((1)kk!(j1+j2Jk)!(j1m1k)!(j2+m2k)!×1(Jj2+m1+k)!(Jj1m2+k)!)
から求めた。ただし、整数 k の和は階乗をとる数がすべて非負である k だけに限られる。この公式は M0 かつ j1j2 の場合に適用されるが、それ以外の場合は関係式
j1m1j2m2|JM=(1)Jj1j2j1m1j2m2|JMj2m2j1m1|JM=(1)Jj1j2j1m1j2m2|JM
から求まる。

2025-02-18

ウッドデッキ塗装5

去年はサボってしまったウッドデッキ塗装を久しぶりに行いました。寒いけど雨の少ない冬にやるのもアリですね。

塗装前

ほうきで掃いて水拭き

塗装後

相変わらず水性ペンキで塗っています。刷毛の再利用もできるし当分これでいいかな。

2025-02-17

13. ウィグナーの D 関数とその応用 vol.4

13.3 ウィグナー-エッカルト型の応用例


前節ではウィグナー D 関数を用いてコンパクトなリー群に関するウィグナー-エッカルトの定理
R,α|QA|R,m=CRRRAmαR||QR||R
を示した。ただし、還元行列要素 R||QR||R
R||QR||R=β,nCRRRBnβR,β|QB|R,n(dimR)(G の体積)
で定義される。今節ではウィグナー-エッカルトの定理の応用例として次の3つを取り上げる。
  1.  光子の吸収・放出についての選択則
  2.  ハドロン・スペクトル
  3.  ケース-ガシオロウィッツ-ワインバーグ-ウィッテン (Case-Gasiorowicz-Weinberg-Witten) の定理

13.3.1 光子の吸収・放出の選択則


2.1節で解説したように、水素原子のハミルトニアンは
H0=p22me2r
で近似できる。ここで、pm は電子の運動量と質量、r は原子核からの距離を表す。H0 は球対称なので
[La,H0]=0
を満たす。ただし、La (a=1,2,3) は角運動量演算子である。固有状態 |α のエネルギー固有値を Eα とする。つまり、H0|α=Eα|α とおく。(13.80)から LaH0|α=H0La|α=EαLa|α が分かる。これは、La|α|α と同じエネルギーを持つことを意味する。同様に、任意の状態の回転 eiLaθa|α も同じエネルギーを持つことが分かる。すなわち、エネルギー準位は角運動量の表現を成す多重項に分類できる。

 ある規約表現 R に属す状態 |α を考える。あるパラメータ θ を用いて別の状態 |β=eiLθ|α を定義する。状態 |β は表現 R の外に出ることはないので、特定の θ に対してβ|eiLθ|α0 となる。例えば、状態を適当に規格化すると β|eiLθ|α=α|eiLθeiLθ|α=1 が満たされる。もし |βR であれば、群の要素の合成が成り立たず、β|eiLθ|α=0 となる。よって、縮退状態は回転群の既約表現で与えられることが分かる。SU(2) 群の既約表現と縮退状態について詳しくは1.3節も参照されたい。

 ここで、ハミルトニアン H0 に相互作用項を導入する。放射ゲージを導入すると外部電磁場のスカラー・ポテンシャル ϕ とベクトル・ポテンシャル A はゲージ条件
ϕ=0,    A=0
に従う。このとき、相互作用項を含むハミルトニアンは
H=(p+eA)22me2r2=H0+e2m(pA+Ap)+e22mA2
と書ける。原子核内の電子に働くクーロン力に比べて外部電磁場の影響は小さいので A2 の項を無視すると、相互作用項は
Hint|Ψ=e2m(pA+Ap)|Ψ=emAp|Ψ
と表せる。ただし、ゲージ条件 A=0 を用いた。電子の運動量のみ考慮すると電流密度は J=emp とおけるので、上式から相互作用項 Hint は基本的に Hint=emAp=AJ で与えられることが分かる。ハイゼンベルク方程式を用いると、演算子 1mp
1mp=˙x=i[H0,x]
と書ける。よって、Hint の行列要素は
α|Hint|β=ieα|A(H0xxH0)|β
で与えられる。ただし、|α, |β は放射現象の終状態と始状態を表す。上式でベクトル・ポテンシャル A
A=ˆeAei(ωtkx)=ˆeAωeikx
とパラメータ表示される。ここで、Aω=Aeiωt であり、角運動量 ω はエネルギー保存則から
ω=EαEβ=ωαωβ
と決まる。放射ゲージ条件(13.81)より、偏光ベクトル ˆe は関係式
ˆek=0
を満たす。放射現象において光の波長は λ105 cm のオーダーであり、これは原子のサイズ r108 cm に比べて十分大きい。つまり、kx2πλr1031 と見積もることができる。微小因子 kx の1次近似でベクトル・ポテンシャルは
A=ˆeAω(1+ikx)+O((kx)2)
と表せる。ゼロ次近似 eikx1 で行列要素(13.85)は
α|Hint|βα|H(0)int|β=ieˆeaAω(EαEβ)α|xa|β=e˙Aα|x|β = eEα|x|β
と表せる。ただし、E=˙A は電場を表す。このゼロ次近似は古典電気力学での電気双極子近似に対応する。

 水素原子の物理状態は SO(3) 群のユニタリー既約表現で記述される。前節で見たように、これらは |l,m でラベルされる。ただし、l は軌道角運動量量子数、m は磁気量子数を表す。ウィグナー-エッカルトの定理(13.77)を適用すると、行列要素(13.90)は
eEal,m|xa|l,m=eEaC1llamml||x(1)||l
と表せる。ただし、x^{(1)}x^a (a = 0, \pm 1) が l=1 テンソル(13.65)として変換することを表す。これは放射現象の選択則がクレブシュ-ゴルダン係数 C^{1 l l^{\prime *}}_{a m m^\prime} で評価されることを意味する。クレブシュ-ゴルダン係数の慣習的な表示法 \bra j_1 m_1 \, j_2 m_2 | J M \ket を用いると C^{1 l l^{\prime *}}_{a m m^\prime}
\begin{eqnarray}    C^{1 l l^{\prime *}}_{a m m^\prime}    & = & \bra 1 a \, l m | l^\prime m^\prime \ket    \nonumber \\    &=& (-1)^{l^\prime - l -1 } \bra l m \, 1 a | l^\prime m^\prime \ket    \nonumber \\    &=& \del_{m^\prime , m + a} \,(-1)^{l^\prime - l -1 } \bra l m \, 1 a | l^\prime m^\prime \ket     \tag{13.92} \end{eqnarray}
と関係付けられる。これより、明らかに磁気量子数 m の選択則は
    \Delta m \, =  \, 0, \pm 1    \tag{13.93}
となることが分かる。ただし、\Delta m = m - m^\prime である。より厳密には、ゼロとならないクレブシュ-ゴルダン係数 \bra l m \, 1 a | l^\prime m + a \ket は以下の表で与えられる。
\begin{array}{|c|ccc|} \hline   l^\prime & a=1 & a=0 & a=-1 \\ \hline   l+1 & \sqrt{\frac{(l+m+1)(l+m+2)}{(2l+1)(2l+2)}} & \sqrt{\frac{(l+m+1)(l-m+1)}{(2l+1)(l+1)}} & \sqrt{\frac{(l-m+1)(l-m+2)}{(2l+1)(2l+2)}} \\   l & - \sqrt{\frac{(l+m+1)(l-m)}{2l(l+1)}} &  m \sqrt{\frac{1}{l(l+1)}} & \sqrt{\frac{(l+m)(l-m+1)}{2l(l+1)}} \\   l-1 & \sqrt{\frac{(l-m-1)(l-m)}{2l(2l+1)}} & - \sqrt{\frac{(l+m)(l-m)}{l(2l+1)}} & \sqrt{\frac{(l+m-1)(l+m)}{2l(2l+1)}} \\ \hline \end{array}

2025-02-16

夏と冬のパノラマ台

去年の夏と今年の1月に精進湖畔からパノラマ台に行ったので記録しておきます。夏はアメリカから一時帰国中の妹家族と1泊2日の旅行の際に寄りました。





夏は陽が昇ると雲が湧いてくるので早朝登山がおススメです。

 先月にテニス仲間と登ったときは冬らしく富士山くっきり見えました。






4年前に黒岳、御坂山を登ったときに出会った方に教えてもらったパノラマ台。湖畔の駐車場からアクセスもいいのでビギナーにもおすすめです!

2025-02-15

22年振り?! ハーレムのアパートでお世話になったTさんと再会

メールボックスの迷惑フォルダに紛れていたメッセージから久しぶりに懐かしい知人と先日再会しました。1999年に私が単身NYに乗り込んだ時にお世話になったアパートの大家さんです。京都出身のエネルギッシュな方で、すぐに意気投合して、2人で近くの公園でよくサッカーやりました。公園にいた多国籍メンバーと草サッカーやって削りあったりしましたね~。その他、イチローの試合見に行ったり、郊外のバザーに連れて行ってもらったり色々お世話になりました。マンハッタン・ハーレムのど真ん中のアパートで、夜に発砲音があったり、うるさすぎて勉強に集中できなかったので少しして私はクィーンズに引っ越しましたが、その後もたまにチャイナタウンでランチしたりしました。9.11のあと少ししてTさんがタイに行くと言って突然居なくなったのには驚きました。その際、車譲ろうかとおっしゃってもらったのですが、残念ながらあの時私には車を持つ余裕がありませんでした。
 
 その後、ぱったり音信不通になりましたが今年1月末にメッセージを受け取りました。タイ語をマスターして現地で運転手付きの生活をされていたそうですが、クーデターの内乱でどうにもならず、その後、中国の雲南大学で勉強されたそうです。で、今はロスで奥様と子供3人とのセミリタイア生活。今年還暦になるのに6,9,12歳の子供がいらっしゃるとのこと、相変わらずエネルギッシュですね!もちろん大谷のワールドシリーズも現地観戦したそうです。アメリカ、タイ、中国の話以外にもサッカーの話などで盛り上がって楽しい時間を過ごせました。年に一度チケットの安い時期に帰国するそうなのでまた来年会いましょう!

2025-02-14

電気四重極子遷移の選択則に関わるクレブシュ-ゴルダン係数

原子内電子の放射現象に関する選択則は例えばこちらにまとめられているが、その導出にはよく知られているようにウィグナー-エッカルトの定理が適用される。この定理の詳細はここでは割愛するが、結論としては電気双極子遷移あるいは磁気双極子遷移の選択則は、任意の角運動量固有状態 | l, m \ketl=1 の角運動量固有状態 |1, a \ket (a = 0, \pm 1) との合成で定義されるクレブシュ-ゴルダン係数 \bra lm \, 1 a | l^\prime m \!+\! a \ket の有無によって決定できる。ただし、軌道角運動量量子数l はゼロ以上の整数であり、磁気量子数 mm=-l, -  l \! +\! 1, \cdots , l\! -\! 1, l の値をとる。ただし、電子のスピンを含めると l は 半奇数を含む全角運動量量子数 j に置き換えれられる。ゼロとならない \bra lm \, 1 a | l^\prime m \!+\! a \ket は下表で与えられる。
\begin{array}{|c|ccc|} \hline   l^\prime & a=1 & a=0 & a=-1 \\ \hline   l+1 & \sqrt{\frac{(l+m+1)(l+m+2)}{(2l+1)(2l+2)}} & \sqrt{\frac{(l+m+1)(l-m+1)}{(2l+1)(l+1)}} & \sqrt{\frac{(l-m+1)(l-m+2)}{(2l+1)(2l+2)}} \\   l & - \sqrt{\frac{(l+m+1)(l-m)}{2l(l+1)}} &  m \sqrt{\frac{1}{l(l+1)}} & \sqrt{\frac{(l+m)(l-m+1)}{2l(l+1)}} \\   l-1 & \sqrt{\frac{(l-m-1)(l-m)}{2l(2l+1)}} & - \sqrt{\frac{(l+m)(l-m)}{l(2l+1)}} & \sqrt{\frac{(l+m-1)(l+m)}{2l(2l+1)}} \\ \hline \end{array}
この表は例えば


に掲載されている。上の表で l=m=0 (あるいは |a| < 1 かつ |l^\prime - l | < 1) の場合に係数が例外的にゼロとなることに注意しよう。

 一方、電気四重極子遷移の選択則は任意の角運動量固有状態 | l, m \ketl =2 の角運動量固有状態 | 2, a \ket (a = 0, \pm 1 , \pm 2) との合成で定義されるクレブシュ-ゴルダン係数 \bra lm \, 2 a | l^\prime m \!+\! a \ket の有無によって決定できる。よって、上と同様の表があれば簡単に電気四重極子遷移の選択則を確認できるのだが、分かる範囲では見つからなかった。そこで、以下のクレブシュ-ゴルダン係数の一般形 (Racah 公式) からこれらの係数を求めることにした。
\begin{eqnarray}    &&    \bra j_1 m_1 \, j_2 m_2 | J M \ket    \nonumber \\    &=& \del_{M, m_1 + m_2}  \sqrt{\frac{(2J+1)\, (J+j_1 - j_2 )! \, (J - j_1 + j_2 )! \,( j_1 + j_2 -J )! }{(j_1 + j_2 + J + 1 )!}}    \nonumber \\    && \times \sqrt{(J+M)! \, (J-M)! \,(j_1 + m_1 )! \, ( j_1 - m_1 )! \, (j_2 + m_2 )! \, ( j_2 - m_2 )!}    \nonumber \\    && \times \sum_k \left( \frac{(-1)^{k} }{k! \,(j_1 + j_2 -J - k )! \,(j_1 - m_1 - k )! \,( j_2 + m_2 - k )! } \right.    \nonumber \\    && \hspace{4cm} \times \left. \frac{1}{ ( J- j_2 + m_1 + k )! \,( J- j_1 - m_2 + k )! } \right)    \tag{1}\end{eqnarray}
ただし、整数 k の和は階乗をとる数がすべて非負である k だけに限られる。この公式は M \ge 0 かつ j_1 \ge j_2 の場合に適用されるが、それ以外の場合は関係式
\begin{eqnarray}    \bra j_1  - \! m_1 \, j_2  - \! m_2 | J \, - \!M \ket &=& (-1)^{J-j_1 -j_2 } \bra j_1 m_1 \, j_2 m_2 | J M \ket      \tag{2} \\    \bra j_2  m_2  \, j_1  m_1 | J \, M \ket &=& (-1)^{J-j_1 -j_2 } \bra j_1 m_1 \, j_2 m_2 | J M \ket     \tag{3} \end{eqnarray}

2025-02-13

放射現象における磁気双極子遷移と電気四重極子遷移の導出

古典電磁気学の1つの重要なトピックとして電気四重極子あるいは多重極子の導出がある。
静電場の場合はスカラー・ポテンシャルの球面調和関数展開あるいはテイラー展開に相当するので自然に理解できるが、磁場を含める場合はベクトル・ポテンシャルを導入して電気双極子近似(ゼロ次近似)より高次の近似を考える必要がある。このとき1次のオーダーの項から磁気双極子と電気四重極子が出てくるが、これらはぞれぞれ電流ベクトルと位置ベクトルの積の非対称成分と対称成分に相当する。この話は例えばジャクソン電磁気学の第9章




に出ていて一応、大学の古典電磁気学で学習する内容となっている。ただし、このような近似が成立するのは基本的に原子・分子・結晶内などミクロな世界の話なので当然ながら量子電磁気学の領域であり、量子力学的な取り扱いが必要となる。

 一方、量子力学の教科書では大抵、ゼロ次近似での電気双極子相互作用の導出までが紹介され、磁気双極子・電気四重極子は古典電磁気学でやったから大丈夫だよね、という感じで学生もあぁあれね、という感じでそのままになってそれぞれ専門分野に進んで行くことが多い(個人的印象)。最近までそんな記憶もすっかり忘れていたが、久しぶりに原子内電子と外部電磁場の相互作用を考える機会があったので磁気双極子相互作用と電気四重極子相互作用を量子力学的に導出できるか手元の参考書・教科書で調べてみた。




しかしながら、分かる範囲では何処にも載っていなかったので自力で導出してみた。詳細は以下の通り。