侮辱罪

5月18日に，「侮辱罪」の厳罰化法案が衆議院の法務委員会を通過した。

これは「刑法等の一部を改正する法律案」のことであり，その改正法案のほとんどは懲役，禁錮を拘禁刑に統合整理することに関わるものだ。そこに，ネットでの誹謗中傷に起因した自殺事件を口実とした侮辱罪の厳罰化が混入されたものになっている。

改正法案の第一条は，「刑法（明治四十年法律第四十五号）の一部を次のように改正する。第二百三十一条中「拘留又は科料」を「一年以下の懲役若しくは禁錮若しくは三十万円以下の罰金又は拘留若しくは科料」に改める。」である。

元の刑法第二百十一条は，第三十四章「名誉に対する罪」の一部である。改正法律案の第二条の後段の方で，上記下線部の「懲役若しくは禁錮」が「拘禁刑」に置き換えられる。それにしても，若しくはと又はの使い方がややこしい。

第三十四章　名誉に対する罪
（名誉毀き損）
第二百三十条　公然と事実を摘示し、人の名誉を毀き損した者は、その事実の有無にかかわらず、三年以下の懲役若しくは禁錮又は五十万円以下の罰金に処する。
２　死者の名誉を毀損した者は、虚偽の事実を摘示することによってした場合でなければ、罰しない。
（公共の利害に関する場合の特例）
第二百三十条の二　前条第一項の行為が公共の利害に関する事実に係り、かつ、その目的が専ら公益を図ることにあったと認める場合には、事実の真否を判断し、真実であることの証明があったときは、これを罰しない。
２　前項の規定の適用については、公訴が提起されるに至っていない人の犯罪行為に関する事実は、公共の利害に関する事実とみなす。
３　前条第一項の行為が公務員又は公選による公務員の候補者に関する事実に係る場合には、事実の真否を判断し、真実であることの証明があったときは、これを罰しない。
（侮辱）
第二百三十一条　事実を摘示しなくても、公然と人を侮辱した者は、拘留又は科料に処する。
（親告罪）
第二百三十二条　この章の罪は、告訴がなければ公訴を提起することができない。
２　告訴をすることができる者が天皇、皇后、太皇太后、皇太后又は皇嗣であるときは内閣総理大臣が、外国の君主又は大統領であるときはその国の代表者がそれぞれ代わって告訴を行う。

どう考えても，侮辱という定義の曖昧な概念のせいで，政治家やデマゴーグ達への侮辱罪だという喚き声が言論の自由を萎縮させる構図が露骨に透けて見える。維新の親玉がスラップ訴訟で批判を抑え込もうとしている現実を見れば，全く安心できないのだが，どのマスコミからも大きな批判がなくすんなりと通ってしまいそうだ。

なお，この改正案の第二条は刑法，第三条は刑事訴訟法，第四条は刑事収容施設及び被収容者等の処遇に関する法律，第五条は刑事収容施設及び被収容者等の処遇に関する法律（続き），第六条は更生保護法，第七条は更生保護法（続き），第八条は更生保護事業法，第九条は更生保護事業法（続き），第十条は少年院法，第十一条は少年鑑別院法，第十二条は少年鑑別院法（続き）となっている。

P. S. 親告罪の第二百三十二条で，上皇の取り扱いはどうなるのかと心配したところ，天皇の退位等に関する皇室典範特例法の附則で，上皇は天皇の上皇后は皇太后の例によるとされていた。

個人送信（１）

国立国会図書館の「個人向けデジタル資料送信サービス（略称：個人送信）」が5月19日から開始された。

国立国会図書館のデジタルコレクションは，約10年前に，「国立国会図書館デジタル化資料」として古典籍資料（貴重書等），図書，インターネット資料を提供開始したのが源流となっている。その後，博士論文，録音映像資料，や様々なタイプの資料が加わってきた。

著作権など権利関係に問題がないものは，すでにインターネット公開されている。著作権保護期間が満了していない資料や著作権の確認が完了していない資料のうち，絶版などで入手困難なものは，これまで図書館送信資料として，準備の調った最寄りの図書館で閲覧することができた。

今回は，このカテゴリーに個人向けデジタル化資料送信サービスとして加わる。国立国会図書館に登録した個人が，図書館を経由することなく，入手困難資料約153万点を直接見ることができるようになる。ただし，当面，PCやタブレットでの閲覧に限られており印刷することはできない。

条件として，日本国内に居住していることと面倒なハードルをつけていた。さらに，国立国会図書館の利用者登録が必要ときた。早速申し込んでみると，本人確認書類として免許証，保険証，マイナンバーカードなどから2点が要求された。なんとかクリアすると，分かりにくい登録番号が振られ，確認終了まで5営業日待てと。これだから日本の公共セクターのデジタル化は終わっている。

とりあえず仮登録したところで，アクセスできるよとの指示があるが，パスワードではねられてしまう残念な仕様になっているようだ。パスワードを変更してみたが（できた），やはりアクセスできない。「ダメなものはダメ（オーケー，グァン姉妹）」ならば初めからそう言って下さい。

［１］個人向けデジタル化資料送信サービス　利用規約

堀河夜討

NHKの大河ドラマ「 鎌倉殿の13人」は，平家が滅亡して義経が追われる段階に入ってきた。

人形浄瑠璃や歌舞伎の演目と重なるところが多いのでたいへん勉強になる。「義経千本桜」をはじめとして，「ひらかな盛衰記」「一谷嫩軍記」「源平布引滝」「近江源氏先陣館」「鬼一法眼三略巻」「鎌倉三代記」「梶原平三誉石切」「御所桜堀河夜討」などがこのテーマに関連した主な作品であり，文楽劇場でも何度か観ている。

今回は，1185年に頼朝の命を受けた土佐坊昌俊が，京都の六条堀河邸または六条室町亭にいる源義経を夜討するという史実に基づいた話だったが，ドラマでは義経の正妻の郷御前（三浦透子）がからむということになっていた。これが，浄瑠璃作品になると，平時忠の娘が正妻の郷の君になって，弁慶の隠し子が出てくるとか，もっと話がややこしくなったりするわけだ。

写真：堀河夜討の図（勝川春章，神奈川県立歴史博物館より引用）

越中大門

 越中大門は富山県射水市にある北陸本線の駅だ。北陸新幹線の開通後は，並行在来線の廃止にともなって第3セクターのあいのかぜ富山鉄道が運営している。

ところで，Unreal Engine 5 は，エピックゲームズというアメリカのコンピュータゲーム・ソフトウェアの開発会社がつくった，ゲームエンジンである。3次元コンピュータグラフィックス生成用のプログラムを3Dエンジンとよんでいたが，この延長上の概念に相当する。

