Screening in ionic systems: simulations for the Lebowitz length

Simulations of the Lebowitz length, xiL (T, rho), are reported for the restricted primitive model hard-core (diameter a) 1:1 electrolyte for densities rho approximately < 4rho(c) and T(c) approximately < T approximately < 40T(c). Finite-size effects are elucidated for the charge fluctuations in various subdomains that serve to evaluate xiL. On extrapolation to the bulk limit for T approximately > 10T(c) the exact low-density expansions are seen to fail badly when rho > 1/10 rho(c) (with rho(c)a3 approximately = 0.08). At higher densities xiL rises above the Debye length, xiD proportional to square root(T/rho), by 10%-30% (up to rho approximately =1.3rho(c)); the variation is portrayed fairly well by the generalized Debye-Hückel theory. On approaching criticality at fixed rho or fixed T, xiL (T, rho) remains finite with xiL(c) approximately = 0.30a approximately = 1.3xiD(c) but displays a weak entropylike singularity.