Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revisionPrevious revision
Next revision		Previous revision
documentation:standard_operators:coulomb_repulsion [2017/02/27 12:42] – Maurits W. Haverkortdocumentation:standard_operators:coulomb_repulsion [2017/05/23 16:43] (current) – Maurits W. Haverkort
Line 44: Line 44:
 The radial part of the operator ($\frac{\mathrm{Min}[r_i,r_j]^k}{\mathrm{Max}[r_i,r_j]^{k+1}}$) is more difficult and will be cast into parameters: The radial part of the operator ($\frac{\mathrm{Min}[r_i,r_j]^k}{\mathrm{Max}[r_i,r_j]^{k+1}}$) is more difficult and will be cast into parameters:
 \begin{equation} \begin{equation}
-R^{(k)}[\tau_1\tau_2\tau_3\tau_4]=e^2\int_0^{\infty}\int_0^{\infty}\frac{\mathrm{Min}[r_i,r_j]^k}{\mathrm{Max}[r_i,r_j]^{k+1}}R_1[r_i]R_2[r_j]R_3[r_i]R_4[r_j]\mathrm{d}r_i\mathrm{d}r_j.+R^{(k)}[\tau_1\tau_2\tau_3\tau_4]=e^2\int_0^{\infty}\int_0^{\infty}\frac{\mathrm{Min}[r_i,r_j]^k}{\mathrm{Max}[r_i,r_j]^{k+1}}R_1[r_i]R_2[r_j]R_3[r_i]R_4[r_j]r_i^2 r_j^2\mathrm{d}r_i\mathrm{d}r_j.
 \end{equation} \end{equation}
  
Line 63: Line 63:
  
 ### ###
-For the case where $n_1=n_2=n_3=n_4$ and $l_1=l_2=l_3=l_4$, i.e. Coulomb repulsion within one shell one defines:+{{:documentation:standard_operators:coulomb_diagram_lll.png?nolink&200 |}}For the case where $n_1=n_2=n_3=n_4$ and $l_1=l_2=l_3=l_4$, i.e. Coulomb repulsion within one shell one defines:
 \begin{equation} \begin{equation}
 F^{(k)} = R^{(k)}[\tau_1\tau_2\tau_3\tau_4]. F^{(k)} = R^{(k)}[\tau_1\tau_2\tau_3\tau_4].
Line 86: Line 86:
 ===== Two shells, shell occupation conserving ===== ===== Two shells, shell occupation conserving =====
  
-{{:documentation:standard_operators:coulomb_diagram_ll.png?nolink|}} 
  
 ### ###
-The Coulomb repulsion between two shells which does not change the number of electrons is given by a direct term ($l_1=l_3$ and $l_2=l_4$) and an indirect or exchange term ($l_1=l_4$ and $l_2=l_3$). The direct term is given by the Slater integrals:+{{:documentation:standard_operators:coulomb_diagram_ll.png?nolink&400 |}}The Coulomb repulsion between two shells which does not change the number of electrons is given by a direct term ($n_1l_1=n_3l_3$ and $n_2l_2=n_4l_4$) and an indirect or exchange term ($n_1l_1=n_4l_4$ and $n_2l_2=n_3l_3$). We here assume that $n_1l_1\neq n_2l_2$. The direct term is given by the Slater integrals:
 \begin{equation} \begin{equation}
 F^{(k)}=e^2\int_0^{\infty}\int_0^{\infty}\frac{\mathrm{Min}[r_i,r_j]^k}{\mathrm{Max}[r_i,r_j]^{k+1}}R_1[r_i]^2R_2[r_j]^2\mathrm{d}r_i\mathrm{d}r_j, F^{(k)}=e^2\int_0^{\infty}\int_0^{\infty}\frac{\mathrm{Min}[r_i,r_j]^k}{\mathrm{Max}[r_i,r_j]^{k+1}}R_1[r_i]^2R_2[r_j]^2\mathrm{d}r_i\mathrm{d}r_j,
Line 118: Line 117:
  
 ===== General case of 4 different shells ===== ===== General case of 4 different shells =====
- 
-{{:documentation:standard_operators:coulomb_diagram_llll.png?nolink|}} 
  
 +###
 +{{:documentation:standard_operators:coulomb_diagram_llll.png?nolink&200 |}}The Coulomb repulsion in the general case allows four different principle quantum numbers and angular momenta. The radial integral is given as:
 +\begin{equation}
 +R^{(k)}[n_1l_1\:n_2l_2\:n_3l_3\:n_4l_4]=e^2\int_0^{\infty}\int_0^{\infty}\frac{\mathrm{Min}[r_i,r_j]^k}{\mathrm{Max}[r_i,r_j]^{k+1}}R_{n_1l_1}[r_i]R_{n_2l_2}[r_j]R_{n_3l_3}[r_i]R_{n_4l_4}[r_j]\mathrm{d}r_i\mathrm{d}r_j,
 +\end{equation}
 +with $\mathrm{Max}[|l_1-l_3|,|l_2-l_4|] \leq k \leq \mathrm{Min}[l_1+l_3,l_2+l_4]$ and $(l_1+l_3)$, $(l2_+l_4)$ either both even or both odd.
 ### ###
  
 +###
 +In Quanty one can implement these operators as:
 +<code Quanty Example.Quanty>
 +NewOperator("U", NF, IndexUp_1, IndexDn_1, IndexUp_2, IndexDn_2, IndexUp_3, IndexDn_3, IndexUp_4, IndexDn_4, Rk)
 +</code>
 +For $l_1=3$, $l_2=0$, $l_3=2$ and $l_4=1$ one has $k=1$ and one could define:
 +<code Quanty Example.Quanty>
 +OppR1pd = NewOperator("U", NF, IndexUp_1, IndexDn_1, IndexUp_2, IndexDn_2, IndexUp_3, IndexDn_3, IndexUp_4, IndexDn_4, {1})
 +</code>
 ### ###
  
 +###
 +Note that in the general case you need to sum over all possible permutations of $n_1l_1$, $n_2l_2$, $n_3l_3$ and $n_4l_4$. Permuting $n_1l_1$ with $n_2l_2$ and at the same time $n_3l_3$ with $n_4l_4$ will not change the value and form of the operator. If $n_1l_1$ is different from $n_2l_2$ and $n_3l_3$ is different from $n_4l_4$ one can add a factor of two in front of the operator and only add one of the permutations. If one of the $n_1l_1$ is the same as $n_2l_2$ or $n_3l_3$ is the same as $n_4l_4$ a permutation will not lead to a new configuration and the factor of two disappears. If you just sum over all possible $n_il_i$ combinations things go right automatically. 
 +###
  
 ===== Table of contents ===== ===== Table of contents =====
 {{indexmenu>.#1|msort nsort}} {{indexmenu>.#1|msort nsort}}

You are here: Quanty - a quantum many body script language » Documentation » Standard Operators » Coulomb repulsion operator (U)

Print/export