Unreal Engine 5 は2020年に発表された。この機能を利用して，Lorentz Dragoというイタリアの3D環境アーティストが個人でモデリング，テクスチャリング，ライティング，アニメーションのすべてを手がけたものが，越中大門の駅の3DCGアニメーションだ。YouTube で2分49秒のコンテンツだけれど，見た目にはまったく現実の駅を撮影したものといわれてもわからない。途中で昼のシーンが急に夜に変わったところではじめてコンピュータによる合成画像だと気付く。

この技術がメタバースに組み込まれると，視覚上はほとんど現実と区別のつかない仮想世界体験が得られるかもしれない。

越中大門は，子どものころ母の実家の滑川まで里帰りするとき，蒸気機関車が牽引する鈍行列車で通過したなじみ深い駅だ。当時の北陸本線の駅を確認すると次のようになっていて，越中大門はちょうど金沢と滑川の中間あたりにある。

（現ＩＲいしかわ鉄道）金沢−東金沢−森本−津幡−倶利迦羅−（現あいの風とやま鉄道）石動−福岡−西高岡−高岡−越中大門−小杉−呉羽−富山−東富山−水橋−東滑川−滑川

ものごころついてから小学校5-6年ごろまで通ったのだが，西高岡，東富山，東滑川はあまり記憶に残っていない。東滑川はまだ開業していなかったかもしれない。北陸線の電化が進んで蒸気機関車が廃止になるころには，里帰りについていくこともなくなってしまった。

何度かトンネルをくぐるのだが，そのたびに蒸気機関車の煙が入らないよう窓を閉めていた。高岡駅のホームなどで母がアイスクリームを買いに降りるのだが，時間内に帰ってこられるかどうかドキドキしながら待っている小学生であった。

写真：越中大門駅のホーム（あいの風とやま鉄道から引用）

AI自動音声合成

 夜のNHKニュースで，AI自動音声合成によるというキャプションが出るものが登場した。これは，VOICEPEAKを使っているのかと調べると，NHK放送技術研究所が開発した独自の日本語音声合成システムらしい。

今回の開発では、「漢字仮名交じり文」から「仮名文字と韻律記号」を自動的に生成し、それを「系列変換モデル」の入力データとすることで大量のデータを効率的に学習させ、日本語の合成音声の品質を向上させることに成功しました。

また、仮名文字と韻律記号を簡単に編集できるユーザーインターフェースや、口調をニュース調や会話調などに切り替えられる技術も開発し、さまざまな番組の演出要件への対応も可能にしました。

ということで違和感なく聞くことができるレベルになって，NHKの朝のニュースにも採用されはじめた。VRのアバターによるニュースが始まるまであと一歩だ。アナウンサーは放送局のかかえるタレント的存在になっているが，その傾向に拍車がかかる。

なお，VOICEPEAKは29,800円くらいで購入できて商用利用も可能だが，VOICEVOXという無料の中品質音声合成ソフトウェアも登場していた。

［１］音声合成ソフトの進化がすごい！（PC Watch 2022.4.16）

沖縄返還50周年

今日は，1972年5月15日の沖縄の施政権返還から50周年，沖縄旅行の話は以前書いている。

1972年は，3月の3日から5日まで大学入試があり，4月からは大阪での生活がはじまった年だ。阪大物理学科の同窓会でも入学50年を迎えてという話があったがどうなることやら。

というわけで，1972年にあった様々なイベントは全て50年目を迎えることになる。そこで，Wikipediaをざっと調べて，自分の印象に残っているものをリストアップしてみた。千日前デパートビル火災と日中共同声明の記事が日本語版Wikipediaの1972年にはなかったが，英語版の1972年にはあった。

入試前の2月にいろいろな事件が集中していてたいへんだった。大学入試で金沢と大阪を往復した雷鳥の中で，パイオニア10号の記事をみて感動したのをおぼえている。カール・セーガンの発案で，地球外知的生命に向けてのメッセージを記録した金属板が搭載されていたからだ。

2月02日　横井庄一帰国
2月03日　札幌オリンピック開幕

2月21日　ニクソン訪中
2月28日　あさま山荘事件終結
3月01日　ローマクラブ「成長の限界」
3月02日　パイオニア10号打ち上げ
3月07日　連合赤軍リンチ殺人事件発覚（山岳ベース事件）
4月16日　川端康成ガス自殺
5月13日　千日前デパートビル火災
5月15日　沖縄返還・沖縄県発足
5月30日　テルアビブ空港乱射事件
9月05日　ミュンヘンオリンピック事件
9月29日　日中共同声明　日中国交正常化
12月7日　アポロ17号打ち上げ（アポロ計画終了）
12月28日　金日成国家主席就任

写真：パイオニア10号に搭載された知的生命体へのメッセージ金属板

［１］日本復帰への道（沖縄公文書館）

仏滅の周期

 トイレのカレンダーを見ていると，周期性のある仏滅の位置が今月末あたりでズレている。

六曜は，暦に記載される日時・方位などの吉凶、その日の運勢などの事項である暦注のひとつである。先勝（せんしょう）・友引（ともびき）・先負（せんぶ）・仏滅（ぶつめつ）・大安（たいあん）・赤口（しゃっこう）の六種がこの順に毎日繰り返す。

調べてみると，六曜の周期性の乱れは次のルールによる。旧暦の月の朔日が一月と七月は先勝，二月と八月は友引，三月と九月が先負，四月と十月が仏滅，五月と十一月が大安，六月と十二月が赤口と，あらかじめ定まった六曜が朔日に割り当てられており，これから順に並べていくことになる。そこで旧暦の月Mと日Dがわかれば，(M+D-2) mod 6 を計算しよう。この剰余の集合 { 0, 1, 2, 3, 4, 5 } に対して，先ほどの六曜の集合 {先勝, 友引, 先負, 仏滅, 大安, 赤口} を順番に対応させれば，各月各日の六曜が決まる。

今年の5月30日は旧暦の皐月朔日（5月1日）になるので，大安にセットされる。このため仏滅が，新暦の5/7−5/13−5/19/−5/25−5/31 という系列からずれて，6/4になったというわけだ。

あまり役に立たない豆知識シリーズ終了。

LDLコレステロール（１）

人間ドック（２）からの続き

以前から，中性脂肪やLDLコレステロール値にしばしば警告が出ていたが，あまり気にすることはなかった。ところが，今回はいよいよ，要治療領域まであと一歩という説明をドクターから直接受けたので早速以前の結果を見直してみた。 

図：人間ドックにおける2001年以降の血清脂質値

