Potentials
Content
Potentials#
Sarkas supports a variety of potentials both built-in and user defined. Currently the potential functions that are implemented include
You can read more about each of these potentials in the corresponding sections below. All the equations will be given in cgs units, however, for easy conversion, we define the charge
which when substituted in gives the equivalent mks formula.
Electron parameters and thermodynamic formulas are given in here.
Coulomb Potential#
Two charged particle with charge numbers \(Z_a\) and \(Z_b\) interact with each other via the Coulomb potential given by
where \(r\) is the distance between ions, \(e\) is the elementary charge.
Yukawa Potential#
The Yukawa potential, or screened Coulomb potential, is widely used in the plasma community to describe the interactions of positively charged ions in a uniform background of electrons. The form of the Yukawa potential for two ions of charge number \(Z_a\) and \(Z_b\) is given by
where \(\lambda_{\textrm{TF}}\) is the Thomas-Fermi wavelength and \(\kappa\) is the screening parameter. In Sarkas \(\kappa\) can be given as an input or it can be calculated from the Thomas-Fermi Wavelength formula.
Notice that when \(\kappa = 0\) we recover the Coulomb Potential.
Exact Gradient-corrected Screened Yukawa Potential#
The Yukawa potential is derived on the assumption that the electron gas behaves as an ideal Fermi gas. Improvements in this theory can be achieved by considering density gradients and exchange-correlation effects. Stanton and Murillo [Stanton and Murillo, 2015], using a DFT formalism, derived an exact-gradient corrected ion pair potential across a wide range of densities and temperatures.
The exact-gradient screened (EGS) potential introduces new parameters that can be easily calculated from initial inputs. Density gradient corrections to the free energy functional lead to the first parameter, \(\nu\),
where \(\lambda\) is a correction factor; \(\lambda = 1/9\) for the true gradient corrected Thomas-Fermi model and \(\lambda = 1\) for the traditional von Weissaecker model, \(\mathcal I_{-1/2}[\eta_0]\) is the Fermi Integral of order \(-1/2\), and \(\Lambda_e\) is the de Broglie wavelength of the electrons.
In the case \(\nu < 1\) the EGS potential takes the form
with
where the parameter \(b\) arises from exchange-correlation contributions, see below. On the other hand \(\nu > 1\), the pair potential has the form
with
Neglect of exchange-correlational effects leads to \(b = 1\) otherwise
where \(k_{\textrm{F}}\) is the Fermi wavenumber and \(\Theta = (\beta E_{\textrm{F}})^{-1}\) is the electron Degeneracy Parameter calculated from the Fermi Energy.
Quantum Statistical Potentials#
An extensive review on Quantum Statistical Potentials is given in [Jones and Murillo, 2007]. The following module uses that as the main reference.
Quantum Statistical Potentials are defined by three terms
where
is due to the Pauli exclusion principle and it accounts for spin-averaged effects,
is the usual Coulomb interaction between two charged particles with charge numbers \(Z_a,Z_b\), and \(U_{\textrm{diff}}(r)\) is a diffraction term. There are two possibilities for the diffraction term. The most common is the Deutsch potential
The second most common form is the Kelbg potential
In the above equations the screening length \(\Lambda_{ab}\) is the thermal de Broglie wavelength between the two charges defined as
Note that the de Broglie wavelength is defined differently in [Hansen and McDonald, 1981] hence the factor of \(2\pi\) in the exponential.
The long range part of the potential is computed using the PPPM algorithm where only the \(U_{\textrm{coul}}(r)\) term is split into a short range and long range part.
The choice of this potential is due to its widespread use in the High Energy Density Physics community.
Moliere Potential#
Moliere-type potentials have the form
with the contraint
more info can be found in [Wilson et al., 1977]
Lennard Jones#
Sarkas support the general form of the multispecies Lennard Jones potential
where
In the case of multispecies liquids we use the Lorentz-Berthelot mixing rules