Van der Waals loops and the melting transition in two dimensions

Evidence for the existence of van der Waals loops in pressure p versus volume v plots has for some time supported the belief that melting in two dimensions is a first order phase transition. We report rather accurate equilibrium p(v) curves for systems of hard disks obtained from long Monte Carlo simulations. These curves, obtained in the constant volume ensemble, using periodic boundary conditions, exhibit well defined van der Waals loops. We illustrate their existence for finite systems that are known to undergo a continuous transition in the thermodynamic limit. To this end, we obtain magnetization m versus applied field curves from Monte Carlo simulations of the 2D Ising model, in the constant m ensemble, at the critical point. Whether van der Waals loops for disk systems behave in the thermodynamic limit as they do for the 2D Ising model at the critical point cannot be ruled out. Thus, the often made claim that melting in 2D is a first order phase transition, based on the evidence that van der Waals loops exist, is not sound.Comment: 10 pages, 6 Postscript figures (submitted to Phys.Rev.E). For related work, see http://pipe.unizar.es/~jf

Simulations of Two-Dimensional Melting on the Surface of a Sphere

We have simulated a system of classical particles confined on the surface of a sphere interacting with a repulsive $r^{-12}$ potential. The same system simulated on a plane with periodic boundary conditions has van der Waals loops in pressure-density plots which are usually interpreted as evidence for a first order melting transition, but on the sphere such loops are absent. We also investigated the structure factor and from the width of the first peak as a function of density we can show that the growth of the correlation length is consistent with KTHNY theory. This suggests that simulations of two dimensional melting phenomena are best performed on the surface of a sphere.Comment: 4 eps figure

Liquid-Solid Phase Transition of the System with Two particles in a Rectangular Box

We study the statistical properties of two hard spheres in a two dimensional rectangular box. In this system, the relation like Van der Waals equation loop is obtained between the width of the box and the pressure working on side walls. The auto-correlation function of each particle's position is calculated numerically. By this calculation near the critical width, the time at which the correlation become zero gets longer according to the increase of the height of the box. Moreover, fast and slow relaxation processes like $\alpha$ and $\beta$ relaxations observed in supper cooled liquid are observed when the height of the box is sufficiently large. These relaxation processes are discussed with the probability distribution of relative position of two particles.Comment: 6 figure

Elastic moduli, dislocation core energy and melting of hard disks in two dimensions

Elastic moduli and dislocation core energy of the triangular solid of hard disks of diameter $\sigma$ are obtained in the limit of vanishing dislocation- antidislocation pair density, from Monte Carlo simulations which incorporates a constraint, namely that all moves altering the local connectivity away from that of the ideal triangular lattice are rejected. In this limit, we show that the solid is stable against all other fluctuations at least upto densities as low as $\rho \sigma^2 = 0.88$. Our system does not show any phase transition so diverging correlation lengths leading to finite size effects and slow relaxations do not exist. The dislocation pair formation probability is estimated from the fraction of moves rejected due to the constraint which yields, in turn, the core energy E_c and the (bare) dislocation fugacity y. Using these quantities, we check the relative validity of first order and Kosterlitz-Thouless-Halperin-Nelson-Young (KTHNY) melting scenarios and obtain numerical estimates of the typical expected transition densities and pressures. We conclude that a KTHNY transition from the solid to a hexatic phase preempts the solid to liquid first order transition in this system albeit by a very small margin, easily masked by crossover effects in unconstrained ``brute-force'' simulations with small number of particles.Comment: 17 pages, 8 figure

Arecibo pulsar survey using ALPHA: Probing radio pulsar intermittency and transients

We present radio transient search algorithms, results, and statistics from the ongoing Arecibo Pulsar ALFA (PALFA) survey of the Galactic plane. We have discovered seven objects through a search for isolated dispersed pulses. All of these objects are Galactic and have measured periods between 0.4 and 4.7 s. One of the new discoveries has a duty cycle of 0.01%, smaller than that of any other radio pulsar. We discuss the impact of selection effects on the detectability and classification of intermittent sources, and compare the efficiencies of periodicity and single-pulse (SP) searches for various pulsar classes. For some cases we find that the apparent intermittency is likely to be caused by off-axis detection or a short time window that selects only a few bright pulses and favors detection with our SP algorithm. In other cases, the intermittency appears to be intrinsic to the source. No transients were found with DMs large enough to require that they originate from sources outside our Galaxy. Accounting for the on-axis gain of the ALFA system, as well as the low gain but large solid-angle coverage of far-out sidelobes, we use the results of the survey so far to place limits on the amplitudes and event rates of transients of arbitrary origin. © 2009. The American Astronomical Society. All rights reserved.

Enhanced stability of the square lattice of a classical bilayer Wigner crystal

The stability and melting transition of a single layer and a bilayer crystal consisting of charged particles interacting through a Coulomb or a screened Coulomb potential is studied using the Monte-Carlo technique. A new melting criterion is formulated which we show to be universal for bilayer as well as for single layer crystals in the case of (screened) Coulomb, Lennard--Jones and 1/r^{12} repulsive inter-particle interactions. The melting temperature for the five different lattice structures of the bilayer Wigner crystal is obtained, and a phase diagram is constructed as a function of the interlayer distance. We found the surprising result that the square lattice has a substantial larger melting temperature as compared to the other lattice structures. This is a consequence of the specific topology of the defects which are created with increasing temperature and which have a larger energy as compared to the defects in e.g. a hexagonal lattice.Comment: Accepted for publication in Physical Review