この22年の人間ドックの検査値の平均は中性脂肪（TG=トリグリセリド）が TG=113±24 mg/dL（35-149），LDLコレステロール（悪玉）が LDL-C=139±11 mg/dL（70-139），HDLコレステロール（善玉）が HDL-C=61±4 mg/dL （40-）となった。総コレステロールは，TC=LDL-C+HDL-C+TG/5で求まり，TC=225±12 mg/dL（130-219）である。括弧内が正常値だ。なお 善玉を除いたnon-HDLコレステロールは，non-HDL-C = TC − HDL-C となっている。

中性脂肪とHDLコレステロールはほぼ問題ないのだが，これまでのLDLコレステロールの平均値が，許容値の70-139 mg/dL の上限に達し，上回ることの方が多かった。LDLコレステロールが140 mg/dL以上は高コレステロール血症に該当し，180 mg/dL以上が要治療対象だ。加齢や運動量の不足からか，定年後の3-4年でLDL-CやTCが上昇傾向にあるのがよろしくない。

高コレステロール血症の 2-3割は食物や生活習慣が原因で，7-8割がコレステロールをつくる肝臓の機能の問題らしい。肝臓の機能のうち遺伝的な原因が主ならばなかなか対応が難しいわけだ。ネットで検索するとちょっと怪しいサプリページがどんどんと引っかかってくる。

とりあえず，次の検査の3ヶ月後まで，動物性脂肪（飽和脂肪酸）を少しでも減らして，有酸素運動を始めろということか。

［１］脂質異常症（eヘルスネット）

［２］人間ドック学会判定区分（2021年度版）

人間ドック（２）

人間ドック（１）からの続き

昨日は，天理市立メディカルセンターの人間ドック健診日だった。

天理市の国民健康保険加入者に対する人間ドック補助制度がある，今年から条件に該当することになったので，早速申し込んでいた。メディカルセンターの受付時間は8:30から8:45までで，10分前につくと6番目だ。大手前病院では，家を6:00前に出て，7:15に病院に到着してやっと6番目前後だったので，えらい違いである。

天理市立メディカルセンターは，高井病院の福祉医療法人である高清会が運営していて，設備も比較的新しく普通の病院とあまりかわらない。健診を受けている人数がそもそも少ないので対応はよい。タニタの最新式体重計では，体脂肪率だけでなく身体各部の脂肪量や筋肉量まで算出してくれる。脂肪量は平均値なのだが，筋肉量が平均より2段階下である。どおりで懸垂ができなくなるわけだ。

気になった点は，(1) 眼底検査がなかなかうまくいかず，３回取り直した。これは私のせいなのか？(2) 聴力検査はかなりアバウトなのであった。(3) 胃Ｘ線（バリウム）は検査技師の方はとても丁寧で親切だったけれど，発泡剤を補助液なしにそのまま口に入れ，大量のバリウムで直接流し込む方式には閉口した。診断台で身体を支えるのもなかなか困難になってきたので，来年は胃カメラかな。

問診は，超音波検査を担当していたお医者さんに血液検査のデータを見ながら丁寧に説明してもらった。LDLコレステロールが170mg/dLと後一歩で要治療の180mg/dLに達するということで脅された。さっそく看護師の保健指導をうけつつ3ヶ月後の再検査を申し込むことに。

［１］バックス発泡顆粒（「透視開始に際して、造影剤投与開始直前あるいは投与開始後、年齢、胃内容積の個人差、造影の体位に応じて、約100~400mLの炭酸ガスの 発生量に相当する量を、少量の水または、造影剤と共に経口投与する」なので，補助液なしでもかまわないらしい。）

