はじめまして

いらっしゃいませこんにちは、期待値観測するには確率の問題に超絶疎い、量子きのこです



好奇心はあるのにメンタル不足でソレ系の職業に就けなかったので
何か出来ないかと思ったわけじゃないですが
結果的に物理や数学とかの面白さを伝えられたらいいなー
なブログになりました。

このブログをきっかけに、一人でも多くの人が楽しんで物理や数学を学んでくれるとサイワイです。ヽﾟーﾟﾉ




サイエンスファンタジーなアニメの感想も書きます。
あまり伝わらないようなことを毎日ぼちぼちやっています。
伝わらない・人を選ぶ内容なので、ついカッとなって今は公開しています。




 「ブログ村」というブログ紹介サイトにいつもお世話になっております。
ﾖｶﾀ！と思ったら押してネ！
 
 
 




※同姓同名の同一人物「量子きのこ」が、pixivにも生息してます。
よかったら覗いてやってください

調和振動子の昇降演算子

量子力学の本を返したのでこのチャンスを活かして、調和振動子のほうの昇降演算子の、
空で書けるようになった部分を書きます。

まず、調和振動子ポテンシャルのハミルトニアンを無次元化します。
ハミルトニアンの次元はエネルギーなので、
同じくエネルギーの次元を持つ



の2倍でハミルトニアン



を割ります。



hにバーがついたのはプランク定数を2πで割ったディラック定数
kはばね定数、mは質量、xは位置、pは運動量、ωは角振動数で



の関係があるので、無次元化したハミルトニアンは以下のようになります。



そしてこれをガウス素数っぽく半分に因数分解します。


この片方、どっちか忘れましたがプラスがついたほうをaと置くと

マイナスがついたほうはaの複素共役(というかたぶんエルミート共役)a†と書くことができます。


†をどっちのaにつける慣習かは忘れましたが
運動量pにプラスiがついたほうを下降演算子、マイナスiがついたほうを上昇演算子と呼び

波動関数に上昇演算子を左からかけると、量子数1個分だけ固有値がアップした波動関数を
波動関数に下降演算子を左からかけると、量子数1個分だけ固有値がダウンした波動関数を
算出することができます。

ではこの芋づるの端っこはといいますと
一番下の波動関数に下降演算子をかけると恒等的にゼロになることが知られているため、それを利用して微分方程式を解くことで算出できます。


つまり具体的にやってみますと

物理量xとpを微分演算子にして

このように解くことができます。(Cは任意の定数です)

粘性抵抗、慣性抵抗、落下運動、両方、運動方程式4

前回までは、粘性抵抗と慣性抵抗両方が加わった落下運動の運動方程式の解について
粘性抵抗と慣性抵抗片方ずつに係数を制限して、元々どちらかしか考慮していない運動方程式の解と同じになるかを検証していた。また、それぞれの終端速度についても議論していた。


今回はいよいよ、粘性抵抗と慣性抵抗両方を含んだ状態の終端速度について考えてみる。


微分方程式の解は以下のようになるのだったが

    

このとき、時間t→∞の極限で終端速度vはどうなるだろうか。

cothの中身が大きくなると1に収束するので
 

分母分子に をかけて
 

これは、元の微分方程式である
mv'=mg-c1v-c2v^2の加速度(力)v'をv'=0としてvについて解いた2次方程式の解そのものである。
そのうえ、数回前に としたDがまさに判別式であったこともわかるだろう。

また、すべての変数、関数、パラメータが正の実数と仮定すると

2次方程式の判別式Dのルートは必ず正の実数であり
Dのルート√(c1^2+4mgc2)は必ずc1より大きいため、v=-c1-√(c1^2+4mgc2)/(2c2)の解はありえないこともわかる。

粘性抵抗、慣性抵抗、落下運動、両方、運動方程式3

前回はc1≠0、c2=0の条件で解いたが、今度は逆にc1=0、c2≠0の条件で解いてみよう。



となるが、元の運動方程式

 

も解いてみると

 

部分分数分解を用いて

 

(A1+A2)v=0

からA1=-A2

 

 


