An accurate equation of state for the one component plasma in the low coupling regime

An accurate equation of state of the one component plasma is obtained in the low coupling regime $0 \leq \Gamma \leq 1$. The accuracy results from a smooth combination of the well-known hypernetted chain integral equation, Monte Carlo simulations and asymptotic analytical expressions of the excess internal energy $u$. In particular, special attention has been brought to describe and take advantage of finite size effects on Monte Carlo results to get the thermodynamic limit of $u$. This combined approach reproduces very accurately the different plasma correlation regimes encountered in this range of values of $\Gamma$. This paper extends to low $\Gamma$'s an earlier Monte Carlo simulation study devoted to strongly coupled systems for $1 \leq \Gamma \leq 190$ ({J.-M. Caillol}, {J. Chem. Phys.} \textbf{111}, 6538 (1999)). Analytical fits of $u(\Gamma)$ in the range $0 \leq \Gamma \leq 1$ are provided with a precision that we claim to be not smaller than $p= 10^{-5}$. HNC equation and exact asymptotic expressions are shown to give reliable results for $u(\Gamma)$ only in narrow $\Gamma$ intervals, i.e. $0 \leq \Gamma \lesssim 0.5$ and $0 \leq \Gamma \lesssim 0.3$ respectively

Monte Carlo simulations of the screening potential of the Yukawa one-component plasma

A Monte Carlo scheme to sample the screening potential H(r) of Yukawa plasmas notably at short distances is presented. This scheme is based on an importance sampling technique. Comparisons with former results for the Coulombic one-component plasma are given. Our Monte Carlo simulations yield an accurate estimate of H(r) as well for short range and long range interparticle distances.Comment: to be published in Journal of Physics A: Mathematical and Genera

Ionic fluids: charge and density correlations near gas-liquid criticality

The correlation functions of an ionic fluid with charge and size asymmetry are studied within the framework of the random phase approximation. The results obtained for the charge-charge correlation function demonstrate that the second-moment Stillinger-Lovett (SL) rule is satisfied away from the gas-liquid critical point (CP) but not, in general, at the CP. However in the special case of a model without size assymetry the SL rules are satisfied even at the CP. The expressions for the density-density and charge-density correlation functions valid far and close to the CP are obtained explicitely

How Multivalency controls Ionic Criticality

To understand how multivalency influences the reduced critical temperatures, Tce (z), and densities, roce (z), of z : 1 ionic fluids, we study equisized hard-sphere models with z = 1-3. Following Debye, Hueckel and Bjerrum, association into ion clusters is treated with, also, ionic solvation and excluded volume. In good accord with simulations but contradicting integral-equation and field theories, Tce falls when z increases while roce rises steeply: that 80-90% of the ions are bound in clusters near T_c serves to explain these trends. For z \neq 1 interphase Galvani potentials arise and are evaluated.Comment: 4 pages, 4 figure

Discretization Dependence of Criticality in Model Fluids: a Hard-core Electrolyte

Grand canonical simulations at various levels, $\zeta=5$-20, of fine- lattice discretization are reported for the near-critical 1:1 hard-core electrolyte or RPM. With the aid of finite-size scaling analyses it is shown convincingly that, contrary to recent suggestions, the universal critical behavior is independent of $\zeta$ (\grtsim 4); thus the continuum $(\zeta\to\infty)$ RPM exhibits Ising-type (as against classical, SAW, XY, etc.) criticality. A general consideration of lattice discretization provides effective extrapolation of the {\em intrinsically} erratic $\zeta$-dependence, yielding (\Tc^ {\ast},\rhoc^{\ast})\simeq (0.0493_{3},0.075) for the $\zeta=\infty$ RPM.Comment: 4 pages including 4 figure

Scalar Casimir Effect on a D-dimensional Einstein Static Universe

We compute the renormalised energy momentum tensor of a free scalar field coupled to gravity on an (n+1)-dimensional Einstein Static Universe (ESU), RxS^n, with arbitrary low energy effective operators (up to mass dimension n+1). A generic class of regulators is used, together with the Abel-Plana formula, leading to a manifestly regulator independent result. The general structure of the divergences is analysed to show that all the gravitational couplings (not just the cosmological constant) are renormalised for an arbitrary regulator. Various commonly used methods (damping function, point-splitting, momentum cut-off and zeta function) are shown to, effectively, belong to the given class. The final results depend strongly on the parity of n. A detailed analytical and numerical analysis is performed for the behaviours of the renormalised energy density and a quantity `sigma' which determines if the strong energy condition holds for the `quantum fluid'. We briefly discuss the quantum fluid back-reaction problem, via the higher dimensional Friedmann and Raychaudhuri equations, observe that equilibrium radii exist and unveil the possibility of a `Casimir stabilisation of Einstein Static Universes'.Comment: 37 pages, 15 figures, v2: minor changes in sections 1, 2.5, 3 and 4; version published in CQ

Standard Cosmological Evolution in a Wide Range of f(R) Models

Using techniques from singular perturbation theory, we explicitly calculate the cosmological evolution in a class of modified gravity models. By considering the (m)CDTT model, which aims to explain the current acceleration of the universe with a modification of gravity, we show that Einstein evolution can be recovered for most of cosmic history in at least one f(R) model. We show that a standard epoch of matter domination can be obtained in the mCDTT model, providing a sufficiently long epoch to satisfy observations. We note that the additional inverse term will not significantly alter standard evolution until today and that the solution lies well within present constraints from Big Bang Nucleosynthesis. For the CDTT model, we analyse the ``recent radiation epoch'' behaviour (a \propto t^{1/2}) found by previous authors. We finally generalise our findings to the class of inverse power-law models. Even in this class of models, we expect a standard cosmological evolution, with a sufficient matter domination era, although the sign of the additional term is crucial.Comment: 15 pages, 6 figures (1 new figure), new version considers both CDTT and mCDTT models. References added. Accepted by Phys Rev

Universality class of criticality in the restricted primitive model electrolyte

The 1:1 equisized hard-sphere electrolyte or restricted primitive model has been simulated via grand-canonical fine-discretization Monte Carlo. Newly devised unbiased finite-size extrapolation methods using temperature-density, (T, rho), loci of inflections, Q = ^2/ maxima, canonical and C_V criticality, yield estimates of (T_c, rho_c) to +- (0.04, 3)%. Extrapolated exponents and Q-ratio are (gamma, nu, Q_c) = [1.24(3), 0.63(3); 0.624(2)] which support Ising (n = 1) behavior with (1.23_9, 0.630_3; 0.623_6), but exclude classical, XY (n = 2), SAW (n = 0), and n = 1 criticality with potentials phi(r)>Phi/r^{4.9} when r \to \infty

Equation of state in the fugacity format for the two-dimensional Coulomb gas

We derive the exact general form of the equation of state, in the fugacity format, for the two-dimensional Coulomb gas. Our results are valid in the conducting phase of the Coulomb gas, for temperatures above the Kosterlitz-Thouless transition. The derivation of the equation of state is based on the knowledge of the general form of the short-distance expansion of the correlation functions of the Coulomb gas. We explicitly compute the expansion up to order $O(\zeta^6)$ in the activity $\zeta$. Our results are in very good agreement with Monte Carlo simulations at very low density

Phase Diagram of the Two Dimensional Lattice Coulomb Gas

We use Monte Carlo simulations to map out the phase diagram of the two dimensional Coulomb gas on a square lattice, as a function of density and temperature. We find that the Kosterlitz-Thouless transition remains up to higher charge densities than has been suggested by recent theoretical estimates.Comment: 4 pages, including 6 in-line eps figure