理想気体のエントロピー

 （１）熱力学的なエントロピー$S(U,V)$を，$dS=\frac{d'Q}{T}\ $で定義する。これを熱力学第一法則の中に埋め込むと$\ dU=TdS-pdV$，すなわち，$dS=\frac{dU}{T}+\frac{p dV}{T}$ となる。

ここに，理想気体の状態方程式$\ pV = nRT = N k_B T\ $と，単原子分子気体の内部エネルギーの表式，$U = n C_v T = n \frac{3}{2} R T = N \frac{3}{2} k_B T\ $ を代入する。

$dS= N k_B \Bigl( \frac{3}{2} \frac{dT}{T} + \frac{dV}{V} \Bigr )\quad \therefore \int dS = N k_B \Bigl( \frac{3}{2} \int \frac{dT}{T} + \int \frac{dV}{V} \Bigr )$ 

$S = S_0 + N k_B \Bigl( \frac{3}{2} \log \frac{T}{T_0} + \log \frac{V}{V_0} \Bigr) = S_0 +  N k_B \log \frac{T^{3/2}V}{T_0^{3/2}V_0}$

（２）一方，統計力学において，自由粒子の小正準集団のエントロピーは，ボルツマンの原理から，$S = k_B \log W(E) = k_B \log \frac{\partial}{\partial E} \Omega_0(E)\  \delta E$となる。

質量$m$の$N$粒子系のエネルギー$E$までの状態数$\ \Omega_0(E)\ $は，$N$粒子系の$6N$次元の位相空間の体積を$6N$次元細胞の体積$\ h^{3N}\ $と同一粒子が区別できないことによる因子$N!$で割ったものになり，運動量空間での半径$\ \sqrt{2mE}\ $の$3N$次元超球の体積の式を使うと，

$\displaystyle \Omega_0(E) = \frac{1}{h^{3N} N!} \int d{\bm q} \int _{\Sigma p_i^2/2m < E} d{\bm p} = \frac{V^N (2\pi m E)^{3/2}}{h^{3N}\Gamma(3N/2+1)}$

$W(E) = \frac{3N}{2} \frac{V^N}{N! (3N/2)!} \Bigl( \frac{2\pi m E}{h^2} \Bigr)^{3/2} \frac{\delta E}{E} =  \frac{3N}{2} \Bigl(\frac{V}{N}\Bigr) ^N\Bigl( \frac{2\pi m E}{3N/2} \Bigr)^{3N/2} e^{5N/2} \frac{\delta E}{E}$

$\therefore S(E) = k_B \log W(E) = N k_b \Bigl\{\frac{3}{2}\log \Bigl( \frac{4\pi m E}{3 h^2 N}\Bigr)  + \log \frac{V}{N} +\frac{5}{2} \Bigr\}$

$=N k_B  \log \frac{(2\pi m k_B T)^{3/2} V e^{5/2}}{N h^3}$

ここで，$\frac{1}{T} = \frac{\partial S}{\partial U}= \frac{\partial S}{\partial E} = \frac{3 N k_B }{2  E}$より，$\frac{E}{N}=\frac{3}{2}k_B T$を用いた。$S(E)$の式の$\log$の中身は無次元であり，$V/N$があるので，全体として示量変数ではないことが保証されている。

理想気体の熱力学的エントロピーと統計力学的エントロピーは，ともに $N k_B \log T^{3/2}V + const.$の形をしているが，定数部分まで含めて同じかどうかがよく理解できていない。

摩擦のあるカルノーサイクル（３）

 摩擦のあるカルノーサイクル（２）からの続き

摩擦のあるカルノーサイクルでクラウジウスの不等式を説明するためには，前回のように$W'$のなかの摩擦力による仕事$\delta w$として表現するかわりに摩擦力で生じた熱$\delta_{\rm q}=-\delta w < 0$として扱うこともできる。

仕事として表現すると：
等温過程 A$\rightarrow$B：（$Q_{\rm H}>0, \quad \delta w_{\rm AB}>0$）
$W'_{\rm AB}=\int_{V_{\rm A}}^{V_{\rm B}} pdV -\delta w_{\rm AB}=W_{\rm AB} -\delta w_{\rm AB}=Q'_{\rm H}$
等温過程 C$\rightarrow$D：（$Q_{\rm L}<0, \quad \delta w_{\rm DC}>0$）
$W'_{\rm CD}=\int_{V_{\rm C}}^{V_{\rm D}} pdV -\delta w_{\rm DC}=W_{\rm CD} -\delta w_{\rm DC}= Q'_{\rm L}$
熱として表現すると：
等温過程 A$\rightarrow$B：（$Q_{\rm H}>0, \quad \delta q_{\rm H}<0$）
$W'_{\rm AB}=\int_{V_{\rm A}}^{V_{\rm B}} pdV + \delta q_{\rm H} = Q_{\rm H} +\delta q_{\rm H} = Q'_{\rm H}$
等温過程 C$\rightarrow$D：（$Q_{\rm L}<0, \quad \delta q_{\rm L} < 0$）
$W'_{\rm CD}=\int_{V_{\rm C}}^{V_{\rm D}} pdV + \delta q_{\rm L} = Q_{\rm L} + \delta q_{\rm L}= Q'_{\rm L}$

カルノーサイクルにおいては，系に入る熱量を温度でわった，エントロピーに対応する状態量$\frac{Q}{T}$の和が保存していた。すなわち，$\dfrac{Q_{\rm H}}{T_{\rm H}}+\dfrac{Q_{\rm L}}{T_{\rm L}} = n R \log \dfrac{V_{\rm B}}{V_{\rm A}} + n R \log \dfrac{V_{\rm D}}{V_{\rm C}} = 0$

一方，摩擦のあるカルノーサイクルで，出入りする熱量に対して，温度で割ったものの和を考えると，$\dfrac{Q'_{\rm H}}{T_{\rm H}}+\dfrac{Q'_{\rm L}}{T_{\rm L}}=\dfrac{Q_{\rm H}+\delta q_{\rm H}}{T_{\rm H}}+\dfrac{Q_{\rm L}+\delta q_{\rm L}}{T_{\rm L}} = \dfrac{\delta q_{\rm H}}{T_{\rm H}}+\dfrac{\delta q_{\rm L}}{T_{\rm L}} \le 0$

これを一般化すると，可逆過程だけのサイクルについては $\displaystyle{ \oint \dfrac{d'Q}{T_{\rm ex}} = 0}$，不可逆過程を含むサイクルについては，$\displaystyle{\oint \dfrac{d'Q}{T_{\rm ex}} \le 0}$，ここで，$d'Q$は系が受け取る熱量で，$T_{\rm ex}$はその熱量を与えた熱源の温度である。これがクラウジウスの不等式。

図：不可逆過程におけるクラウジウスの不等式

準静的過程がわからない

エントロピーがわからないからの続き

熱力学の入門的教科書を手元に並べて呻吟している。

そういえば，教養課程で物理学科の専門科目として大学１年のときにクラス担任の国富信彦先生が担当したのが「物理学要論？」だった（科目名も忘れてしまった）。そこで，最初につまづいたのが応力テンソルと準静的過程だった。熱力学の初歩のところでは，覆水盆に返らずの話をしながら準静的過程の説明があったので，これは可逆過程なのか不可逆過程なのかどうなっているの？と混乱したのだった。

さて，並べているやさしい教科書は以下のとおり

１．フェルミ熱力学（エンリコ フェルミ，三省堂 1973）

２．熱・統計力学（戸田盛和，岩波書店 1983）

３．熱・統計力学の考え方（砂川重信，岩波書店 1993）

４．熱学入門（藤原邦男・兵藤俊夫，東京大学出版会 1995）

５．ゼロからの熱力学と統計力学（和逹三樹・十河清・出口哲生，岩波書店 2005）

フェルミには準静的過程というワードは出てこない。それに相当するものは熱平衡状態をつないでいく可逆過程で考えるという立場だ。戸田さんは，この本（熱力学）で扱う可逆過程は準静的過程に限定すると注意している。砂川さんは，力学的な過程も可逆過程に含め，可逆過程は必ずしも時間反転してたどる必要がないとしている。準静的過程は可逆過程の集合に含まれる。藤原さんや和逹さんは，可逆過程＝準静的過程としている。

広島大学の戸田昭彦さんは準静的過程と可逆過程に対してもっと細かな議論を展開していた。

摩擦のあるカルノーサイクル（２）

 摩擦のあるカルノーサイクル（１）からの続き

「エントロピーについての理解を図るため，不可逆過程の具体的な例を構成したい。そのためカルノーサイクルの等温過程においてのみピストンに散逸のある抵抗力=摩擦力が働くモデルを考える。この摩擦力は，ピス トンの運動方向と逆向きに作用し，その仕事はピストンに熱として放出され，カルノーサイクル の作業物質である理想気体には影響を及ぼさないものとする。この考察において作業物質の系=理想気体がする仕事は，ピストンを用いて測定されることに留意する。すなわち，ピストンに働く力の総和とピストンの変位の積によって作業物質系が「する」仕事や「される」仕事（＝負の「する」仕事）が定義される」

ということで，前書きをかいて計算をはじめてみたもののなかなか難渋するのであった。

図：摩擦のあるカルノーサイクルの散逸過程

等温過程 A $\rightarrow$ B：（$\delta_{\rm AB}>0$）

$W'_{\rm AB}=\int_{V_{\rm A}}^{V_{\rm B}}p dV - \int_{d_{\rm A}}^{d_{\rm B}} f dx = nRT_{\rm H} \log \frac{V_{\rm B}}{V_{\rm A}} - \delta w_{\rm AB} = Q_{\rm H} - \delta w_{\rm AB} = Q'_{\rm H}$

断熱過程 B $\rightarrow$ C：

$W'_{\rm BC}=\int_{V_{\rm B}}^{V_{\rm C}}p dV = \int_{\rm B}^{\rm C} -dU = U_{\rm B}-U_{\rm C} = n C_V (T_{\rm H}-T_{\rm L})$

等温過程 C $\rightarrow$ D：（$\delta_{\rm DC}>0$）

$W'_{\rm CD}=\int_{V_{\rm C}}^{V_{\rm D}}p dV - \int_{d_{\rm C}}^{d_{\rm D}} f dx = nRT_{\rm L} \log \frac{V_{\rm D}}{V_{\rm C}} - \delta_{\rm DC} = Q_{\rm L} - \delta w_{\rm DC}= Q'_{\rm L}$

断熱過程 D $\rightarrow$ A：

$W'_{\rm DA}=\int_{V_{\rm D}}^{V_{\rm A}}p dV = \int_{\rm D}^{\rm A} -dU = U_{\rm D}-U_{\rm A} = n C_V (T_{\rm L}-T_{\rm H})$

したがって，この摩擦のあるカルノーサイクルの効率は次の式で与えられる。

$W'_c = W'_{\rm AB}+W'_{\rm BC}+W'_{\rm CA}+W'_{\rm AD}$

$= nRT_{\rm H}\log \frac{V_{\rm B}}{V_{\rm A}} -\delta_{\rm AB}  + nRT_{\rm L} \log \frac{V_{\rm D}}{V_{\rm C}} - \delta_{\rm DC} = Q_{\rm H} -\delta_{\rm AB} + Q_{\rm L} - \delta_{\rm DC}$

$\eta' = W'_c/Q'_{\rm H}=1-\frac{-Q'_{\rm L}}{Q'_{\rm H}} \approx 1-\frac{-Q_{\rm L}}{Q_{\rm H}} \bigl( 1+\frac{\delta w_{\rm DC}}{-Q_{\rm L}} \bigr) \bigl(1+\frac{\delta w_{\rm AB}}{Q_{\rm H}} \bigr)$

$\le 1-\bigl( \frac{-Q_{\rm L}}{Q_{\rm H}}\bigr) =  1-\frac{T_{\rm L}}{T_{\rm H}} = \eta_c$

こうして，摩擦のあるカルノーサイクルの効率$\eta'$は，カルノーサイクルの効率$\eta_c$より小さくなる。しかし，このままでは，$\int \frac{d'Q}{T} \le 0$の説明にうまくつながらない。

