XMCDを(再度)作る。ピークの位置を求める。step関数をXASから引いて積分する。境目を決めてXMCDを積分する。edgeを以下の関数を使って、求める。b要は、各エネルギー領域を指定してあげて、その中のデータ点を用いて最小二乗法で直線を定めるというもの。eVくらいで設定した。プラスhelicityのXASにのみ適用してみる。[eV]")plt.tight_layout()plt.show()うまく引けている感じがする。規格化後のXASとXMCD。norm_mピークの位置を、次の関数で求める。eVの範囲にあるので、そのように設定。740.54...d状態以外の寄与を差し引くためのstep関数を次のように用意。]になるように調整した。(この関数は、(遷移金属原子における)d状態以外の寄与のためのモデルであり、それらは一般に軌道モーメントは小さい。)求めたピークの位置を使って、step関数を作る。)plt.show()それらしく階段関数が置けている。curve_fitの使い方は以下のサイトを参照。おっぱいそん!階段関数を差し引いたXAS、そしてXMCDを次のように積分する。ここではSimpson法を用いる。18.914...ちょっと後ろの方でフニャッとしているが、こんなもんだろう。境目は、XASの一回微分がゼロになるところを持って来た。-4.5744...差分による近似的に微分した配列が得られるnp.diffの使い方は、以下のサイトを参照。ruleに必要なXAS及びXMCDの積分値を求めることが出来た。ちなみに、適当に作ったスペクトルは以下のスクリプトを使用。backローレンツ関数の大きさ等は以下を参照。定数ポテンシャルを与える球対称電子密度について。化学ポテンシャルの温度依存性。Sommerfeld展開。OPW(直交化平面波法)月別アーカイブはてなブログをはじめよう!koideforestさんは、はてなブログを使っています。ブログを報告するスターをつけました引用をストックしましたストック一覧を見る閉じる引用するにはまずログインしてくださいログイン閉じる引用をストックできませんでした。再度お試しください閉じる限定公開記事のため引用できません。