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

On quantum effects near the liquid-vapor transition in helium

The liquid-vapor transition in He-3 and He-4 is investigated by means of path-integral molecular dynamics and the quantum virial expansion. Both methods are applied to the critical isobar and the critical isochore. While previous path-integral simulations have mainly considered the lambda transition and superfluid regime in He-4, we focus on the vicinity of the critical point and obtain good agreement with experimental results for the molar volume and the internal energy down to subcritical temperatures. We find that an effective classical potential that properly describes the two-particle radial distribution function exhibits a strong temperature dependence near the critical temperature. This contrasts with the behavior of essentially classical systems like xenon, where the effective potential is independent of temperature. It is conjectured that, owing to this difference in behavior between classical and quantum-mechanical systems, the crossover behavior observed for helium in the vicinity of the critical point differs qualitatively from that of other simple liquids

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

Comment on ``Scaling Laws for a System with Long-Range Interactions within Tsallis Statistics''

In their recent Letter [Phys. Rev. Lett. 83, 4233 (1999)], Salazar and Toral (ST) study numerically a finite Ising chain with non-integrable interactions decaying like 1/r^(d+sigma) where -d <= sigma <= 0 (like ST, we discuss general dimensionality d). In particular, they explore a presumed connection between non-integrable interactions and Tsallis's non-extensive statistics. We point out that (i) non-integrable interactions provide no more motivation for Tsallis statistics than do integrable interactions, i.e., Gibbs statistics remain meaningful for the non-integrable case, and in fact provide a {\em complete and exact treatment}; and (ii) there are undesirable features of the method ST use to regulate the non-integrable interactions.Comment: Accepted for publication in Phys. Rev. Let

Stabilization of colloidal suspensions by means of highly-charged nanoparticles

We employ a novel Monte Carlo simulation scheme to elucidate the stabilization of neutral colloidal microspheres by means of highly-charged nanoparticles [V. Tohver et al., Proc. Natl. Acad. Sci. U.S.A. 98, 8950 (2001)]. In accordance with the experimental observations, we find that small nanoparticle concentrations induce an effective repulsion that prevents gelation caused by the intrinsic van der Waals attraction between colloids. Higher nanoparticle concentrations induce an attractive potential which is, however, qualitatively different from the regular depletion attraction. We also show how colloid-nanoparticle size asymmetry and nanoparticle charge can be used to manipulate the effective interactions.Comment: Accepted for publication in Physical Review Letters. See also S. Karanikas and A.A. Louis, cond-mat/0411279. Updated to synchronize with published versio

Monte Carlo cluster algorithm for fluid phase transitions in highly size-asymmetrical binary mixtures

Highly size-asymmetrical fluid mixtures arise in a variety of physical contexts, notably in suspensions of colloidal particles to which much smaller particles have been added in the form of polymers or nanoparticles. Conventional schemes for simulating models of such systems are hamstrung by the difficulty of relaxing the large species in the presence of the small one. Here we describe how the rejection-free geometrical cluster algorithm (GCA) of Liu and Luijten [Phys. Rev. Lett 92, 035504 (2004)] can be embedded within a restricted Gibbs ensemble to facilitate efficient and accurate studies of fluid phase behavior of highly size-asymmetrical mixtures. After providing a detailed description of the algorithm, we summarize the bespoke analysis techniques of Ashton et al. [J. Chem. Phys. 132, 074111 (2010)] that permit accurate estimates of coexisting densities and critical-point parameters. We apply our methods to study the liquid--vapor phase diagram of a particular mixture of Lennard-Jones particles having a 10:1 size ratio. As the reservoir volume fraction of small particles is increased in the range 0--5%, the critical temperature decreases by approximately 50%, while the critical density drops by some 30%. These trends imply that in our system, adding small particles decreases the net attraction between large particles, a situation that contrasts with hard-sphere mixtures where an attractive depletion force occurs.Comment: 11 pages, 10 figure

Monte Carlo Renormalization Group Analysis of Lattice ϕ4\phi^4 Model in D=3,4D=3,4

We present a simple, sophisticated method to capture renormalization group flow in Monte Carlo simulation, which provides important information of critical phenomena. We applied the method to $D=3,4$ lattice $\phi^4$ model and obtained renormalization flow diagram which well reproduces theoretically predicted behavior of continuum $\phi^4$ model. We also show that the method can be easily applied to much more complicated models, such as frustrated spin models.Comment: 13 pages, revtex, 7 figures. v1:Submitted to PRE. v2:considerably reduced redundancy of presentation. v3:final version to appear in Phys.Rev.

Quantum spin chains with site dissipation

We use Monte Carlo simulations to study chains of Ising- and XY-spins with dissipation coupling to the site variables. The phase diagram and critical exponents of the dissipative Ising chain in a transverse magnetic field have been computed previously, and here we consider a universal ratio of susceptibilities. We furthermore present the phase diagram and exponents of the dissipative XY-chain, which exhibits a second order phase transition. All our results compare well with the predictions from a dissipative $\phi^4$ field theory