Optimized energy calculation in lattice systems with long-range interactions

We discuss an efficient approach to the calculation of the internal energy in numerical simulations of spin systems with long-range interactions. Although, since the introduction of the Luijten-Bl\"ote algorithm, Monte Carlo simulations of these systems no longer pose a fundamental problem, the energy calculation is still an O(N^2) problem for systems of size N. We show how this can be reduced to an O(N logN) problem, with a break-even point that is already reached for very small systems. This allows the study of a variety of, until now hardly accessible, physical aspects of these systems. In particular, we combine the optimized energy calculation with histogram interpolation methods to investigate the specific heat of the Ising model and the first-order regime of the three-state Potts model with long-range interactions.Comment: 10 pages, including 8 EPS figures. To appear in Phys. Rev. E. Also available as PDF file at http://www.cond-mat.physik.uni-mainz.de/~luijten/erikpubs.htm

Critical properties of the three-dimensional equivalent-neighbor model and crossover scaling in finite systems

Accurate numerical results are presented for the three-dimensional equivalent-neighbor model on a cubic lattice, for twelve different interaction ranges (coordination number between 18 and 250). These results allow the determination of the range dependences of the critical temperature and various critical amplitudes, which are compared to renormalization-group predictions. In addition, the analysis yields an estimate for the interaction range at which the leading corrections to scaling vanish for the spin-1/2 model and confirms earlier conclusions that the leading Wegner correction must be negative for the three-dimensional (nearest-neighbor) Ising model. By complementing these results with Monte Carlo data for systems with coordination numbers as large as 52514, the full finite-size crossover curves between classical and Ising-like behavior are obtained as a function of a generalized Ginzburg parameter. Also the crossover function for the effective magnetic exponent is determined.Comment: Corrected shift of critical temperature and some typos. To appear in Phys. Rev. E. 18 pages RevTeX, including 10 EPS figures. Also available as PDF file at http://www.cond-mat.physik.uni-mainz.de/~luijten/erikpubs.htm

Universal finite-size scaling analysis of Ising models with long-range interactions at the upper critical dimensionality: Isotropic case

We investigate a two-dimensional Ising model with long-range interactions that emerge from a generalization of the magnetic dipolar interaction in spin systems with in-plane spin orientation. This interaction is, in general, anisotropic whereby in the present work we focus on the isotropic case for which the model is found to be at its upper critical dimensionality. To investigate the critical behavior the temperature and field dependence of several quantities are studied by means of Monte Carlo simulations. On the basis of the Privman-Fisher hypothesis and results of the renormalization group the numerical data are analyzed in the framework of a finite-size scaling analysis and compared to finite-size scaling functions derived from a Ginzburg-Landau-Wilson model in zero mode (mean-field) approximation. The obtained excellent agreement suggests that at least in the present case the concept of universal finite-size scaling functions can be extended to the upper critical dimensionality.Comment: revtex4, 10 pages, 5 figures, 1 tabl

Test of renormalization predictions for universal finite-size scaling functions

We calculate universal finite-size scaling functions for systems with an n-component order parameter and algebraically decaying interactions. Just as previously has been found for short-range interactions, this leads to a singular epsilon-expansion, where epsilon is the distance to the upper critical dimension. Subsequently, we check the results by numerical simulations of spin models in the same universality class. Our systems offer the essential advantage that epsilon can be varied continuously, allowing an accurate examination of the region where epsilon is small. The numerical calculations turn out to be in striking disagreement with the predicted singularity.Comment: 6 pages, including 3 EPS figures. To appear in Phys. Rev. E. Also available as PDF file at http://www.cond-mat.physik.uni-mainz.de/~luijten/erikpubs.htm

Shape of crossover between mean-field and asymptotic critical behavior in a three-dimensional Ising lattice

Recent numerical studies of the susceptibility of the three-dimensional Ising model with various interaction ranges have been analyzed with a crossover model based on renormalization-group matching theory. It is shown that the model yields an accurate description of the crossover function for the susceptibility.Comment: 4 pages RevTeX + 3 PostScript figures. Uses epsf.sty and rotate.sty. Final version; accepted for publication in Physics Letters

The heat capacity of the restricted primitive model electrolyte

The constant-volume heat capacity, C_V(T, rho), of the restricted primitive model (RPM) electrolyte is considered in the vicinity of its critical point. It is demonstrated that, despite claims, recent simulations for finite systems do not convincingly indicate the absence of a divergence in C_V(T, rho)--which would point to non-Ising-type criticality. The strong qualitative difference between C_V for the RPM and for a Lennard-Jones fluid is shown to result from the low critical density of the former. If one considers the theoretically preferable configurational heat-capacity density, C_V/V, the finite-size results for the two systems display qualitatively similar behavior on near-critical isotherms.Comment: 5 Pages, including 5 EPS figures. Also available as PDF file at http://pallas.umd.edu/~luijten/erikpubs.htm

Medium-range interactions and crossover to classical critical behavior

We study the crossover from Ising-like to classical critical behavior as a function of the range R of interactions. The power-law dependence on R of several critical amplitudes is calculated from renormalization theory. The results confirm the predictions of Mon and Binder, which were obtained from phenomenological scaling arguments. In addition, we calculate the range dependence of several corrections to scaling. We have tested the results in Monte Carlo simulations of two-dimensional systems with an extended range of interaction. An efficient Monte Carlo algorithm enabled us to carry out simulations for sufficiently large values of R, so that the theoretical predictions could actually be observed.Comment: 16 pages RevTeX, 8 PostScript figures. Uses epsf.sty. Also available as PostScript and PDF file at http://www.tn.tudelft.nl/tn/erikpubs.htm

Universality and the five-dimensional Ising model

We solve the long-standing discrepancy between Monte Carlo results and the renormalization prediction for the Binder cumulant of the five-dimensional Ising model. Our conclusions are based on accurate Monte Carlo data for systems with linear sizes up to L=22. A detailed analysis of the corrections to scaling allows the extrapolation of these results to L=\infinity. Our determination of the critical point, K_c=0.1139150 (4), is more than an order of magnitude more accurate than previous estimates.Comment: 6 pages LaTeX, 1 PostScript figure. Uses cite.sty (included) and epsf.sty. Also available as PostScript and PDF file at http://www.tn.tudelft.nl/tn/erikpubs.htm

Do crossover functions depend on the shape of the interaction profile?

We examine the crossover from classical to non-classical critical behaviour in two-dimensional systems with a one-component order parameter. Since the degree of universality of the corresponding crossover functions is still subject to debate, we try to induce non-universal effects by adding interactions with a second length scale. Although the crossover functions clearly depend on the range of the interactions, they turn out to be remarkably robust against further variation of the interaction profile. In particular, we find that the earlier observed non-monotonic crossover of the effective susceptibility exponent occurs for several qualitatively different shapes of this profile.Comment: 7 pages + 4 PostScript figures. Accepted for publication in Europhysics Letters. Also available as PDF file at http://www.cond-mat.physik.uni-mainz.de/~luijten/erikpubs.htm