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

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

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

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

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.

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.

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

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

Characterization of relativistic electron bunch duration and travelling wave structure phase velocity based on momentum spectra measurements on the ARES linac at DESY

The ARES linac at DESY aims to generate and characterize ultrashort electron bunches (fs to sub-fs duration) with high momentum and arrival time stability for the purpose of applications related to accelerator R&D, e.g. development of advanced and compact diagnostics and accelerating structures, test of new accelerator components, medical applications studies, machine learning, etc. During its commissioning phase, the bunch duration characterization of the electron bunches generated at ARES has been performed with an RF-phasing technique relying on momentum spectra measurements, using only common accelerator elements (RF accelerating structures and magnetic spectrometers). The sensitivity of the method allowed highlighting different response times for Mo and Cs2Te cathodes. The measured electron bunch duration in a wide range of machine parameters shows excellent agreement overall with the simulation predictions, thus demonstrating a very good understanding of the ARES operation on the bunch duration aspect. The importance of a precise in-situ experimental determination of the phase velocity of the first travelling wave accelerating structure after the electron source, for which we propose a simple new beam-based method precise down to sub-permille variation respective to the speed of light in vacuum, is emphasized for this purpose. A minimum bunch duration of 20 fs rms, resolution-limited by the space charge forces, is reported. This is, to the best of our knowledge, around 4 times shorter than what has been previously experimentally demonstrated based on RF-phasing techniques with a single RF structure. The present study constitutes a strong basis for future time characterization down to the sub-fs level at ARES, using dedicated X-band transverse deflecting structures.Comment: 17 pages, 11 figures. To be submitted to Physical Review Accelerators and Beam