博士前期課程講義 平成29年度 13回目

2017年7月11日。

今日の内容は、ランダム・ウォーク、レビー・フライト、DLA。

ランダム・ウォークについては、ランダムな高分子鎖の広がりないしはコンフォーメーションの模型として(習っている人は)習っているでしょう。回転半径(慣性半径かも知れません)が重合度nに比例する話です。第nステップでの位置をR(n)として、<R2(n)>-<R(n)>2=na2(aは格子間隔)を示した後、粒子の存在確率u(r,t)が連続極限(長時間経過後)に拡散方程式になることをまず示す。その後、空間フーリエ変換を使って拡散方程式をとく。確率分布u(r,t)がわかったので、時刻tに置けるr2の平均が計算でき、それがtに比例することがわかる。もちろん、これは<R2(n)>-<R(n)>2=na2と同じ意味の微視的な式。これから、ランダム・ウォークの軌跡のフラクタル次元はdw=2(空間の次元dによらない)ことが導出できる。

レビー・フライトは、ランダム・ウォークでステップあたりの移動距離が一定となっていた正弦を外し、べき分布としたもの。1<fw≦2ステップの長さsがSより小さい確率をP(S)=(S/a)fwとする(もちろんS≧a)。fw=1の場合は弾道模型なることをコメントし、fw>2nの場合はステップ長一定の場合のと透過である(「中心極限定理により」という難解なもの、とのこと)のコメントも。支配方程式をランダム・ウォークの場合と同じように連続近似して求め、それをフーリエ変換を使って解く。ただし、1次元の場合に限定。ポイントは、最終的に∂u~(q,t)/∂t ~ -qfw

u(q,t)を通じて、u~(q,t)=u~(q,0)exp(-Ddiffqfwt)となること。ランダム・ウォークの場合はu~(q,t)=u~(q,0)exp(-Ddiffq2t)であった。これからdw=fwが導かれる(積分してn ~ Rfwを示すのは困難)。異常拡散(平均二乗変位∝tβ)の話をして次へ。

DLAの実験例は既に挙げてある。少し補足をした後に、DLAのアルゴリズムラプラス方程式に基づいた解析について話。その後、DLAのフラクタル次元の定式を与えたToyoki and Honda (Phys. Rev. Lett.)の理論の紹介をして終了。

一杯まで掛った(30秒未満しか時間が余らなかった)。途中で休憩を入れて、測定中のサンプルのセッティングをしたかったが、とても無理。フーリエ変換を使って拡散方程式を解くことは、他の授業でも需要があるようで、講義終了後に「教えて下さい」という院生あり。

講義2 平成29年度 第9回目

2017年7月7日。

後半の第2回目。前回の最後に定義した化学ポテンシャルについて。節タイトルは、化学ポテンシャルの性質となっているが、その前半は相平衡の条件。化学ポテンシャルが熱力学ポテンシャルをモル数で微分したものと等しいことは、前回は単なる関係式だったが、今日はその意味を。私は、等温定積系の二相平衡をヘルムホルツエネルギー最小から導出することをやっている。二相の圧力が等しいという相平衡の条件の一つが、二相界面をどの位置に持って来たら圧力の釣り合いが実現できるかという、力学的平衡の条件であることがわかるから。物質移動に関する平衡の条件のみが出てくる、等温定圧系でギブスエネルギーを最小にする方法では、力学的平衡の条件が二相の圧力が等しいことを「同様に」で説明しづらい。それならば、エネルギーの移動に関する平衡の条件も出てくる、内部エネルギー最小でやればいいでないか、と言うことになるが、TV一定はイメージが簡単だが、SV一定は少し困難がある。従って、エネルギーの移動に関する平衡の条件が二相の温度が等しいということは、同様で済ます。元も、孤立系のエントロピー最大ならば、困難はないが、何のために熱力学ポテンシャルを定義したのか、と言うことになってしまう。

ギブス関係式、オーラーの関係式、ギブス・デュエム関係式については、内部エネルギーUがS,V,nの一次の同時式であることからスタート。教科書のように、等温定圧でのギブス・デュエム関係式Σndμ=0では終らない。単成分についてdμ=-SmdT+VmdPで終り、次の理想気体の化学ポテンシャルへつなげる。添え字のmはモル量(モルエントロピー、モル体積)であることを表す。

