パッと探したところ と出来る理由が見当たらなかったので、真面目に計算してみた。 （もちろん、ハミルトニアンの並進対称性から明らかなのはわかっているが、真面目な方法でやって同じ答えになるはずである）Diracのbraket表示を使って、演算子の方程式から…

※訂正: マイナスの境界条件で内向き波になるように訂正(2019/03/08)Hartree原子単位系以下、これを証明する。フーリエ変換を使って、微分方程式を解く。 平面波の角運動量展開 パリティ対称性 より、被積分関数はに対して偶関数。よって、 球ベッセル関数は…

実験結果か何かで、散乱振幅の角度分布が得られているとする。 ここからサイト行列を復元する方法を考える。 サイト行列が求まれば、そこからArgand diagramを作って、（準）共鳴準位があるとかないとか、追加で情報が得られる（もし角度分布のエネルギー依…

仕事で図を作る必要があったので、練習がてらまとめた。位相シフトなどの概念にについては、wikipedia等を参照して頂きたい。 散乱振幅 - Wikipedia 簡単に言うと、「ポテンシャルが無い時の波は、ポテンシャルによってどれくらい変調されたか？」を表す量で…

以前、第二量子化についての記事を書いたが、誤って消してしまったので、再度投稿。任意の波動関数を、基底関数で展開し、その展開係数をとする。 これに対応する場の演算子は、消滅演算子を用いて、次のように書ける。したがって、第二量子化は展開係数を量…

以前、pythonにおける水素様原子のためのLaguerre陪関数について確認した。 koideforest.hatenadiary.comせっかくなので、pythonで規格化まで含めた水素様原子波動関数のスクリプト（原子単位系）を作った。 from scipy.special import assoc_laguerre from …

ほぼ自分用のメモ from scipy.special import spherical_jn, spherical_yn, sph_harm from sympy import * from sympy.physics.wigner import clebsch_gordan, gaunt # spherical Bessel function def jl( l, rho ): return spherical_jn( l, rho ) # spheri…

四つ足で御馴染みのメタンの分子軌道が、固体におけるLiebフラットバンドが現れる構造と同じということで、テンションが上がった。以下の簡単な飛び移りだけを考えたハミルトニアン行列の固有値と固有ベクトルを求める。 対角成分（各原子上でのエネルギー）…

N粒子系における二体ポテンシャルを、粒子ペアに番号を付けて、その間の相互作用の和（）として定義する。二体ポテンシャルとしたが、結局は相互作用に何かしら番号が付けられて、全体がその和として与えらえるのであれば、いろんな応用が考えられる（と思う…

を原点から測った距離とし、原点から見て球対称な関数をとする。 原点から見て位置にあるサイトがあるとする。 サイトを原点に取り直した任意の位置ベクトルをと定義する。 やりたいことは、をサイトから見ると球対称ではないので、サイトから見たときの球対…

例えば球面調和関数これはある意味のエンタングルメントを断ったと言えるか？ 添字が一個なので、ベクトルに拡張すれば、よって、ある種の行列化によって変数分離が可能である。ポイントは、もともと行列ではなくただの数だったので、「左に横ベクトル・右に…

Gaunt積分の自分的メモ。 Clebsch–Gordan coefficients - Wikipedia 3-j symbol - WikipediaGaunt積分（多分一般的な呼び名ではない。が、自分の分野ではよく出てくる。） Gaunt係数（Wikipediaとかに載っているのはこっち。） Clebsh-Gordan係数と3j記号の…

sympy moduleを一部使ってゴニョゴニョ計算させて、それを元にnumpy moduleで処理（例えば逆行列計算）させようとした時に虚数単位が引っかかって怒られることがあった。 sympyの代数表記を数値表記に戻す方法として"N()"や"~.evalf()"があるが、これではsym…

反変ベクトルに対する変換行列と共変ベクトルに対する変換行列は、計量テンソルで以下のように結び付けられる。しかし、やっぱりわかったようなわからないような、という感じに苛まれる。 ここでと定義すると、以下のような関係であることがわかる。つまり、…

計量テンソルを以下のように定義。それで以下の4元ベクトルの内積が、ローレンツ変換で不変であることが特殊相対論の特徴であった。で、これはアインシュタインの規約を用いているから実際には和である。 しかし、頭ではわかっているのだが、心がどうも胡散…