（注）ここではサイクルに入る熱量をすべて正，サイクルが外部にする仕事を全て正にとる。

摩擦のあるカルノーサイクル（１）

エントロピーがわからないからの続き

熱力学の第二法則と仮に導入したエントロピーの違いを明らかにしたい。普通の教科書ではクラウジウスの原理やトムソンの原理によるわけだが，そのためには，不可逆過程の考察が必要となる。その一番簡単な例は，力学でもなじみのある摩擦現象だと思う。

摩擦現象は力学的に細かく詰めて考えると何だか複雑で面倒なことになるが，とりあえず，摩擦現象は運動や仕事が熱に変換されて，力学的エネルギーが熱エネルギーとして環境中に散逸する過程だと考えることにする。環境中の熱エネルギーが直接集まってきて力学的エネルギーになるような現象は，巨視的には観察されないので，摩擦現象は不可逆過程である。

そこで，これまで練習してきたカルノーサイクルに摩擦を導入すれば，理解がしやすいのではないかと考えた。どの教科書をみてもそんな具体的な議論はされていない。クラウジウスやトムソンにしたがって，より抽象的な熱機関（可逆機関，不可逆機関）の組み合わせでの議論が進んでいくわけだ。

というわけで，より具体的な摩擦のあるカルノーサイクルを構成してみることに（続く）。

エントロピーがわからない

カルノーサイクルからの続き

エントロピーについての熱力学的な導入の論理がすっきりしないと，統計力学の授業が進めにくい。もちろん天下りでボルツマンの原理を導入してしまえばあとは計算だけになる。でも，それでは熱力学との関係もうやむやになりそうだ。

熱力学の第一法則で，$dU=d'Q+d'W=d'Q-pdV$のd'Qの部分も状態量の組み合わせで書けるとありがたい。$p$は示強変数，$V$は示量変数であり，その積がエネルギーの次元を持つ示量変数になっている。使える状態量として示強変数である温度$T$があるので，これに相補的な示量変数で温度との積がエネルギーの次元を持つ状態量をエントロピー$S$として導入して，$d'Q=TdS$とおくことにする。