理想気体の化学ポテンシャルは、ギブス・デュエム関係式dμ=-SmdT+VmdPから出てくる(∂μ/∂P)T=Vm積分するのみ。するとμ(T,P2)=μ(T,P1)+RTln(P2/P1)は、即座に出て来る。理想混合気体のμi(T,P,{xi})=μi*(T,P)+RT ln xiは、G(T,P,{xi})=G*(T,P)+ΔGmix; G*=Σniμi*(T,P), ΔGmix=RTΣniln xiから。

質量作用の法則のところは、化学量論係数を定義し、反応進行度を導入し、それを用いてギブスエネルギー最小を論じるだけ。化学ポテンシャルの式まで変形したら、後はμi(T,P,{xi})=μi*(T,P)+RT ln xiを適用したら、質量作用の法則になる。

6分オーバーしてしまった。指がこわばって板書に難があったことは確か。「エントロピー、増えようが減ろうが知ったことかい」と言う川柳(?)を紹介したのは余分だったかもしれないが、印象に残って欲しいとは思う。

○合事前審査

2017年7月6日。

博士後期過程の学生の指導資格のことを〇合という。

研究部で統一された基準ほぼギリギリで出された方と私の2名が〇合「事前審査」に通った。まだ「事前審査」の段階で、事務に提出した書類を学科に戻して審査(いわば予備審査かな)した後、研究部での審査をもって完了となる。

私の所属している学科では、申請調書の提出前に、いわば事前審査を行っている。かつて「教授に該当する」というのが〇合の前提だった。それを引きづっていて、私の所属学科では事前審査を通った前例はなく、〇合准教授が他学科から移ってきたものを除けば、〇合の合格者はゼロであった。正確には、学部と研究部が別組織だった頃に、研究部の審査で〇合にパスした准教授はいた - A学科の材料系は材料系の研究部に所属、A学科の情報系は情報系の研究部に所属の頃の話・・・具体的には、情報系にはレター=速報の概念がなく、例えばPhysical Review Lettersでもショートペーパー扱いになってしまうが、そういうことをせずに分野標準で評価していた頃の話。

上の例だと、情報系の基準を杓子定規に適用して、レター(短報でなく、フルペーパー(full length paperと言った方がいいかも)の速報も含む)を排除した上で論文数をカウントすることをやって、他学科ならばらくらく〇合合格の業績の教員が〇合に通らない不公平がおきる。研究部の上部の方が不公平さをなくしたい意向を示している、と言う背景が生じた。

講義2 平成29年度 第8回目

2017年7月5日。

講義2の8回目は、後半の1回目。第二法則で熱力学の本質は終っている。後半は自由エネルギーを用いて第二法則を書き換えること。そして数学的に扱いやすい形にして、様々な応用へ。本質的には新概念は、ない。

熱力学ポテンシャルと言う語を必ず紹介している。一つは、自由エネルギー最小が経項状態を与えるという、力学的釣り合い=力学的ポテンシャル最小と同様な意味があること。もう一つは計から取り出せる仕事の最大が熱力学ポテンシャルの差になるという、力学的な系の最大仕事と同様な意味のあること。後者を先に説明し、前者を孤立系におけるΔS≧0と同じように説明した。

熱力学の関係式の説明。どの変数を制御した実験を行い、どの量をどの変数の関数としてそくていするか、ということを毎回言うようにしている。

最後に開放系において系に物質が流入したときの仕事の表現と言う観点で科学ポテンシャルを導入。熱力学ポテンシャルをモル数で偏微分する定義については、次節(次の時間)にやるところと関係している。次回に詳細にやります。待てない人は、教科書に書いてある程度のところを自分で読んでおいて下さい。

 

続きを読む

博士前期課程講義 平成29年度 12回目

2017年7月4日。

フラクタル次元の解析法で既にやっているものお除いて紹介。スケール変換法は、ダブって説明。視野拡大法を説明するために必要。視野拡大法の「発展形」として回転半径法。高分子については、まさにぞれそのものが出てくる。DLAなどの成長するパターンについては強力。一つのバターンがあって、それに回転半径法を適用することはできない。無理やり適用した形の例が載っている教科書の例は紹介。密度相関関数法は、一部は済み。フラクタルバターン以外についても、相関関数は計算して見るといい。フラクラル以外に適用できるものとして、アイランドのサイズ分布の解析がある。

今日は、体がだるい。胸が苦しい。まさに梅雨の天候だから、だるいのかも知れない。気の利いたことは、一切やらずに終了。

講義2 平成29年度 第7回目

2017年6月30日。

