- #1
Asban
- 7
- 0
Hello,
I'm new in the realm of numerical analysis.
I need to find the local probability density of energy P(E) of a single electron on site i from a self-consistent equation for the energy (mean field equation of energy).
E_i=\epsilon_i+\sum_j(\frac{1}{1+e^{E_j/T}}-\frac{1}{2})\frac{e^2}{r_{ij}}
\epsilon_i - energy at site i with random values between -1/2 to 1/2
E_j - average energy of site j
E_i - average energy of site j
r_{nn} - average distance of nearest neighbors
\frac{e^2}{r_{nn}} - coulomb potential energy between nearest neighbors
There are 10,000 sites
The sites were uniformly distributes on a 2D sample
The values of temperature and r_{nn} are \frac{e^2}{r_{nn}T}=20
Now the only explanation of how probability density ( P(E) ) was found is:
1."...starting with a random set of energies and evolving them iteratively within the mean-field model"
2."solving the equations for many instances and averaging over them"
I understand the physics of this equation and how they got to it and I'm searching for a good refrence that explains this general numerical method in details or someone that can explain that method.
Most regards
Ofek
I'm new in the realm of numerical analysis.
I need to find the local probability density of energy P(E) of a single electron on site i from a self-consistent equation for the energy (mean field equation of energy).
E_i=\epsilon_i+\sum_j(\frac{1}{1+e^{E_j/T}}-\frac{1}{2})\frac{e^2}{r_{ij}}
\epsilon_i - energy at site i with random values between -1/2 to 1/2
E_j - average energy of site j
E_i - average energy of site j
r_{nn} - average distance of nearest neighbors
\frac{e^2}{r_{nn}} - coulomb potential energy between nearest neighbors
There are 10,000 sites
The sites were uniformly distributes on a 2D sample
The values of temperature and r_{nn} are \frac{e^2}{r_{nn}T}=20
Now the only explanation of how probability density ( P(E) ) was found is:
1."...starting with a random set of energies and evolving them iteratively within the mean-field model"
2."solving the equations for many instances and averaging over them"
I understand the physics of this equation and how they got to it and I'm searching for a good refrence that explains this general numerical method in details or someone that can explain that method.
Most regards
Ofek