時間依存シュレーディンガー方程式（波束の散乱）をExcel VBAでルンゲクッタ8次で計算してGIFアニメに。
だいぶ前に書いた話をGIFアニメにしてみる。シュレーディンガー方程式も空間微分を差分化すると普通の常微分方程式になるので、ルンゲクッタ8次のDormand-Prince法で計算できる。
具体的には空間の拡散項は2次の差分
(φi+1 - 2*φi + φi-1) / (dx*dx)
で近似。中央に井戸型ポテンシャルを置いて、左から波束が入ってくるときの計算。
これは”計算物理”にでてたもので、もともとはシッフの例題にあるものだそうです。
量子力学の本で私が勉強したシッフとかメシアとか、朝永さんとかは最近の人は読まないんだろうなあ、、、JJサクライさんのは今でもOK?
トンネル効果の計算もしてみた。それはこちら。
https://sci.tea-nifty.com/blog/2013/01/excelvba8gif-07.html
« Sパラメータのスミスチャートを描くExcelシートを作った。アップロードしときます。 | トップページ | 布亀の救急箱のCMがなぜか大晦日にたくさん流れていますが、あの赤ちゃんは迅くん。 »
「学問・資格」カテゴリの記事
- 高周波（RF・マイクロ波・ミリ波・5G）関連ニュース2020年9月17日 IEEE microwave magazine10月号はマイクロ波顕微鏡特集、 Microwave Journalはミリ波パッケージ。Vayyarの4Dイメージングレーダの仕組み、MIPIのA-PHY、など。(2020.09.17)
- 高周波（RF・マイクロ波・ミリ波・5G）関連ニュース2020年8月17日 Microwave Journal 8月号はドローンで衛星通信をテスト、220GHzの半導体テストシステム。DesignConで導体の粗さモデルと因果律について解説。Qualcommとワシントンポストの5Gクラス、Samsung Galaxy Watch3分解(2020.08.17)
- 高周波（RF・マイクロ波・ミリ波・5G）関連ニュース2020年8月10日 Microwave Magazine 9月号はシリコンフォトニクスをマイクロ波・ミリ波通信に使う、コロナと5G、IEEE RFICでimecが低消費電力60GHzレーダを発表、Qorvoの5G向けフィルタ解説、など。(2020.08.11)
- Python＋Scipyでルンゲクッタ8次のDOP853(Dormand＆Prince)を使う（その16) 時間依存のシュレーディンガー方程式のトンネル効果を計算してGIFアニメに。(2020.07.20)
- Python＋Scipyでルンゲクッタ8次のDOP853(Dormand＆Prince)を使う（その15) 蔵本モデルを計算してGIFアニメに。 (2020.07.18)
« Sパラメータのスミスチャートを描くExcelシートを作った。アップロードしときます。 | トップページ | 布亀の救急箱のCMがなぜか大晦日にたくさん流れていますが、あの赤ちゃんは迅くん。 »
コメント