今日で6月も終りですね。熱力学第二法則も今日で終わりです。今日は、エントロピーの計算を行います。その後に第二法則に付随する事項をやります。前半のレポート課題も配布する。嫌なことも通知しなければないらい・・・本年度は、出席の不足している学生にポートフォリオを提出させて(試験を受けるに値するだけの学習を尾行っているかチェックして)、試験の受験資格の認定を行うことはしない。もう、一杯なので、他の職務に支障が出てしまうから。

まずは、第一法則のところの演習問題の回収。

dS=d'qrev/Tを印象付ける話で、可逆的食物摂取の話。冗談ぽい話がツボにはまって、理解を深めたり、記憶に残ったりすることがある。生身の人間が講義する意義。今日の私の体温は36.5℃でした。いつもは36.1-2℃なので、膠原病の関節症状があって、熱が少し高いのかもしれません。これは、十数年前の学生が言ったことですんので、冗談で言ったのかもしれません。今日は36.5℃で〇〇カロリーを摂取したんで、これだけのエントロピーを増やしました。普段だと、もう少し体温が低いので同じカロリー摂取でもエントロピー増加は少し多いことになります。ポイントは、「可逆」です。可逆的にカロリーを摂取した場合はそうなります。可逆的食物摂取と言う話に進むと、飛んでもないことになります・・・。摂取したカロリー排出すれば元に戻ったことになるのでしょうか。あなた自身、系は元に戻ったのかもしれませんが、環境は元に戻りませんね。

そんな話をした後に、エントロピーの計算を。まずは、理想気体の等温膨張。次いで温度変化に伴うエントロピー変化。相転移に伴うエントロピー変化は、相転移熱の説明が薄かったので、それを補うように話した。第一法則のところの順番だと、次は反応エントロピーとなるが、それは別の章で詳しくやることになる。混合のエントロピーは、少し仰々しく。理想気体の混合の場合は、等温膨張と同じだが、仰々しく「混合のエントロピー」といって、モル分率を使った式を書くのがためになる。つまり、理想気体でなくても、分子の(大きさ)相互作用がなければ同じものが出てくる。

エントロピーの分子論的な意味を説明し、第三法則について話した後、マクスウェルデーモンの話もする。皆さん、神様の話し、待っていたでしょう。反省してもうエントロピーを増やすことはしませんので、神様、私が無駄に増やしたエントロピーを元に戻して下さい、の話。可逆的カロリー摂取/食物摂取と同じで、印象付けにしかならない人もいれば、高度な話として捉えられる学生もいるかもしれない。あ、その前に「可逆」について。医学などでは、例えば関節が変形してもそれが元に戻せたり、胃に穴が空いても元に戻せたりする場合に「可逆」、そうでない場合に「不可逆」っていいますね。体に不可逆変化が起きてしまうと、今後の人生が大変になりますよ。環境に変化を残したくないと言って、体が不可逆な変化をしてしまわないように。さて、神様の話。物理学では、神様でなくてデーモンって言葉を使います。まじめに議論されてきた話です。気体の自由膨張を考えましょう。不可逆変化ですね。ところが、気体が入っていた容器と真空だった容器をつなぐ管を分子1個が通る程度の細さにし、分子が右向きに動いていたらシャッターを閉め、左向きに動いていたらシャッターを開けるという操作をしたとしましょう。元の状態に戻ってしまいますね。このような操作をするデーモンをマクスウェルデーモンって呼びます。昔はデーモンだったかもしれませんが、今は計算機制御で分子を識別してシャッターを開け閉めする装置はできそうですね。デーモンということで計算機の電力は問題にならない程度に理想化したとしましょう。その場合、自由膨張を行う前の状態に戻ったのですから、系のエントロピーは減少しますね。最新の研究によると。計算機が情報処理を行ったことに伴うエントロピー変化と系のエントロピー変化を合わせると、トータルではエントロピーが減少することはないということです。管とシャッターの系にとっては、計算機は外界だったのですね。

博士前期課程講義 平成29年度 11回目

2017年6月28日。

昨日の講義には、難解なハウスドルフ次元の話が含まれる。講義1の採点を今月中に終らせようと悪戦苦闘しているなかで、どこまでできるか。正確にやりつつ、ある程度はわかった気にもさせる必要がある。

ベクトル空間の次元から入り、位相次元を説明。その後、ハウスドルフ測度などの数学の準備を行った後にハウスドルフ次元へ。相似次元を再定義。最後にボックスカウント次元の定義。ハウスドルフ次元が従来次元を含むことは説明した。ボックスカウント次元がハウスドルフ次元の簡略版の意味を持っていることも説明した。

実は、本日は代休。臀部皮膚膿瘍/粉瘤の切除。