もしこれができれば，状態量空間中の点を$A$，基準点を$O$として，$S(A)=\int_O^A \frac{d'Q}{T}$は状態量になる。この積分が状態量であるということは，平衡状態Aのみに依存して積分の経路にはよらないはずである。

そこで，カルノーサイクルの断熱過程で実際にこの量を計算してみれば，断熱過程ではエントロピー$S$が一定になる。つまり，カルノーサイクルというのは，エントロピーと温度を2軸とする状態図において，等温線と等エントロピー線に囲まれた長方形領域になる。

ここまでの議論は，準静的過程＝可逆過程について成り立つ話である。不可逆過程だとどうなるのか。肝腎の熱力学の第二法則との関係がついていないわけなのでさらなる検討が必要だ。

クローニッヒ・ペニーモデル（４）

クローニッヒ・ペニーモデル（３）からの続き

これまでの結果を用いて，具体的な数値を代入したグラフを作成する。与えられた$k_n$の値に依存して，$-1 \le \cos(k_n a) \le 1$の条件から，1電子のエネルギー$E$がとりうる範囲についての条件が定まる。このとき，電子が取り得るエネルギー領域を許容帯（allowed band），取りえないエネルギー領域を禁制帯（forbidden band）とよぶ。相互作用のない電子の多粒子系では，これらの離散化された許容帯に電子が充填されることになる。

Mathematicaによって，この様子を確認してみれば次のようになる。

a = 2; b = 0.04; hc = 1973; m = 1.022*10^6; p = 2*1.05017;
v = 2*hc^2*p/(1.022*10^6*a*b)
(hc)^2/(m a^2)
200.001
0.952233

\[Alpha][e_] := Sqrt[m*e]*(a - b)/hc;
\[Beta][e_] := Sqrt[m*(v - e)]*b/hc;
\[Gamma][e_] := Sqrt[m*e]*a/hc;
X[e_] := Cos[\[Alpha][e]] Cosh[\[Beta][ e]] + ((\[Beta][e]/b)^2 - (\[Alpha][e]/(a - b))^2)/(2*\[Alpha][ e]/(a - b)*\[Beta][e]/b) Sin[\[Alpha][e]] Sinh[\[Beta][e]]

Plot[{ X[e], 1, (Cos[\[Gamma][e]] + p*Sin[\[Gamma][e]]/\[Gamma][e]), -1, Cos[\[Gamma][e]]}, {e, 0, 200}, PlotRange -> {-1.5, 3.5}, PlotStyle -> {Red, Gray, Blue, Gray, Orange}]
Table[FindRoot[X[e] == 1, {e, (hc a n )^2/(2 m) }], {n, 1, 5}]
{{e -> 2.95059}, {e -> 37.6031}, {e -> 45.1073}, {e -> 150.412}, {e -> 158.267}}
Table[ FindRoot[X[e] == -1, {e, 0.9*(hc a n )^2/(2 m) }], {n, 1, 6}]
{{e -> 9.40077}, {e -> 16.0036}, {e -> 84.6069}, {e -> 92.3771}, {e -> 242.886}, {e -> 242.886 + 0. I}}

d0 = Plot[{10^6 (e - 2.95059), 10^6 (e - 9.40077), 10^6 (e - 16.0036), 10^6 (e - 37.6031), 10^6 (e - 45.1073), 10^6 (e - 84.6069), 10^6 (e - 92.3771), 10^6 (e - 150.421), 10^6 (e - 158.267)}, {e, 0, (4 Pi)^2}, PlotStyle -> Table[{Gray, Dotted}, 9], PlotRange -> {0, 13}];
f0 = Plot[{0, Pi, 2 Pi, 3 Pi, 4 Pi}, {e, 0, (4 Pi)^2}, PlotStyle -> Table[{Green, Dotted}, 5]];
g0 = Plot[\[Gamma][e], {e, 0, (4 Pi)^2}, PlotStyle -> {Orange, Dashed}];
g1 = Plot[ArcCos[X[e]], {e, 0, Pi^2}];
g2 = Plot[2 Pi - ArcCos[X[e]], {e, Pi^2, (2 Pi)^2}];
g3 = Plot[2 Pi + ArcCos[X[e]], {e, (2 Pi)^2, .98 (3 Pi)^2}];
g4 = Plot[4 Pi - ArcCos[X[e]], {e, .98 (3 Pi)^2, .95 (4 Pi)^2}];

Show[{d0, f0, g0, g1, g2, g3, g4}, PlotRange -> {-1, 15}]

図１：境界条件から決まる（energy-cos(ka)）のグラフ

図２：図１から決まる分散関係（energy-ka）のグラフ

ポテンシャルの幅$b$を格子周期$a$の2%に設定したため，$b \rightarrow 0$の近似がよく当てはまる。そこで，$X(pa)$と$\cos \gamma(e) + P \sin \gamma(e)/\gamma(e)\ $の誤差は1%以下となり，グラフ上では区別できなくなっている。

なお，$e=\frac{(\hbar k)^2}{2m}=\frac{(\hbar c)^2}{2mc^2 a^2}\cdot (ka)^2 = 0.952233 (ka)^2\ $であることに注意する。

クローニッヒ・ペニーモデル（３）

クローニッヒ・ペニーモデル（２）からの続き

1次元ポテンシャルに周期性があるときに，ブロッホの定理から$\psi(x)=e^{ikx}\varphi(x)$と表わせて，$\psi(x+a)=e^{ik(x+a)}\varphi(x+a)=e^{ika} e^{ikx}\varphi(x)=e^{ika}\psi(x)$が成り立つ。このときの波動関数は運動量演算子の固有状態なのだろうか？違います。前回やったように，このハミルトニアンは有限の並進操作に対して不変だけれど，運動量に対応する無限小並進操作については不変ではないから。

ところで，この長さ$L=N a$の1次元周期ポテンシャルモデルの両端を同一視する周期境界条件をつけると（$N$はポテンシャルステップの数＝原子数，$a$はポテンシャルの周期＝原子間隔），$\psi(L)=\psi(0) \quad \psi(L)=e^{i k a \cdot N}\psi(0) \quad \therefore e^{i k a N}=1$

これから$k$に対する条件，$k_n = \frac{2\pi n}{a N}\quad (n=0,\pm 1, \pm 2 \cdots)\ $が得られる。$k_n$は量子数 $n$ で特徴づけられるこの状態の波数という意味をもつ。

前回得られた境界条件は，系のエネルギーを$E$，ポテンシャルの深さと幅を$V_0, b$，ポテンシャル周期を$a$として，$p=\frac{\sqrt{2mE}}{\hbar}$，$q=\frac{\sqrt{2m(V_0-E)}}{\hbar}$とおくと，$\cos k a  =  \cos p(a-b) \cosh qb + \frac{q^2-p^2}{2 p q} \sin p(a-b) \sinh qb$ である。これは，与えられた$k = k_n$に対して，系のエネルギーを決定する式になる。

(1) $b \rightarrow 0$ の極限では$p_n=k_n$となり，$E_n=\frac{\hbar k_n^2}{2m} = \frac{2 \hbar^2 \pi^2 n^2}{m a^2 N^2}$となる。

(2) 次に，$V_0 b$を一定に保ちながら，$b \rightarrow 0,\ V_0 \rightarrow \infty$とするδ関数型極限を考える。このとき，$\sinh qb \rightarrow qb$であり，右辺第2項は，$\frac{(q^2-p^2)ba}{2} \frac{\sin p(a-b)}{p a}$となる。最終的に，$\cos k a  =  \cos p a + \frac{m c^2 V_0 b a}{(\hbar c)^2} \frac{\sin pa}{pa}$ という近似式が得られる。

例えば，$a=2$ Å，$b=0.04$ Å，$mc^2 = 0.511 \times 10^6$ eV，$\hbar c =1973$ eV Å， $V_0=100$ eVとすると， 無次元のポテンシャル強度パラメータは，$\frac{m c^2 V_0 b a}{(\hbar c)^2}=1.05$となる。

日本國憲法前文

憲法記念日だが，ウクライナへのロシアの侵略戦争が続いている。先の見えない円安で疲弊した我々の心の隙に，右翼デマゴーグ達の好戦的なイデオロギーが陽に陰に染み込んでいく。こんなときは，日本國憲法の前文を写経して心を鎮めるしかない。若者は「最高法規の意志～ 憲法の本質と改正の動向 ～」をチラ見する方が役に立つかも。

　　　日本國憲法 

　日本國民は、正當に選擧された國會における代表者を通じて行動し、われらとわれらの子孫のために、諸國民との協和による成果と、わが國全土にわたつて自由のもたらす惠澤を確保し、政府の行爲によつて再び戰箏の惨禍が起ることのないやうにすることを決意し、ここに主權が國民に存することを宣言し、この憲法を確定する。そもそも國政は、國民の嚴肅な信託によるものであつて、その權威は國民に由來し、その權力は國民の代表者がこれを行使し、その福利は國民がこれを享受する。これは人類普遍の原理であり、この憲法は、かかる原理に基くものである。われらは、これに反する一切の憲法、法令及び詔勅を排除する。

　日本國民は、恒久の平和を念願し、人間相互の關係を支配する崇高な理想を深く自覺するのであつて、平和を愛する諸國民の公正と信義に信頼して、われらの安全と生存を保持しようと決意した。われらは、平和を維持し、專制と隷從、壓迫と偏狹を地上から永遠に除去しようと努めてゐる國際社會において、名譽ある地位を占めたいと思ふ。われらは、全世界の國民が、ひとしく恐怖と缺乏から免かれ、平和のうちに生存する權利を有することを確認する。

　われらは、いづれの國家も、自國のことのみに專念して他國を無視してはならないのであつて、政治道徳の法則は、普遍的なものであり、この法則に從ふことは、自國の主權を維持し、他國と對等關係に立たうとする各國の責務であると信ずる。

　日本國民は、國家の名譽にかけ、全力をあげてこの崇高な理想と目的を達成することを誓ふ。

クローニッヒ・ペニーモデル（２）

クローニッヒ・ペニーモデル（１）からの続き

前回，周期ポテンシャル中で正のエネルギーを持った電子の運動について考えた。自由電子とはいえ，金属中に束縛されているのだからポテンシャルの上端に対して負のエネルギーを持った電子を考えなければならなかった。

上智大学名誉教授の清水清孝さん（元有馬研）が，電子情報通信学会の知識ベース知識の森の12群5編量子力学・電子物理・相対論を執筆していて，そこにバンド理論入門についての話もあった。

そこで，前回のポテンシャルの符号を $-V_0 \rightarrow V_0$としたモデルで$0<E<V_0$の場合を考える。すなわち，ポテンシャルは，つぎの形を繰り返したものになる。
領域 Ⅰ：$V(x) =  \ 0 \ \cdots \ (0 < x < a-b)\ $
領域 Ⅱ：$V(x) = V_0 \ \cdots \ (-b < x < 0)$

自由電子のエネルギーは正だが，ポテンシャルの高さより小さいので，領域Iでは解は平面波，領域IIでは指数関数になる。$p=\frac{\sqrt{2mE}}{\hbar}, q=\frac{\sqrt{2m(V_0-E)}}{\hbar}$とおいて一般解は，

領域 Ⅰ：$\psi_{\rm I}(x) = A e^{i p x} + A' e^{-i p x} \ \cdots \ (0 < x < a-b)\ $

領域 Ⅱ：$\psi_{\rm II}(x) =  B e^{q x} + B' e^{- q x} \ \cdots \ (-b < x < 0)\ $

前回のように波動関数とその導関数が，領域Iと領域IIの境界で連続である条件から，

$\psi_{\rm I}(0) = \psi_{\rm II}(0);\ \psi_{\rm I}'(0) = \psi_{\rm II}'(0)$，$\psi_{\rm I}(a-b) = e^{i k a }\psi_{\rm II}(-b);\ \psi_{\rm I}'(a-b) = e^{i k a }\psi_{\rm II}'(-b)$

$A+A'=B+B'$
$i p(A-A')=q(B-B')$
$A\ e^{ip(a-b)}+A'e^{-ip(a-b)}=(B\ e^{-qb}+B'e^{qb})e^{ika}$
$i pA\ e^{ip(a-b)}-i pA'e^{-ip(a-b)}=(qB\ e^{-qb}-qB'e^{qb})e^{ika}$

