Equilibrium fluid-solid coexistence of hard spheres

We present a tethered Monte Carlo simulation of the crystallization of hard spheres. Our method boosts the traditional umbrella sampling to the point of making practical the study of constrained Gibbs' free energies depending on several crystalline order parameters. We obtain high-accuracy estimates of the fluid-crystal coexistence pressure for up to 2916 particles (enough to accommodate fluid-solid interfaces). We are able to extrapolate to infinite volume the coexistence pressure [p(co)=11.5727(10)k(B)T/σ(3)] and the interfacial free energy [γ({100})=0.636(11)k(B)T/σ(2)].     

Lattice-spin mechanism in colossal magnetoresistive manganites

We present a single-orbital double-exchange model, coupled with cooperative phonons (the so called breathing modes of the oxygen octahedra in manganites). The model is studied with Monte Carlo simulations. For a finite range of doping and coupling constants, a first-order metal-insulator phase transition is found, which coincides with the paramagnetic-ferromagnetic phase transition. The insulating state is due to the self-trapping of every carrier within an oxygen octahedron distortion.     

A measure of conductivity for lattice fermions at finite density

We study the linear response to an external electric field of a system of fermions in a lattice at zero temperature. This allows to measure numerically the Euclidean conductivity which turns out to be compatible with an analytical calculation for free fermions. The numerical method is generalizable to systems with dynamical interactions where no analytical approach is possible.     

Review of Recent Developments in the Random-Field Ising Model

A lot of progress has been made recently in our understanding of the random-field Ising model thanks to large-scale numerical simulations. In particular, it has been shown that, contrary to previous statements: the critical exponents for different probability distributions of the random fields and for diluted antiferromagnets in a field are the same. Therefore, critical universality, which is a perturbative renormalization-group prediction, holds beyond the validity regime of perturbation theory. Most notably, dimensional reduction is restored at five dimensions, i.e., the exponents of the random-field Ising model at five dimensions and those of the pure Ising ferromagnet at three dimensions are the same.     

Vibrational spectrum of topologically disordered systems

The topological nature of the disorder of glasses and supercooled liquids strongly affects their high-frequency dynamics. In order to understand its main features, we analytically studied a simple topologically disordered model, where the particles oscillate around randomly distributed centers, interacting through a generic pair potential. We present results of a resummation of the perturbative expansion in the inverse particle density for the dynamic structure factor and density of states. This gives accurate results for the range of densities found in real systems.     

Measures of critical exponents in the four-dimensional site percolation

Using finite-size scaling methods we measure the thermal and magnetic exponents of the site percolation in four dimensions, obtaining a value for the anomalous dimension very different from the results found in the literature. We also obtain the leading corrections-to-scaling exponent and, with great accuracy, the critical density.     

First-order transition in a three-dimensional disordered system

We present the first detailed numerical study in three dimensions of a first-order phase transition that remains first order in the presence of quenched disorder (specifically, the ferromagnetic-paramagnetic transition of the site-diluted four states Potts model). A tricritical point, which lies surprisingly near the pure-system limit and is studied by means of finite-size scaling, separates the first-order and second-order parts of the critical line. This investigation has been made possible by a new definition of the disorder average that avoids the diverging-variance probability distributions that plague the standard approach. Entropy, rather than free energy, is the basic object in this approach that exploits a recently introduced microcanonical Monte Carlo method.     

Phase diagram of a polydisperse soft-spheres model for liquids and colloids

The phase diagram of soft spheres with size dispersion is studied by means of an optimized Monte Carlo algorithm which allows us to equilibrate below the kinetic glass transition for all size distributions. The system ubiquitously undergoes a first-order freezing transition. While for a small size dispersion the frozen phase has a crystalline structure, large density inhomogeneities appear in the highly disperse systems. Studying the interplay between the equilibrium phase diagram and the kinetic glass transition, we argue that the experimentally found terminal polydispersity of colloids is a purely kinetic phenomenon.