Note230 水素原子の数値計算

先日は水素原子の量子力学による厳密解を眺めてみたがそのように解析的に計算できるのは極めて限られている。実際は数値計算に頼らざる終えない。そこで今日は単純な数値計算を行ってエネルギーを求めてみようと思う。厳密解によれば基底状態(n=1)では
https://cdn-ak.f.st-hatena.com/images/fotolife/c/cat_falcon/20190805/20190805233320.jpg
なので概ねこの値が出てくるはずだ。
https://cdn-ak.f.st-hatena.com/images/fotolife/c/cat_falcon/20190805/20190805233328.jpg
後は積分数値計算するのだがら面倒なのでそのまま(賢い計算では無い)
https://cdn-ak.f.st-hatena.com/images/fotolife/c/cat_falcon/20190805/20190805233331.jpg
としてパソコンのCPUに強引に任せてしまう。求める積分
https://cdn-ak.f.st-hatena.com/images/fotolife/c/cat_falcon/20190805/20190805233339.jpg
モンテカルロ法で力任せに計算すると
https://cdn-ak.f.st-hatena.com/images/fotolife/c/cat_falcon/20190805/20190805233348.jpg
この結果は実際の値と良く一致している。