ここで，$\alpha = p(a-b)$，$\beta = q b$と置くと，

$\left( \begin{array}{cccc} 1 & 1 & -1 & -1 \\ i p & -i p & -q & q \\ e^{i\alpha} & e^{-i\alpha} & -e^{-\beta}e^{ika} & -e^{\beta}e^{ika} \\ i p e^{i\alpha} & - i p e^{-i\alpha} & -q e^{-\beta}e^{ika} & q e^{\beta}e^{ika} \end{array} \right)  \left( \begin{array}{c} A\\ A' \\ B \\ B' \end{array} \right) = 0$

4元連立方程式が自明でない解を持つためには，行列式が0であり，列の加減により等価な行列式を求める（3行目から$e^{ika}$を，2列目から$i$を外にくくり出した）。

$i e^{ika} \cdot \left| \begin{array}{cccc} 1 & 0 & -1 & 0 \\ 0 & p & 0 & q \\ \cos \alpha \ e^{-ika} & \sin \alpha \  e^{-ika}& -\cosh \beta & - \sinh \beta \\ -p\ \sin \alpha & p\ \cos \alpha & q\ \sinh \beta \ e^{ika} & q\ \cosh \beta \ e^{ika} \end{array} \right|  = 0$

これから4x4行列式の値$|M|$を計算する。

$|M|=  p\ q\ e^{i k a} \left | \begin{array}{cc} -\cosh \beta  & - \sinh \beta  \\ \sinh \beta & \cosh \beta \end{array} \right | + q \left | \begin{array}{cc}  \sin \alpha e^{-i k a } & - \cosh \beta \\  p \cos \alpha & q \sinh \beta e^{ika} \end{array} \right |$

$\qquad \quad - q\ p\ e^{-i k a}\ \left | \begin{array}{cc} \cos \alpha  &  \sin \alpha \\ - \sin \alpha & \cos \alpha  \end{array} \right | - p \left | \begin{array}{cc} - \sinh \beta  & \cos \alpha e^{-ika} \\ q \cosh \beta e^{ika} & -p \sin \alpha \end{array} \right |$

$\quad = -2 p q \cos ka + q( q \sin \alpha \sinh \beta + p \cos \alpha \cosh \beta ) -p (p \sin \alpha \sinh \beta -q \cos \alpha \cosh \beta )$

$\quad = -2 p q \cos ka +2 p q \cos \alpha \cosh \beta +(q^2-p^2) \sin \beta \sinh \beta = 0$

最終的に得られる関係式は，次のとおりである。

$\cos ka  = \cos \alpha \cosh \beta -\frac{p^2-q^2}{2 p q} \sin \alpha \sinh \beta$

$\qquad = \cos p(a-b) \cosh qb -\frac{p^2-q^2}{2 p q} \sin p(a-b) \sinh qb$

