Rejection-free Geometric Cluster Algorithm for Complex Fluids

We present a novel, generally applicable Monte Carlo algorithm for the simulation of fluid systems. Geometric transformations are used to identify clusters of particles in such a manner that every cluster move is accepted, irrespective of the nature of the pair interactions. The rejection-free and non-local nature of the algorithm make it particularly suitable for the efficient simulation of complex fluids with components of widely varying size, such as colloidal mixtures. Compared to conventional simulation algorithms, typical efficiency improvements amount to several orders of magnitude

Identification of the critical temperature from non-equilibrium time-dependent quantities

We present a new procedure able to identify and measure the critical temperature. This method is based on the divergence of the relaxation time approaching the critical point in quenches from infinite temperature. We introduce a dimensionless quantity that turns out to be time-independent at the critical temperature. The procedure does not need equilibration and allows for a relatively fast identification of the critical temperature. The method is first tested in the ferromagnetic Ising model and then applied to the one-dimensional Ising spin glass with power-law interactions. Here we always find a finite critical temperature also in presence of a uniform external field, in agreement with the mean-field picture for the low temperature phase of spin glasses.Comment: 6 pages, 10 figure

Critical aging of a ferromagnetic system from a completely ordered state

We adapt the non-linear $\sigma$ model to study the nonequilibrium critical dynamics of O(n) symmetric ferromagnetic system. Using the renormalization group analysis in $d=2+\epsilon$ dimensions we investigate the pure relaxation of the system starting from a completely ordered state. We find that the average magnetization obeys the long-time scaling behavior almost immediately after the system starts to evolve while the correlation and response functions demonstrate scaling behavior which is typical for aging phenomena. The corresponding fluctuation-dissipation ratio is computed to first order in $\epsilon$ and the relation between transverse and longitudinal fluctuations is discussed.Comment: 5 pages, revtex

Depinning Transition of a Two Dimensional Vortex Lattice in a Commensurate Periodic Potential

We use Monte Carlo simulations of the 2D one component Coulomb gas on a triangular lattice, to study the depinning transition of a 2D vortex lattice in a commensurate periodic potential. A detailed finite size scaling analysis indicates this transition to be first order. No significant changes in behavior were found as vortex density was varied over a wide range.Comment: 5 pages, 8 figures. Revised discussion of correlation length exponent using a more accurate finite size scaling analysis. New figs. 5 and 6. Old figs. 6 and 7 now figs. 7 and

Corrections to Scaling for the Two-dimensional Dynamic XY Model

With large-scale Monte Carlo simulations, we confirm that for the two-dimensional XY model, there is a logarithmic correction to scaling in the dynamic relaxation starting from a completely disordered state, while only an inverse power law correction in the case of starting from an ordered state. The dynamic exponent $z$ is $z=2.04(1)$.Comment: to appear as a Rapid commu. in Phys. Rev.

Measuring the equation of state of a hard-disc fluid

We use video microscopy to study a two-dimensional (2D) model fluid of charged colloidal particles suspended in water and compute the pressure from the measured particle configurations. Direct experimental control over the particle density by means of optical tweezers allows the precise measurement of pressure as a function of density. We compare our data with theoretical predictions for the equation of state, the pair-correlation function and the compressibility of a hard-disc fluid and find good agreement, both for the fluid and the solid phase. In particular the location of the transition point agrees well with results from Monte Carlo simulations.Comment: 7 pages, to appear in EPL, slightly corrected versio

Dynamical Critical Phenomena in three-dimensional Heisenberg Spin Glasses

Spin-glass (SG) and chiral-glass (CG) orderings in three dimensional (3D) Heisenberg spin glass with and without magnetic anisotropy are studied by using large-scale off-equilibrium Monte Carlo simulations. A characteristic time of relaxation, which diverges at a transition temperature in the thermodynamic limit, is obtained as a function of the temperature and the system size. Based on the finite-size scaling analysis for the relaxation time, it is found that in the isotropic Heisenberg spin glass, the CG phase transition occurs at a finite temperature, while the SG transition occurs at a lower temperature, which is compatible with zero. Our results of the anisotropic case support the chirality scenario for the phase transitions in the 3D Heisenberg spin glasses.Comment: 9 pages, 19 figure

Monte Carlo Simulations of Short-time Critical Dynamics with a Conserved Quantity

With Monte Carlo simulations, we investigate short-time critical dynamics of the three-dimensional anti-ferromagnetic Ising model with a globally conserved magnetization $m_s$ (not the order parameter). From the power law behavior of the staggered magnetization (the order parameter), its second moment and the auto-correlation, we determine all static and dynamic critical exponents as well as the critical temperature. The universality class of $m_s=0$ is the same as that without a conserved quantity, but the universality class of non-zero $m_s$ is different.Comment: to appear in Phys. Rev.

Dynamic structure factor of the Ising model with purely relaxational dynamics

We compute the dynamic structure factor for the Ising model with a purely relaxational dynamics (model A). We perform a perturbative calculation in the $\epsilon$ expansion, at two loops in the high-temperature phase and at one loop in the temperature magnetic-field plane, and a Monte Carlo simulation in the high-temperature phase. We find that the dynamic structure factor is very well approximated by its mean-field Gaussian form up to moderately large values of the frequency $\omega$ and momentum $k$. In the region we can investigate, $k\xi \lesssim 5$, $\omega \tau \lesssim 10$, where $\xi$ is the correlation length and $\tau$ the zero-momentum autocorrelation time, deviations are at most of a few percent.Comment: 21 pages, 3 figure

Hexatic Order and Surface Ripples in Spherical Geometries

In flat geometries, two dimensional hexatic order has only a minor effect on capillary waves on a liquid substrate and on undulation modes in lipid bilayers. However, extended bond orientational order alters the long wavelength spectrum of these ripples in spherical geometries. We calculate this frequency shift and suggest that it might be detectable in lipid bilayer vesicles, at the surface of liquid metals and in multielectron bubbles in liquid helium at low temperatures. Hexatic order also leads to a shift in the threshold for the fission instability induced in the later two systems by an excess of electric charge.Comment: 5 pages, 1 figure; revised version; to appear in Phys. Rev. Let