565 research outputs found
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
- …