初期条件t=0でv=0を与えるとG2=-1となって


ここで、式を整理してvの式にする
両辺にv+√(mg/c2)をかけると
 
ここで、右辺の分母分子にをかけると

これは双曲線関数ハイパボリックタンジェントなので


を得、見事に両者は一致する。


また、終端速度については

t→∞でv=√(mg/c2)になるが、

これについても元の運動方程式

mv'=mg-c2v^2のv'=0の条件で解くことで、単なる代数方程式として
終端速度v=√(mg/c2)を簡単に得ることができる。

粘性抵抗、慣性抵抗、落下運動、両方、運動方程式2

前回、粘性抵抗と慣性抵抗両方を受ける落下運動について運動方程式を解いた。

以下のような結果になったが
 


G=m/c2、B=c1/c2、A=mg/c2

このようにB,A,G,Dを定義したことを踏まえて、元の格好に戻すと

 

このような姿になる。

ここで、c1やc2をそれぞれゼロとしたとき、

元の微分方程式の解と一致しているかどうかを確かめてみよう。

c1≠0、c2=0の場合は

 

となり、

運動方程式

 

の解と一致しているはずだ。

早速解いてみると

 

 

 

 

初期条件としてt=0でv=0を与えると任意の定数G2=-mg/c1となるので

 

 

 

一方、 こちらはcothの中身を展開し

大きい分母の分母分子にexp(+)-exp(-)をかけると


となって見事一致する。

ちなみに、終端速度もちゃんと一致していることを確かめておこう。

t→∞でのvはv→mg/c1に収束するが

これは実は、元の運動方程式

mv'=mg-c1vのv'=0の条件で解けば、微分方程式ではなくただの代数方程式として
0=mg-c1v
c1v=mg
v=mg/c1

と、簡単に解けてしまう。

粘性抵抗、慣性抵抗、落下運動、両方、運動方程式

質量をm、重力加速度をg、速度v、時間t
速度の1乗に比例する空気抵抗の係数をc1、速度の2乗に比例する空気抵抗の係数をc2とすると

運動方程式は以下のように書ける。
速度は高さy方向の下向きをプラスに定義する。

 

ここで、v^2の係数を1にしたいので、両辺を-c2で割ると

 

ここで、式の簡単化のためG=m/c2、B=c1/c2、A=mg/c2とすると

 

dvの分母に右辺を持ってくると以下のようになる。

 

両辺を積分してvの式にすると、速度vについての運動が解けるのだが、vの2乗が邪魔でうまく積分できないので、部分分数分解を用いる。


 

という恒等式を作り、A1、A2、Dがいくつになるのか定めることにする。

まず、 になるようにDを定めたい。

 

なので、 であることがわかる。なぜいきなりDという文字を仮定したのかはあとでわかるが、これは2次方程式の判別式Dの意味である。


次にA1とA2を定める。
 という恒等式からA1とA2を定めたいので

v(A1+A2)=0
からA1=-A2が導かれるので
B/2(A1+A2)もゼロである。

D(A1-A2)=1にA1=-A2を大入して、

 を得る。


そうすると元の微分方程式は

 
なので積分ができ

 

となる。G2は任意の積分定数。

対数同士の引き算は対数の中の割り算として引っ込めることができるため

 

両辺に2Dを掛け算して
 
(ここでG2はまだ任意の定数なので、2DG2を新たなG2として上書きする)


両辺の指数をとって

 
(ここでもexp(G2)を新たなG2として上書きしている)

左辺の分母分子に2を掛け算してやると
 


初期条件
t=0でv=0を与えて任意の定数G2を定める。

 

微分方程式の解は
 

となる。


これをvの式に直す。
両辺に(B+2D)と(2v+B+2D)を掛け算する。

 

展開すると


全部左辺に移行し、vを含むか含まないかで整理すると






  

 

ここで1/(2v)の式に着目してみることにする。




第1項は約分できるので
 

問題は第2項だが

の分母分子にexp(x)を掛け算すると

となって、これは双曲線関数ハイパボリックタンジェントtanhの逆数ハイパボリックコタンジェントcothであることがわかる。
 