（付）この式は，$E>V_0$の場合は次のように変形できるので，このまま使うことができる。$q=\frac{\sqrt{2m(V_0-E)}}{\hbar} \rightarrow i \frac{\sqrt{2m(E-V_0)}}{\hbar} \equiv i\bar{q}\ $として，$\sinh ix = i \sin x,\ \cosh ix = \cos x$を用いると，

$\cos ka  =  \cos p(a-b) \cos \bar{q}b -\frac{p^2 + \bar{q}^2}{2 p \bar{q}} \sin p(a-b) \sin \bar{q} b\ $となる。これは（１）で得られた式と同等なもの。

図：クローニッヒ・ペニーモデル

クローニッヒ・ペニーモデル（１）

クローニッヒ・ペニーモデルは周期性を持った１次元井戸型ポテンシャルのモデルであり，結晶のバンド構造の定性的な特徴を説明することができる。

図：クローニッヒ・ペニーモデルの設定（Wikipediaより引用）

電子の質量を$m$として，1粒子の1次元シュレーディンガー方程式は，$-\frac{\hbar^2}{2m}\frac{d^2}{dx^2}\psi(x) +V(x) \psi(x) = E \psi(x)$である。ポテンシャルは周期性を持ち，$n$を整数として$V(x+n a)=V(x)\ $であり，ブロッホの定理から波動関数は，$\psi(x) = e^{i k x} \varphi(x)$ かつ $\varphi(x + n a)=\varphi(x)$を満たす。

ポテンシャルは，つぎの形を繰り返したものになる。

領域 Ⅰ：$V(x) =  \ 0 \ \cdots \ (0 < x < a-b)\ $

領域 Ⅱ：$V(x) = -V_0 \ \cdots \ (-b < x < 0)$

波動関数の周期性からは，$\psi(x+a) = e^{i k (x+a)}\varphi(x+a) = e^{i k (x+a)}\varphi(x) = e^{i k a} \psi(x) \ $が成り立つ。その導関数は$\psi'(x+a) = e^{i k a}\psi '(x)$となる。なお，$\psi'(x) = i k e^{i k x}\varphi(x) + e^{i k x}\varphi'(x) = i k \psi(x) + e^{i k x}\varphi'(x)$である。

自由電子のエネルギーは正であり，ポテンシャルが定数であるため，どちらの領域でも解は平面波になる。$p=\frac{\sqrt{2mE}}{\hbar}, q=\frac{\sqrt{2m(E+V_0)}}{\hbar}$として一般解は，

領域 Ⅰ：$\psi_{\rm I}(x) = A e^{i p x} + A' e^{-i p x} \ \cdots \ (0 < x < a-b)\ $

$\qquad \quad = e^{i k x} (A e^{i (p-k) x} + A' e^{-i (p+k) x} )$

領域 Ⅱ：$\psi_{\rm II}(x) =  B e^{i q x} + B' e^{-i q x} \ \cdots \ (-b < x < 0)\ $

$\qquad \quad = e^{i k x} (B e^{i (q-k) x} + B' e^{-i (q+k) x})$

波動関数とその導関数が，領域Iと領域IIの境界で連続であるという条件を書く。
$\psi_{\rm I}(0) = \psi_{\rm II}(0);\ \psi_{\rm I}'(0) = \psi_{\rm II}'(0)$，$\psi_{\rm I}(a-b) = e^{i k a }\psi_{\rm II}(-b);\ \psi_{\rm I}'(a-b) = e^{i k a }\psi_{\rm II}'(-b)$

$A+A'=B+B'$
$p(A-A')=q(B-B')$
$A\ e^{ip(a-b)}+A'e^{-ip(a-b)}=B\ e^{-iqb+ika}+B'e^{iqb+ika}$
$pA\ e^{ip(a-b)}-pA'e^{-ip(a-b)}=qB\ e^{-iqb+ika}-qB'e^{iqb+ika}$

ここで，$\alpha = p(a-b)$，$\beta = q b$と置く。

$\left( \begin{array}{cccc} 1 & 1 & -1 & -1 \\ p & -p & -q & q \\ e^{i\alpha} & e^{-i\alpha} & -e^{-i\beta}e^{ika} & -e^{i\beta}e^{ika} \\ p e^{i\alpha} & - p e^{-i\alpha} & -q e^{-i\beta}e^{ika} & q e^{i\beta}e^{ika} \end{array} \right)  \left( \begin{array}{c} A\\ A' \\ B \\ B' \end{array} \right) = 0$

波動関数の係数に対するこの4元連立方程式が自明でない解を持つためには，行列式が0でなければならない。この条件をつかって，2列と1列の和と差，4列と3列の和と差から等価な行列式を求めれば，次のようになる（3行目から$e^{ika}$を外にくくり出した）。

$e^{ika} \cdot \left| \begin{array}{cccc} 1 & 0 & -1 & 0 \\ 0 & p & 0 & q \\ \cos \alpha \ e^{-ika} & i \sin \alpha \  e^{-ika}& -\cos \beta & -i \sin \beta \\ i p\ \sin \alpha & p\ \cos \alpha & i q\ \sin \beta \ e^{ika} & q\ \cos \beta \ e^{ika} \end{array} \right|  = 0$

これから4x4行列式の値$|M|$を計算する。

$|M|= p\ q\ e^{i k a} \left | \begin{array}{cc} -\cos \beta  & - i \sin \beta  \\ i \sin \beta & \cos \beta \end{array} \right | + q \left | \begin{array}{cc} i \sin \alpha e^{-i k a } & - \cos \beta \\ p \cos \alpha & i q \sin \beta e^{ika} \end{array} \right |$

$\qquad \quad - q\ p\ e^{-i k a}\ \left | \begin{array}{cc} \cos \alpha  & i \sin \alpha \\ i \sin \alpha & \cos \alpha  \end{array} \right | - p \left | \begin{array}{cc} -i \sin \beta  & \cos \alpha e^{-ika} \\ q \cos \beta e^{ika} & i p \sin \alpha \end{array} \right |$

$\quad = -2 p q \cos ka + q( -q \sin \alpha \sin \beta + p \cos \alpha \cos \beta ) -p (p \sin \alpha \sin \beta -q \cos \alpha \cos \beta )$

$\quad = -2 p q \cos ka +2 p q \cos \alpha \cos \beta -(p^2+q^2) \sin \beta \sin \beta = 0$

最終的に得られる関係式は，次のとおりである。

$\cos ka  = \cos \alpha \cos \beta -\frac{p^2+q^2}{2 p q} \sin \alpha \sin \beta$

$\qquad = \cos p(a-b) \cos qb -\frac{p^2+q^2}{2 p q} \sin p(a-b) \sin qb$

4行4列の行列式は，Mathematicaを使えば手軽に計算できるのだけれど，手計算でもなんとかなる場合があるということを学ぶ。