Critical exponents and unusual properties of the broken phase in the 3d-RP2 antiferromagnetic model
We present the results of a Monte Carlo simulation of the antiferromagnetic RP^2 model in three dimensions. 
With finite-size scaling techniques we accurately measure the critical exponents and compare them with those of 
O(N) models. We are able to parameterize the corrections-to-scaling. The symmetry properties of the broken 
phase are also studied.
Mean field and computer simulation study of a nematogenic lattice model including three-body interactions
Abstract By now, nematogenic lattice models have been extensively studied in the literature dealing with liquid crystals; they usually involve cylindrically symmetric (uniaxial) particles and pairwise additive interaction potentials. On the other hand, quantum mechanical perturbation theory shows that interatomic or intermolecular potentials are only approximately pairwise additive; the pairwise additivity approximation (PAA) of molecular interactions has been extensively and systematically used in the statistical mechanics of condensed matter, and has proven rather successful. As an attempt to move beyond the PAA in the simulation of mesogenic systems, we have considered here a nematogenic lattice model consisting of uniaxial particles, whose centres of mass are associated with a simple-cubic lattice, and whose interaction potentials consist of a pairwise additive dispersion (Nehring–Saupe) term, restricted to nearest neighbours (and already studied in the literature), plus a short-range triplet-additive one, i.e., the Kielich–Stogryn generalization of the Axilrod–Teller–Muto formula for three atoms. The model has been studied by Mean Field theory and Monte Carlo simulation; the three-body term was found to produce a recognizable quantitative effect on the nematic ordering transition.
Solving φ1,24 field theory with Monte Carlo
We study lattice g0φ4 field theory for all g0 and fixed renormalized mass M in one and two dimensions using Monte Carlo techniques. We calculate the dimensionless renormalized coupling constant gR = gRM4−d, where d is the dimension of space—time, at fixed small values of the lattice spacing a for various g0 and lattice sizes. Our results are in quantitative agreement with the analyses of high temperature and strong coupling series which rely on extrapolation from large to small lattice spacing.