ここでようやく2vの式に戻すと

トンネル障壁にパルス波をぶち込む

Excelでマイクラ風フルル

2次元井戸型ポテンシャル。解析解につきこのような形のポテンシャルに

4次元空間＋1次元時間状態図の熱力学「カルノーサイクル」

3Dカルノーサイクル状態図(Excelおもちゃ)

シークバーで見れるバージョンはpixivにあります。


べき乗関数の宝庫である熱力学の状態図、
両対数方眼紙でべき乗関数は直線になりますが、
変数が3つ以上あるので、とりあえずエントロピーを除いて、3D対数ブロックにしたらどうなるか試してみました。

下にDL用Excelファイル(マクロなし)へのリンクを貼りましたので、自粛期間中に暇だったら遊んでみてください。(いつもながら対象年齢高めですみません)
シート保護の関係でサンプルgif動画と異なる表示がされています。




DLしたら、「編集を可能にする」にしてからご使用ください。

3DカルノーサイクルExcelファイル（ダウンロード）

身勝手の極意～井戸型ポテンシャル編～

昨日布団に入ってから思いついたことを今朝起きた直後に思い出したのでメモしておきます。


おとといのブログの、さらに簡略版を思いつきました。
量子数が偶数の場合と奇数の場合(波動関数が奇関数と偶関数)を場合分けしながらいっぺんに計算する方法です


脳に過剰な負荷がかかるため、あまりお勧めはしません。笑

有限深さ井戸型ポテンシャルの波動方程式を、対称性を用いて簡易に解く

有限深さ井戸型ポテンシャルは以前解いてブログに載せたのですが
そのときは力技で解いていました。

もし、偶数番目の量子数で波動関数が奇関数、奇数番目の量子数で偶関数になるとわかっているなら、どれくらい簡略化されるのか
これについてまだブログには載せていませんでした。



幅2aの井戸型ポテンシャルの、中のポテンシャルはV=0、外のポテンシャルはV=V0とし、
井戸の中心をx=0とします。(-a≦x≦a)


x<0の場合は、偶関数なら波動関数f(-x)=f(x)とし、奇関数ならf(-x)=-f(x)とするわけです。


境界条件はx=±aでf(±a)とその微分f'(±a)が連続であることですが
x<0を考慮しなくて済むので、

4本の連立方程式が2本まで簡略化できます。
つまり、4次ではなく2次の行列方程式を解けばいいだけになるので、かなり楽です。

波数をk、積分定数をA,Bとして(iは虚数単位)
f(x)=Aexp(ikx)+Bexp(-ikx)

という波動関数の境界条件を解けばいいのですが、
左右対称か左右反対称どちらかの解だけになるので、波動関数fは三角関数
井戸の中はf(x)=Asinkxかf(x)=Bcoskxのどちらかだけになりますし、
井戸の外はf(x)=Cexp(-k1x)になります。


量子数が奇数の場合、波動関数は偶関数なので
f(a)=Bcoska=Cexp(-k1a)
f'(a)=-kBsinka=-k1Cexp(-k1a)

量子数が偶数の場合、波動関数は奇関数なので
f(a)=Asinka=Cexp(-k1a)
f'(a)=kAcoska=-k1Cexp(-k1a)

これを解けばいいだけになります。



太線から下は、無限深さ井戸型ポテンシャルへの近似で、検算です。
k1=∞になるので
n番目の量子数がそのまま量子数nとなります。
量子数nが偶数だと波動関数が奇関数になり
量子数nが奇数だと波動関数が偶関数になることが確認できました

ターレいいよねターレ

ブラシネスターレットトラック
首を何回回しても、元々断線してるから断線しないデュラハンなターレットトラック

モータじゃないんだから、ブラシがあっても電子制御していいもん！

やはり現実という世界線はフィクションの1つでしかないのか

俺の地元に、俺なんかよりもはるかに個性的な人がいるなんて思ってもみなかった。
負けたよ。完敗だ。悔しいなあ

たかが地元でこれか…世界って広いんだな