半径Rの構造クラスターをASEで作る。 例として単体アルミニウムを使う。 import math from ase.spacegroup import crystal from ase.visualize import view # parameters LC = 4.05 R = 7. NSL = math.ceil( 2. * R / LC ) # make bulk structure aluminium …

Python module のASE を使って、吸着構造を最適化。 ここでは簡単な有効媒質理論を使うが、calculatorを選ぶことで、AbinitとかVASPとかも使えるらしい。 Introduction: Nitrogen on copper — ASE documentation 元サイトのチュートリアルに少し追加して、吸…

#このページは誤っています。Path operator を次のように用意しておく。誤り（一般的なPath operatorの定義ではない。） 誤りはサイト間の移動、はサイト上での繰り込まれた相互作用（site T matrix）を表すとする。 後のためにを一応、site potential を使…

参考文献："Nonequilibrium Green's Functions Approach to Inhomogeneous Systems", Karsten Blazer and Michael Bonitz (Springer 2013). 該当箇所：2.1.1 Keldysh Contourこの手のもので自分が最初に詰まってしまうのは、 ハミルトニアンが時間に依存して…

零点振動している調和振動子と同じエネルギーを持つ平面波とで波動関数を比較する。 最後ら辺、サボりました。ここで、として自由平面波のエネルギーを代入して、波数と有効距離を結びつけると、 綺麗に有効距離の逆数になる。これで求めたを持つ平面波（）…

Wikipediaには日頃からよくお世話になっているが、自分で確認するのは大事だなと感じた事件。水素原子の波動関数の可視化とかカッコイイなぁと思い、ネットサーフィン中に以下のサイト様を訪問。 scipyの特殊関数 - 篠突く雨の日記自分も早速コピーしてやっ…

「スレーター型軌道の数値積分は大変」と呪文のように聞かされて続けて心にハードルが出来てしまったので、ここらで打破することを試みた。 前回、多次元のモンテカルロ積分をまとめた。 koideforest.hatenadiary.com 今回はスクリプトを組んで、どれくらい…

モンティ・ホール問題 モンティ・ホール問題 - Wikipedia モンティとはテレビの司会者的な存在らしい。 ルールは以下の通り。(1) 3つのドア (A, B, C) に（景品、ヤギ、ヤギ）がランダムに入っている。 (2) プレーヤーはドアを1つ選ぶ。 (3) モンティは残り…

そもそも分子分母が多項式の時点でやる気がなくなる。すごくやりたくない。苦手意識が強い。 でも多分それは処方箋が頭に入っていないだけ。 今回は対処法を簡単に考察する。例えば、球ベッセル関数を球ノイマン関数で割ったものにおける原点近傍の振舞が知…

q-number（演算子）をハットで表そうとしたらめっちゃズレるので、c-number（演算子じゃない普通の数）は小文字、q-numberは大文字で表すことにする。 とし、と書けるとするならば、はどんな形になるか？ということを考える。 多分これを摂動展開と呼んでい…

クーロン散乱の記事を読んでいて、よく考えると???となった位相の問題。 要するに、 (?)となるか？というもの。 いや成らんやろ。がどんなに大きかろうが、結局が重要なのであって、によってもたらされる位相のズレを解消することは出来ない。 一つ有り得る…

今日、朝起きてからずっと多電子原子波動関数の作成プログラムやっているけど、なかなかうまく行かない。 コードは引き続き下記のサイトを参照。 http://www-cms.phys.s.u-tokyo.ac.jp/~naoki/CIPINTRO/CIP/atom.htmlとりあえずs軌道は大丈夫。 でもp軌道に…

束縛状態の波動関数の計算に抵抗がありまくりんぐなので、作ってみました。 普段は連続状態を計算していて、境界条件って言っても解析解につなぐだけなので、「正しいエネルギー（固有値）を見つける」というプロセスに心のハードルが無限の井戸型ポテンシャ…

ただの愚痴であるが、位相空間が苦手である。 大抵、ミクロカノニカル分布から入ると思うが、その説明に 位相空間の導入 謎のエネルギー範囲の設定 エルゴード仮説 等重率の原理 が使われる。 その中で一番状態数の多い状態が平衡に対応するとし、変分をかま…

今更ながらIPythonデータサイエンスクックブックを購入した。もちろん私費である。せっかくなので一次元のシュレーディンガー方程式でも解こうかと思った。 一階微分方程式に対しては、かくあきさんのサイトで紹介されているローレンツアトラクターのソース…