時間依存シュレーディンガー方程式（波束の散乱）をExcel VBAでルンゲクッタ8次で計算してGIFアニメに。
だいぶ前に書いた話をGIFアニメにしてみる。シュレーディンガー方程式も空間微分を差分化すると普通の常微分方程式になるので、ルンゲクッタ8次のDormand-Prince法で計算できる。
具体的には空間の拡散項は2次の差分
(φi+1 - 2*φi + φi-1) / (dx*dx)
で近似。中央に井戸型ポテンシャルを置いて、左から波束が入ってくるときの計算。
これは”計算物理”にでてたもので、もともとはシッフの例題にあるものだそうです。
量子力学の本で私が勉強したシッフとかメシアとか、朝永さんとかは最近の人は読まないんだろうなあ、、、JJサクライさんのは今でもOK?
トンネル効果の計算もしてみた。それはこちら。
http://sci.tea-nifty.com/blog/2013/01/excelvba8gif-07.html
« Sパラメータのスミスチャートを描くExcelシートを作った。アップロードしときます。 | トップページ | 布亀の救急箱のCMがなぜか大晦日にたくさん流れていますが、あの赤ちゃんは迅くん。 »
「学問・資格」カテゴリの記事
- #仮面ライダービルド 第40話の黒板の数式は、、、一瞬しか映らずわからなかった、、、ハミルトニアンと、L,Rが見えたので量子もつれ関連の何か？→全然違いました。密度行列繰り込み群か！(2018.06.17)
- #仮面ライダービルド 第40話の話数を表す数式はペル方程式x^2-77y^2=1の最小解がx=351,y=40ということ。(2018.06.17)
- #AKB総選挙 2018年の結果がジップの法則（というかべき乗則）に乗るかを検証。一位は松井珠理奈さんでした。3位、7位、15位で関数形が変わる。(2018.06.16)
- Powerpointで数式入力をLaTeX形式で行う。WordはデフォルトでできるけどPowerpointはちょっと初期設定が必要。(2018.06.13)
- #仮面ライダービルド 第39話の黒板の数式はワインバーグ・サラム理論（＆ヒッグス機構）？(2018.06.10)
トラックバック
この記事のトラックバックURL:
http://app.cocolog-nifty.com/t/trackback/512682/56425297
この記事へのトラックバック一覧です： 時間依存シュレーディンガー方程式（波束の散乱）をExcel VBAでルンゲクッタ8次で計算してGIFアニメに。:
« Sパラメータのスミスチャートを描くExcelシートを作った。アップロードしときます。 | トップページ | 布亀の救急箱のCMがなぜか大晦日にたくさん流れていますが、あの赤ちゃんは迅くん。 »
コメント